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1.  INTRODUCTION 

The  frictional  resistance  of  a  surface  and  its  correlation  to  the 
surface  geometry  is  of  interest  to  workers  in  a  wide  range  of  dis¬ 
ciplines.  Manufacturers  of  extruded  materials,  for  example,  are 
interested  in  the  flow  of  the  material  through  the  extrusion  nozzle. 

The  shape  and  surface  geometry  of  the  nozzle  will  affect  the  subsequent 
mechanical  properties  of  the  extruded  product. 

The  flow  around  the  relatively  complex  shape  of  a  ship's  hull  is 
dependent  on  both  the  form  of  the  hull  and  on  its  surface  finish. 
Ideally,  if  surfaces  could  be  made  hydro-  or  aero-dynamically  smooth 
then  the  naval  architect's  job  of  designing  a  ship  form  with  minimum 
pressure  drag  and  minimum  shear  drag  would  provide  the  optimum  solution 
in  terms  of  power  required  to  propel  the  vessel  at  set  conditions. 

The  process  of  transforming  the  architect's  idealised  ship  form 
to  the  manufactured  product  superimposes  a  surface  finish  onto  the  ship 
form.  This  practical  surface  finish,  even  when  the  ship  is  brand-new, 
fails  to  satisfy  the  requirement  of  hydrodynamically  smooth  condition. 

The  study  of  the  relationship  between  the  surface  finish  and  the 
overall  performance  of  the  ship  hull  (see  Appendix  1)  has  shown  there 
to  be  potentially  very  rewarding  reductions  in  power  consumption  with 
the  introduction  of  quality  control  on  surface  finish. 

This  investigation  sets  out  to  correlate  specific  parameters  of 
surface  finish  to  significant  increases  in  the  frictional  drag  for  the 
case  of  open-channel  flov;  over  nominally  flat  surfaces.  The  work  is 
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specifically  related  to  the  problem  of  the  frictional  drag  of  ship  hull 
surfaces  and  the  surfaces  investigated  include  a  4.8  metre  long  replica 
surface  of  part  of  a  ship  hull. 

In  an  investigation  of  the  viscous  shear  and  pressure  forces 
resisting  the  motion  of  a  ship  hull  it  is  important  to  know  what  the 
enveloping  fluid  ’sees'  of  the  hull  surface.  On  a  typical  in-service 
ship  hull  the  external  surface  is  made  up  of: 

1.  The  roughness  of  the  component  sections  of  the  hull  surface  - 
pai:-  ted  steel  plates  whose  surface  is  made  up  of  wavelengths  from 
microns  to  metres.  The  short  wavelengths  are  made  up  from:  (a)  the 
substrate  surface  finish;  (b)  the  method  of  paint  application;  (c)  the 
physical  properties  of  the  paint;  and  (d)  the  paint  composition.  The 
longer  wavelengths  are  the  waviness  and  errors  of  form  of  the  hull 
plates. 

2.  A  spectrum  of  surface  wavelengths  produced  from  the  fabri¬ 
cation  techniques  used  in  constructing  the  hull  surface.  The  short 
wavelengths  are  made  up  of  joints,  bulges  and  dents.  The  longer 
surface  wavelengths,  of  up  to  hundreds  of  metres,  make  up  the  overall 
shape  or  curvature  of  the  hull  and  the  errors  of  form  of  the  hull 
shape . 

3.  The  smaller  scale  roughness,  wavelengths  from  microns  to 
several  metres,  can  be  substantially  altered  by:  (a)  chemical  ana 
mechanical  damage;  (b)  the  attachment  of  rigid  marine  growths  such  as 
barnacles;  (c)  the  attachment  of  compliant  growths  of  slimes  and  weeds. 
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Thus,  the  surface  of  a  typical  ship  hull  comprises  a  large  number 
of  surface  features.  Each  of  the  surface  features  on  the  hull  is  bound 
to  effect  both  a  viscous  shear  component  and  a  pressure  (form)  component 
to  the  frictional  ’-esistance  of  the  surface. 

The  positioning  of  a  surface  feature  on  the  hull  will  determine 
its  contribution  to  the  total  resistance  of  the  hull  surface,  since 
the  outline  shape  of  the  hull  will  produce  a  variety  of  flow  regimes 
around  the  hull.  Hence  for  a  full  analysis  of  the  frictional  resis¬ 
tance  of  a  ship  hull  surface  a  large  selection  of  surface  textures  and 
flow  regimes  must  be  examined  and  the  interaction  between  the  different 
flow  regimes  and  surfaces  investigated. 

The  momentum  balance  of  the  flow  provides  the  fundamental  equa¬ 
tions  for  the  determination  of  the  flow  parameters  whilst  a  statistical 
analysis  of  the  surface  geometry  is  used  to  characterise  the  surface. 

The  existence  of  a  logarithmic  relationship  for  a  region  of  the  flow 
near  to  a  solid  boundary  can  be  used  as  an  alternative  means  of 
determining  the  flow  parameters.  Whilst  the  logarithmic  'lav?*  (see 
Appendix  2) ,  in  its  simplest  form,  is  valid  for  two  dimensional  flow 
over  hy drodyn ami cally  smooth  surfaces,  leading  researchers  have  modi¬ 
fied  it  to  include  surface  roughness  effects.  Hence  a  number  of 
empirical  equations  have  been  formed  with  various  roughness  function 
terms  in  an  attempt  to  correlate  surface  frictional  resistance  with 
specific  roughness  parameters  such  as  mean  apparent  amplitude  and  root- 
mean-square  roughness.  Although  the  basic  logarithmic  relationship 
has  been  shovyn  to  exist  using  the  principles  of  dynamical  similarity 
and  other  criteria  the  various  off-shoot  correlations  have  no  real 
theoretical  basis  and  therefore  their  use  in  this  investigation  is 
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2-  TECHNIQUES  FOR  THE  EVALUATION  OF  THE  FRICTIONAL  DRAG  OF 
SURFACES 


2*1  The  choice  of  flow  system 

An  appropriate  method  of  investigating  the  relation  between  the 
surface  frictional  resistance  of  a  ship  hull  surface  and  its  geomatry 
would  be  to  take  flow  measurements  around  the  hull  of  a  ship  whilst  she 
is  in  service  and  then  correlate  the  resulting  flow  characteristics 
with  the  corresponding  surface  parameters  measured  whilst  the  vessel  is 
in  dry  dock. 

The  practicalities  of  obtaining  surface  friction  values  around  a 
ship  hull  present  economic  and  instrumentation  problems. 

A  number  of  v;orkers  have  carried  out  ship  trials  from  which  the 
contribution  the  hull  roughness  Eakes  to  the  hull's  total  drag  has  been 
determined  from  an  analysis  of  the  ship's  propulsion  power  requirements. 
There  is  the  inherent  difficulty  of  separating  the  form  (pressure- 
based)  drag  from  the  roughness  (viscous-based)  drag.  This  stems  from 
the  absence  of  flow  measurements  on  hvdrodynamically  smooth  hulls 
(Granville  1956) . 

Kore  economically  viable  drag  measurements  on  scaled-down  models 
of  ships  suffer  from  the  difficulties  of  the  subsequent  scaling-up  of 
the  ship  model  results  to  full-size.  This  involves  checking  the 
ship-model  correlation  using  full-scale  ship  trials  (see  Appendix  )). 

Apart  from  the  methods  cf  evaluating  surface  frictional  resistance 
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The  removal  of  the  long  surface  wavelengths  when  going  from  the 
curved  body  of  a  ship  to  the  parallel  'lov?  situation  on  a  flat  plate  or 
in  a  pipe  assumes  that  these  relatively  long  surface  wavelengths,  or 
features,  make  no  contribution  to  the  frictional  resistance  of  the 
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surface.  We  are  therefore  faced  with  the  problem  of  deciding  when  a 
surface  wavelength  is  such  that  it  no  longer  contributes  to  the  fric¬ 
tional  resistance  of  the  surface. 

As  yet  there  is  no  theory  to  call  upon  to  distinguish  this  limiting 
value  of  surface  wavelength  and  therefore  it  would  be  presumptious  to 
neglect  surface  features  above  or  below  any  arbitrary  value.  Also,  in 
attempting  to  make  a  comparison  between  hull  frictional  resistance  and 
piate/pipes  frictional  resistance  we  are  dealing  with  three  fundamen¬ 
tally  different  flow  patterns  each  of  which  will  be  affected  differently 
by  the  surfaces  under  investigation. 

In  the  pipe,  the  flow  is  constrained  on  all  sides  and  therefore 
the  shear  layer  is  confined  to  the  diameter  of  the  pipe,  in  comparison 
to  flew  around  a  ship  where  there  are  no  constraining  boundaries  and 
the  boundary  layer  is  free  to  develop  around  the  greater  part  of  the 
submerged  hull. 

In  plate  flow  if  we  assume  bne  plate  to  be  infinitely  wide,  the 
flow  will  bear  a  closer  resemblance  to  the  flow  around  a  ship,  but  the 
flow  will  be  affected  by  the  water-air  interface.  If  the  plate  com¬ 
prises  the  floor  of  a  channel,  then  there  will  be  three-dimensional 
effects  due  to  the  channel  sides.  Therefore,  in  order  to  reduce  the 
number  of  differences  in  the  comparison  of  the  ship  and  laboratory 
flow  regimes  to  a  minimum,  a  flow  channel  is  required  which  is  as  wide 
and  as  deep  as  possible  and  even  then  account  must  be  taken  of  the 
inherent  modifications  to  the  flow  (Grass  1871} . 

Using  laboratory  tests  in  pipes  and  on  plates  it  is  only  possible 
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The  flow  system  used  in  the  present  ii.  estigation  consists  of  a 
parallel  glass-sided,  variable  slope  water  channel.  The  channel  sides 
are  glass  panels,  1.5  metres  long,  carefully  butted  together  and  the 
joints  filled.  The  channel’s  variable  slope  facility  along  with  an 
adjustable  outlet  weir  enable  a  wide  variation  in  the  flow  configura¬ 
tion  within  the  open  top  channel.  The  flow  conditions  were,  in  fact, 
restricted  to  sub-critical  (i.e.  tranquil)  mode  to  obviate  the 
phenomena  of  shock  waves  forming  at  the  sensing  probes  and  oblique 
shocks  generated  from  the  imperfections  of  the  channel  sides.  The 
channel  is  approximately  20  m  long  and  30.5  cm  wide  and  can  accommodate 
a  plate  or  replica  of  a  ship  hull  30,5  cm  wide  up  to  20  m  long. 

A  deturbulating  tank  with  a  streamlined  converging  section  is  bolted 
to  the  channel  at  inlet  and  is  used  to  produce  reasonably  uniform  flow 
at  inlet,  free  of  swirl  and  flow  stratification  -  ideally  one-dimen¬ 
sional  potential  flow. 

The  floor  of  the  channel  provides  the  base  for  the  working 
surface,  and  is  made  up  of  box  sections  30.5  cm  wide  by  3.05  m  long 
which  are  securely  bolted  together  and  supported  at  three  points,  - 
on  rollers  at  two  jacking  points  and  at  one  hinge  point.  A  schematic 
diagram  of  the  flew  system  is  shown  in  Figure  2-1. 


to  determine  it3  suitability  as  applied  to  a  study  of  the  frictional 
resistance  of  ship  hull  surfaces.  It  is  realised  that  the  flow  around 
the  hull  of  a  ship  is  very  different  to  the  confined  flow  of  a  water 
channel,  and  a  comparison  is  only  valid  if  account  is  taken  of  these 
differences  and  any  results  from  the  flow  measurements  adjusted 
accordingly. 


2 . 3  Components  of  the  flow  system  | 

Two  pumps  are  available  to  provide  the  fluid  motion  through  the  4 

open  channel.  The  first  of  these,  pump  A,  is  capable  of  producing  4 

'  3 
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channel  Reynolds  numbers,  Rc  of  up  to  a  maximum  of  3  x  10  ,  where  Rc  =  ;  1 

!  | 

j  ^ 

Uex/v.  Comparing  this  Reynolds  number  with  typical  ships  Reynolds  :  1 

I 

number,  a  modem  sea-going  vessel  such  as  a  400  metre  long  tanker  3 

« 
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cruising  at  15  knots  (7.7  m/sec)  develops  a  Reynolds  number  of  3  x  10  .  4 

I 

The  ships  Reynolds  number  is  given  by  Rs  =  UsLs/v  where  Us  =  speed  of  1 

I 

I 

the  ship  and  Ls  =  the  length  of  the  ship.  j 


A  second,  larger  capacity  pump  B,  loaned  by  BSRA  has  been 
installed  into  the  flow  system  and  this  has  been  used  to  develop 
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Reynolds  numbers  of  up  to  10  . 


The  maximum  usable  value  of  R^.  from  pump  B  is  somewhat  less  than 
this,  in  the  region  of  7  x  10®.  This  is  because  at  the  higher 
Reynolds  numbers  the  channel  flow  becomes  supercritical  which  is 
undesirable  for  the  present  flow  instrumentation .  Apart  from  the 
limiting  value  of  Reynolds  numbc's  due  to  the  onset  of  critical  flow 
there  is  also  a  limiting  depth  of  flow  below  which  the  flow  effects 


due  to  the  water-air  interface  become  significant. 


The  limiting  depth 
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of  flow  in  the  channel  was  found  to  be  in  the  region  of  6  inches  (15  cm). 

The  two  pumps  are  operated  separately,  the  inactive  pump  is  iso¬ 
lated  from  the  flow  system  using  valves  I  and  II  (Figure  2-1).  Pump  B, 
the  main  pump,  is  a  roto-dynamic  pump  of  the  radial  flow  type  with  a 
36  cm  (internal  diameter)  inlet  and  an  18  cm  (internal  diameter)  outlet. 
In  the  present  system  the  pump  outlet  is  connected  into  the  15  cm 
(internal  diameter)  supply  line  to  the  flov;  channel.  A  Dali  flow  tube 
is  positioned  in  the  supply  line  with  its  pressure  tappings  connected 
to  a  mercury  underwater  manometer.  From  the  supply  line  the  water 
enters  the  flow  channel  via  the  deturbulating  tank  and  a  series  of 
honeycombs . 

Pump  B  operates  under  conditions  of  a  negative  static  suction  head, 
drawing  water  from  a  below  ground  level  pit,  2  metres  deep  with  a  water 
capacity  of  7.5  m3  (1.7  thousand  gallons).  The  pump  when  commissioned 
could  deliver  4.0  m3/rain.  at  960  rpm,  but  the  pump  has  subsequently 
been  derated  by  reducing  the  size  of  the  impeller,  to  approximately 
3.5  m3/min.  at  960  rev/min. 

The  smaller  capacity  pump  A  is  an  axial  flow  extraction  pump  with 
the  pump  head  submerged  in  the  water  pit.  The  pump  is  capable  of 
delivering  1.2  m3/min.  at  2900  rev/min. 

For  both  pumps  it  was  found  necessary  to  introduce  two  honeycomb 
sheets  into  the  deturbulating  tank  in  order  to  remove  large  scale 
eddies  and  swirl  which  was  otherwise  entering  the  channel  and  pro¬ 
ducing  non-uniform  flow  at  the  channel  inlet.  The  swirl  probably 
originates  from  two  ninety  degree  bends  in  the  pipe  line  preceding  the 
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deturbulating  tank.  One  of  the  honeycomb  flow  straignteners  is  posi¬ 
tioned  horizontally  over  the  deturbulating  tank  inlet,  and  the  other 
vertically  across  the  channel  inlet.  Such  an  arrangement  produced  the 
desirable  flow  conditions  at  the  inlet  to  the  channel. 


A  continuous  cycle  filtration  unit  is  incorporated  in  the  system 
in  order  to  keep  the  water  at  a  high  level  of  purity  required  to  operate 
the  flow  instrumentation.  For  this,  a  portable  extraction  pump  draws 
water  from  the  sump.  The  water  is  then  allowed  to  cascade  by  gravity 
through  a  series  of  different  mesh  size  filter  gauzes  back  into  the 
sump.  The  gauzes  are  periodically  removed  and  cleaned. 


2.4  Techniques  for  measuring  Cf 


There  are  several  techniques  available  for  the  measurement  of  the 
local  coefficient  of  friction  of  a  surface.  Only  a  few  of  these  are 
suitable  for  the  particular  experimental  set-up  used.  These  can  be 
divided  into  four  main  types: 

1.  Pressure  measurements. 

2.  Force  measurements. 

3.  Heat  transfer  measurements. 

4.  Optical  methods. 


2.4.1  Pressure  measurements 


Devices  such  as  Preston  or  Stanton  tubes  can  be  used  to  measure 
the -surface  friction  directly.  These  are  probably  the  most  widely 
used  of  the  manometric  methods.  Both  are  based  or  the  universality  of 
the  velocity  profile  in  the  inner  wall  region.  Preston's  (1954) 
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method  of  determining  Cf  was  to  use  a  calibrated  pitot  tube  placed  in 
contact  with  the  wall.  The  dynamical  pressure  p  exerted  on  the  tip  of 
the  tube  depends  on  P/V/T^,  the  outer  radius  r  and  the  inner  radius  ri 
of  the  tube.  The  pressure  measured  by  the  total-head  tube  in  conjunc¬ 
tion  with  the  surface  static  pressure  measured  at  the  same  station  along 
the  plate  was  calibrated  by  Preston  in  terms  of  the  local  surface-shea.' 
stress. 

There  are  several  other  types  of  pressure  devices  which  can  be 
used  directly  or  indirectly  to  determine  the  surface  friction.  Devices 
such  as  the  Preston  tube  can  only  be  used  on  relatively  smooth  flat 
surfaces  and  provide  an  obstruction  to  the  flow  with  subsequent  alter¬ 
ations  in  the  flow  pattern  around  the  area  of  measurement.  Pitot- 
static  tubes  can  be  used  to  measure  mean  velocity  profiles  in  the 
boundary  layer  which  can  then  be  used  in  the  computation  of  Cf.  The 
velocity  measurements  are  absolute  but  they  cannot  be  made  near  to  the 
work  surface  since  the  resulting  gap  between  the  tube  and  the  surface 
will  act  as  a  constriction  to  the  flow,  producing  secondary  flow 
effects  which  may  be  measured  by  the  tube  itself. 

In  a  fully  developed  turbulent  pipe  flow  the  friction  loss  or 

pressure  drop,  simply  related  to  the  wall  shear  stress  x  t  can  be 

a) 

functionally  related  to  the  average  velocity  U  of  the  pipe  flow,  pipe 
diameter  D,  fluid  properties  of  density  p  and  kinematic  viscosity  v, 

and  a  geometrical  description  of  the  rough  wall  k,  kj,  k^, . 

Musker  et  al  (1976)  used  measurements  from  pipe  flow  to  predict  local 
surface  friction  coefficients  for  a  ship  hull  surface,  assumed  equiva¬ 
lent  to  a  plane  surface.  The  pipe  flow  average  velocity  is  chosen  so 
that  dynamic  similarity  between  the  flow  in  the  inner  regions  of  the 
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pipe  flow  and  that  on  the  ship  hull  can  be  assumed.  The  law  of  the 
wall  combined  with  the  law  of  the  wake  (Coles  1956)  can  be  used  to 
represent  both  the  boundary  layer  flov;  on  a  flat  plate  without  pressure 
gradient  and  the  fully  developed  turbulent  pipe  flov;  for  smooth  sur¬ 
faces.  The  effect  of  surface  roughness  on  the  mean  velocity  distri¬ 
bution  is  assumed  to  be  given  by 


where  C  is  the  law  of  the  wall  constant  which  depends  on  kuT/v, 
k  is  the  ‘typical’  roughness  height. 

Alternatively,  the  mean  velocity  distribution  can  be  written  as 


where  X  is  a  roughness  function  primarily  dependent  on  the  roughness 
Reynolds  number  kuT/v. 


From  measurements  of  the  static  pressure  drop  and  the  mean  vel¬ 
ocity  profiles  in  the  pipe  flov;  the  roughness  function  x  •  given  by 


X 


i°ge 
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‘-e~ClRk  +  C2R  J 


(Musker  at  al  1976) 


where  B  =  law  of  the  wall  constant,  R^  =  kuT/v,  and  and  C2 
are  empirical  constants  determined  from  the  roughness  function, 
is  calculated  for  a  given  surface.  The  'working  section  of  the  pipe  is 
lined  with  twelve  identical  61  cm  long  cylindrical  replicas  of  the 


surface  under  investigation. 


Thus,  by  using  dynamical  similarity  which  requires  that  the 
roughness  Reynolds  numbers  for  the  pipe  and  flat  plate  are  identical, 
the  full  scale  flat  plate  surface  friction  coefficients  can  be  cal¬ 
culated  knowing  the  roughness  function  for  a  given  surface  geometry. 


2.4.2  Force  measurements 


The  local  surface  friction  can  be  measured  directly  on  a  specimen 
surface  mounted  movably  in  a  sector  of  the  principal  surface.  The 
arrangement  is  given  the  name  floating-element  sur  oe  friction  balance 
(Smith  &  Walker  1957,  Reda  1974,  Driftmyer  1975  and  Karlsson  1978) . 


The  floating  element  is  used  by  Karlsson  (1978)  to  measure  local 
surface  friction  values  in  a  flat  plate  turbulent  boundary  layer.  The 
measurements  are  part  of  an  investigation  into  the  effect  of  irregular 
ship  hull  roughness  on  the  frictional  resistance  of  full  scale  ships. 


The  floating  element  is  located  on  the  centre  line  of  a  wind 
tunnel  wall  covered  with  a  positive  replica  of  a  rough  ship  hull  sur¬ 
face.  the  floating  surface  is  located  flush  with  the  surrounding 
surface  using  a  dial  gauge.  Pressure  tappings  are  located  on  the  up¬ 
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stream  and  downstream  edges  of  the  element  to  allow  the  measurement  of 
the  pressure  in  the  gaps  ('v  0.3  mm  wide). 


The  flow  exerts  a  force  F  on  the  element  surface  and  the  element 
is  displaced  in  the  flow  direction.  F  is  given  by 


F  ~  VA 


where  A  is  the  area  of  the  floating  surface. 


A  micrometer  screw  is  used  to  restore  the  element  to  its  unloaded 


position,  which  is  indicated  by  a  differential  transformer.  If  there  is 
no  pressure  difference  around  the  edge  of  the  element,  the  restoring 
force  Fr  will  equal  the  force  exerted  by  the  flow.  Pressure  differences 
in  the  gap  around  the  element  will  produce  an  additional  force  Fp 
given  by 


Fp  =  Apbh 


where  b  is  the  width  of  the  floating  element,  h  is  the  height  of  the 
floating  element  edge.  Thus,  resolving  the  forces  in  the  direction  of 
the  fluid  motion,  gives 


F  =  FR  +  Fp  (Figure  2-2) 


The  local  surface  friction  coefficient  is  obtained  from  the 


equation 


Cf=^ 


The  reference  velocity  Ue  was  recorded  for  all  the  measurements 
using  a  pitot-static  tube  of  Prandtl's  design  and  an  accurate  single¬ 
limb  micromanometer  (Preston  1954).  Thus,  readings  taken  from  the 
floating  element  together  with  the  velocity  U£  and  appropriate  values 
of  atmospheric  pressure,  wind  tunnel  temperature  etc.,  were  used  to 
calculate  the  local  surface  friction  coefficient  as  a  function  of  unit 
Reynolds  number  UQ/v 


Roughness  functions  (Hama  1954)  were  also  determined  for  each  of 
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the  replicated  surfaces  investigated. 

Although  the  boundary- layer  flow  on  a  rotating  disc  is  three- 

dimensional  in  constrast  to  the  two-dimensional  flow  on  a  flat  plate 

or  axi-symmetric  through  a  pipe,  the  rotating  disc  provides  a  convenient 

means  of  measuring  surface  friction  effects  particularly  at  high  speeds. 

The  torque  or  moment  M  of  one  side  of  a  rotating  disc  in  an  unbounded 

fluid  medium  can  be  assumed  to  depend  upon  the  angular  velocity  w,  the 

radius  R,  fluid  properties  of  density  p  and  kinematic  viscosity  v,  and 

the  geometry  of  the  disc's  surface  k,  k. ,  k~  . 

X  ^ 

By  dimensional  analysis  the  moment  coefficient  Cm  is  a  function  of 
two  dimensionless  groups,  a  Reynolds  number  Re  and  a  relative  roughness 

'V. 

R/k  for  a  particular  surface  geometry  R.  Hence 
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For  a  smooth  surface  AB  =  0  and  (A3) '  =  0. 
Hence, 


=  A  /T  log  R  Jc~  +  /T  B,  „ - - — 

^  /ST  e  e  m  / btT  1,0  /io~ 


~  A  [Jl  1o%  A" 

/  8it  18 


Thus,  correlating  plots  of  l/*^-  against  log  Re  &  for  smooth 

m  m 

and  rough  surfaces,  the  local  friction  at  the  edge  of  the  disc  is  given 
by 


{l  -  [2A  +  {AB)e' 


A  similar  analysis  can  be  used  to  relate  the  torque  or  a  rotating 
drum  in  an  unbounded  fluid  to  its  surface  geometry.  This  technique  is 
presently  used  by  workers  at  Newcastle  University  to  investigate  the 
frictional  properties  of  ship  hull  surfaces. 


2.4.3  Heat  transfer  measurements 


Instruments  based  on  a  heat  transfer  principle  such  as  the  hot¬ 
wire  anemometer  have  a  small  size  of  the  sensing  element  producing 
negligible  disturbances  to  the  flow,  good  response  to  high-frequency 
fluctuations  and  suitability  for  electron' c  instrumentation.  The 
instrument'  can  be  used  within  a  few  thousandths  of  an  inch  of  a  solid 
surface.  The  detecting  element  consisting  of  a  very  fine  short  metal 
wire  or  film,  only  a  few  microns  in  diameter.  The  element  is  heated 
by  ar  electrical  current  and  when  placed  in  a  flowing  fluid  is  cooled 
by  the  fluid,  causing  the  temperature  to  drop  and,  consequently,  the 
electrical  resistance  of  the  wire  to  diminish.  Hence  the  anemometer 
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can  be  used  to  record  fluctuating  velocities  of  very  high  frequency 
and  in  the  case  of  evaluating  Cf  can  be  used  to  measure  the  moan  velo¬ 
city  distribution  down  to  very  near  to  the  solid  surface.  The  main 
disadvantage  of  the  anemometer  is  that  it  is  comparative  rather  than 
absolute  and  thus  its  accuracy  depends  on  the  method  of  calibration. 
Ideally  the  detecting  element  should  be  calibrated  in  the  same  flow 
situation  as  that  in  which  it  is  to  be  used  so  as  to  reduce  the  number 
of  variables  affecting  the  calibration.  The  heat- transfer  character¬ 
istics  of  the  sensing  element  can  change  if  it  becomes  contaminated 
with  dirt  or  if  in  the  case  of  water  flow,  bubbles  fora  on  it.  Hence 
the  instrument  can  only  be  operated  in  a  'clean'  environment.  When  the 
sensor  is  placed  in  close  proximity  to  a  solid  boundary  having  a  diff¬ 
erent  thermal  conductivity  than  the  surrounding  fluid,  the  heat  loss 
from  the  sensing  element  will  vary  (Wills  1962) . 


The  standard  anemometer  unit  has  one  sensing  element  which  can 
measure  velocity  fluctuations  in  two  dimensions,  but  measuring  probes 
with  two  sensing  elements  perpendicular  to  one  another  are  available 
for  measurement  of  the  turbulence  intensity  in  three  dimensions. 

This  measurement  is  essential  for  the  correct  computation  of  the  flow 
parameters.  For  an  approximate  analysis  of  the  flow  the  turbulence 
term  is  generally  assumed  to  be  negligible  but  for  cases  such  as  flow 
over  'rough'  surfaces  the  turbulence  term  becomes  increasingly  signif¬ 
icant  in  any  subsequent  analysis. 


2.4.4  Optical  methods 


Included  in  this  section  is  the  measurement  of  fluid  velocities 
using  the  Doppler  shift  of  a  sea*  tered  laser  beam,  first  used  and 


suggested  by  Yah  and  Cummins  { 196/) .  Light  scattered  by  a  moving 
particle  causes  a  Doppler  shift  in  frequency,  giving  the  following 
approximate  relation,  when  the  velocity  is  very  small  compared  with  the 
velocity  o£  light 

f  =  (n/A0)  U.  (ej)  -  (eg) . 

n  i  sc  in 


where  n  is  the  refraction  index  '“f  the  fluid,  IK  is  the  fluid  velocity 
vector,  Xg  the  vacuum  wavelength  and  (sj)^  and  (ejK^  are  the  unit 
vectors  in  the  direction  of  the  scattered  beam  and  of  the  incident 
beam  respective lv . 


The  experimental  set  up  is  shown  in  Figure  2-3.  The  light  from 
the  incident  beam  scattered  by  the  fluid  particles,  in  the  direction  of 
the  photomultiplier,  combines  with  tne  non- scattered  reference  beam. 

The  photomultiplier  acts  as  ar.  amplifier  and  as  a  heterodyne  receiver. 
Such  a  receiver  mixes  twe  signals  of  different  frequencies  an**  produces 
an  output  signal  with  a  frequency  equal  to  the  difference  in  the 
frequencies  of  the  input  signals.  Hence  the  receiver  generates  a 
current  whose  AC  component  is  equal  to  the  difference,  Doppler,  fre¬ 
quency.  When  U  is  the  velocity  of  the  fluid  at  the  point  of  inter¬ 
section  of  the  two  beams,  the  difference  frequency  is 


fD  =  •  (2n/>.0)  U  Sin  (6^2) 


I  ;1 


where  angle  0^  is  measured  in  the  fluid. 


The  advantage  that  optical  methods  have  over  all  the  other  tech¬ 
niques  available  is  that  they  do  not  physically  interfere  with  the 
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tlow  and  also  they  are  absolute.  The  instrument  is  ideally  suited  to 
measuring  high-frequency  velocity  fluctuations.  Information  about  the 
turbulence  can  be  obtained  by  utilising  a  spectrum  analyser  whose  output 
is  scanned  and  integrated. 

2.4.5  The  measurement  of  Cf  in  this  investigation 

It  can  be  seen  from  the  preceding  notes  that  tv a  two  most  versatile 
and  appropriate  techniques  for  the  measurement  of  the  coefficient  of 
frictional  resistance,  or  for  that  matter  any  other  flow  parameter,  are 
the  Hot-v/ire/ film  anemometer  and  the  Laser  Doppler-Shift  instrument. 

The  Laser  Doppler-Shift  method  is  more  attractive  than  the  Hot-wire/film 
anemometer  in  that  it  requires  no  calibration  and  dees  not  interfere 
with  the  flow.  In  fact,  the  Laser  method  would  be  ideally  suited  to  the 
task  of  calibrating  the  anemometer  sensing  elements  in  the  working 
fluid.  However,  such  an  instrument  v;as  not  available  for  the  present 
investigation.  Hence  the  hot-film  anemometer  'with  a  single  cylindrical 
senscr  is  used  to  record  the  velocity  distributions  across  the  boundary 
layer  adjacent  to  the  v;ork  surface.  The  anemometer  was  calibrated,  in 
'fully*  developed  smooth  turbulent  flow,  on  the  axis  of  a  horizontal 
glass  flow  tube  connected  to  a  constant-head  water  tank  (Chapter  3} . 


The  resulting  boundary  layer  velocities  are  analysed  using  the 
streamwise  momentum  integral,  given  by 
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IS 


Thus,  using  the  computed  values  of  0  and  H  and  graphically 

determined  values  of  dU  /dx  and  d0/d*  over  the  work  surface,  a 

e 

series  of  values  of  are  determined  from  the  two-dimensional 

momentum  integral  equation. 


Corresponding  values  of  Cf(2)  are  computed  (see  Appendix  4) 
assuming  the  validity  of  a  logarithmic  region  in  the  boundary  layer 
velocity  distribution.  Values  of  Cf(3)  are  also  computed  by  ficting 
the  boundary  layer  velocity  distribution  to  a  logarithmic  plus  a  wake 
relationship  given  by 


=  -  log 
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FIG.  2-2  THE  FLOATING  ELEMEN 
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FIG.  2-3  THE  LASER-DQPPLER  ANEMOMETER 


laser 

source 


lens 


'V 

splitter 

plate 


mirror 

<x 

npiftrar 


neutral 
density  filter 


■  tuw 

regime 


rranspare 
sides  of  f  iov 
apparatus 


\  yC  scattered 
N>Cr—  beam 

reference 
i  .  rxbearn 


Trrrmrrr? 

mirror 


W\ 

X\ 


II  Dhotomui 


3 . 1  Introduction 

Preliminary  measurements  of  the  flow,  where  no  great  accuracy  in 
the  velocity  readings  near  the  solid  boundaries  is  required,  were 
carried  out  using  a  pitot-static  tube  and  a  propeller  type  current  meter. 
Flow  velocities  in  the  case  of  the  pitot-static  tube  were  calculated 
using  readings  of  stagnation  and  static  pressures,  resulting  in  an 
absolute  measure  of  velocity. 

The  pitot-static  tube  is  connected  to  a  water  manometer  open  to  the 
atmosphere,  the  difference  in  the  two  pressures  is  measured  in  terms  of 
the  height  difference  in  the  manometer  readings.  The  set  up  was  such 
that  the  manometer  readings  for  th'  stagnation  and  static  pressures 
could  be  measured  separately  using  a  depth  gauge  from  the  open  top  of 
the  manometer  tube.  Arrangements  were  made  to  purge  the  pitot-static 
tube  and  manometer  of  air  locks  and  obstructions  using  an  air  suction- 
line  and  a  pressurised  water  line  (Figure  3-1).  The  pitot-static  tube 
is  held  in  position  by  a  carriage  arrangement  which  allows  it  movement 
in  the  x,  y  and  z  directions.  A  serious  disadvantage  of  using  the 
pitot-static  tube  in  a  liquid  environment  is  the  length  of  time  that  must 
be  allowed  for  the  manometer  reading  to  settle  down  to  an  equilibrium 
position  before  a  reading  can  be  taken. 

An  alternative  to  the  pitot-static  tube  is  a  water  current  meter. 

The  current  meter  used  consisted  of  a  propeller  unit,  comprising  the 
propeller,  the  shaft,  the  rotating  sleeves  and  ball-bearings  and  a 


meter  body,  containing  the  electrical  impulse  generating  unit  plus 


an  electrical  connection  to  an  electrical  counter  which  registers  each 
revolution  of  the  propeller.  A  stop-watch  is  used  to  register  the 
number  of  revolutions  per  second,  n,  and  the  mean  flow  velocity  computed 
from  the  calibration 

U  =  K.n  +  Aj 

The  constants  K  and  Aj  are  predetermined  for  the  particular  propeller 
used.  The  resulting  velocity  profiles  were  used  to  determine  the 
approximate  extent  of  the  side  wall  shear  layers. 

The  flow  analysis  used  necessitated  measuring  the  boundary  layer 
velocity  distribution  very  close  to  the  work  surface.  Hence  the  meas¬ 
uring  probe  used  must  be  as  small  as  possible.  This  is  especially  so 
in  the  region  near  the  surface  where  the  gap  between  the  probe  and 
surface  could  act  to  accelerate  the  flow  through  the  constriction  so 
formed  and  hence  produce  a  change  in  the  flow  pattern  in  the  vicinity 
of  the  measuring  probe.  Both  the  current  meter  and  pitot-static  tube 
cannot  be  relied  upon  tc  give  reliable  velocity  measurements  close  to 
a  solid  boundary  due  to  their  relative  size.  These  instruments  also 
have  the  disadvantage  that  they  do  not  have  an  electrical  output  which 
can  be  electronically  manipulated  to  provide  the  required  fluid  para¬ 
meters.  The  hot-film  anemometry  has  the  advantage  of  an  electrical 
signal  and  the  measuring  probe  is  relatively  small. 

3.2  The  constant  temperature  anemometer 


The  use  of  hot  sensors  for  the  measurement  of  particle  velocity 
in  fluid  flow  relies  on  a  relation  similar  to  the  formula  derived  bv 
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King  (1914)  governing  convective  heat  transfer  from  a  cylinder  in  an 
incompressible  potential  flow.  According  to  King,  the  rate  of  heat  loss 
Q  and  flow  velocity  U  are  related  through 

Q  =  (Aj  +  Bj  /u)  (T  -  T0) 

where  A^  and  are  empirical  constants,  U  is  the  flow  velocity,  T  is 
the  wire  temperature  and  Tq  is  the  fluid  temperature  for  thermal  equili¬ 
brium  conditions.  The  rate  of  heat  loss  from  the  sensor  must  equal  the 
heating  power  produced  in  the  wire  by  the  electric  current  I. 

Incorporating  the  temperature  differences  between  the  wire  and 
fluid  in  the  empirical  constants  A  and  B,  the  heat  loss  formula  takes 
the  form 

IJR/R  -  RQ  =  A  +  B  (U)n 

where  R  is  the  operating  resistance  of  the  hot-film,  RQ  is  the  resist¬ 
ance  of  the  hot-film  at  fluid  temperature  and  values  of  A,  B  and  n  are 
determined  to  give  the  best  fit  to  the  measured  data  within  a  selected 
velocity  interval. 

The  anemometer  output  is  a  bridge  voltage  V,  and  the  squared 
voltages  V2  and  VQ2 ,  where  VQ  is  the  bridge  voltage  at  zero  velocity, 
are  linearly  related  tc  the  heat  loss  of  the  sensor  at  the  velocity  in 
question  and  at  zero  velocity  respectively. 

An  analog  computer  is  used  to  linearise  the  anemometer  so  that  the 
meaurements  of  fluctuating  velocities  and  velocity  profiles  are  not 


affected  by  amplitude  distortion.  The  analog  computer  used  has  the 
following  transfer  function 


V  =  K(V. 
out  m 


V.  ) 
ino 


m 


Thus,  the  non  linear  relationship, 


VB2  =  A  +  BUn 

between  the  output  voltage  V  from  the  hot  wire  or  hot-film  anemometer 
ana  the  measured  flow  velocity  (U)  can  be  linearised  to  read 


with  good  approximation  down  to  velocities  of  20  cm/sec. 

3.2.1  Theory  of  operation  of  the  CTA 

Figure  3-2  shows  a  circuit  diagram  of  the  Disa  CTA  unit  used. 

The  sensing  device  used  consists  of  a  thin  electrically  heated  metal 
film  fused  to  a  quartz  fibre  for  support.  The  fibre  plus  ^ilm  has  a 
thin  insulating  coat  of  quartz  and  is  suspended  between  the  two  probe 
supports . 

The  hot  film  probe  forms  one  arm  of  a  Wheatstone  bridge  circuit, 
the  other  arm  consisting  of  a  'non'  inductive  variable  resistance  box. 
The  bridge  is  powered  by  an  amplifier  whose  output  voltage  is  con¬ 
trolled  by  the  bridge  imbalance.  The  system  operates  as  a  so-called 
null-seeking  device.  That  is,  if  fcr  instance,  the  probe  resistance 
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is  less  than  the  resistance  of  rhe  adjustable  bridge  arm,  then  an  error 
voltage  will  be  produced  across  the  amplifier  inputs.  This  voltage  is 
applied  across  a  very  high  gain  amplifier  and  back  to  the  bridge  top  so 
that  the  current  through  the  bridge  increases,  causing  the  probe  to  heat 
up  until  its  resistance  is  equal  to  the  value  of  the  variable  resistance 
box,  whereupon  the  bridge  is  again  in  balance. 


Any  flow  over  the  probe  will  try  to  cool  it,  resulting  in  a  small 
change  in  the  probe  resistance  which  in  turn  produces  an  error  voltage 
across  the  bridge.  This  voltage  is  amplified  in  the  servo  amplifier 
and  fed  back  to  the  bridge  top,  causing  the  bridge  current  to  increase 
and  the  probe  temperature  to  return  to  its  original  value.  The  voltage 
which  is  fed  to  the  bridge  top  to  maintain  the  probe  temperature  can  be 
related  to  fluid  velocity  by  calibration. 


3.2.2  The  probe  operating  resistance 


The  probe  resistance  is  determined  using  the  following  procedure: 


1.  The  probe  cable  and  support  are  terminated  using  a  shortening 
probe  and  the  bridge  is  adjusted  for  balance  with  the  feedback  loop 
open  by  adjusting  a  potentiometer  which  is  connected  in  series  with  the 
decade  resistance  box. 


2.  The  shortening  probe  is  then  replaced  by  the  actual  probe  to 
be  used  and  again  the  bridge  is  adjusted  for  balance,  this  time  by 
adjusting  the  resistance  box.  The  resistance  reading  on  the  decade 
box  corresponds  to  the  cold  resistance  of  the  probe  Rp  (assuming  a 
bridge  ratio  of  1:1)  whex. 


29 


Rd  =  R  +  R 


R  =  Sensor  Resistance  (at  ambient  temperature) 

O 

R  =  Resistance  of  probe  leads 

L 


3.  The  required  operating  resistance  (R^)  for  the  probe  can  now  be 
set  on  the  decade  resistance  using  the  formula 


Rt  =  R  +  R  (1  +  a  (T  -  T0)) 

O  J_i 


where  is  the  sensor  resistance  at  temperature  T0  (ambient  temp¬ 
erature),  a  =  temperature  coefficient  of  resistance  of  the  sensor 
material  and  T  =  required  operating  temperature  of  sensor  (<150°C). 
The  bridge  is  now  balanced. 


4.  Turning  the  anemometer  control  switch  to  operate,  which  closes 
the  feedback  loop,  causes  the  probe  to  heat  up  until  its  resistance 
reaches  the  required  value  and  the  anemometer  is  ready  for  use. 


3.2.3  Frequency  response 


The  anemometer  is  checked  to  make  sure  the  instrument  is  not 
oscillating  due  to  imbalance  of  the  reactance  on  the  two  sides  of  the 
bridge.  Moreover  since  the  anemometer  is  designed  to  measure  rapidly 
fluctuating  velocities,  the  frequency  response  of  the  system  is  also 
checked. 


The  ideal  way  to  check  system  stability  and  frequency  response 
would  be  to  place  the  probe  in  a  flow  and  to  subject  it  to  a  step 
change  in  velocity.  By  observing  the  anemometer  output  on  an  oscillo- 
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scope  the  impulse  response  of  the  system  could  be  checked.  In  practice 
we  simulate  this  by  feeding  an  electrical  square  sine  wave  signal  to  the 
servo  amplifier  output.  This  causes  a  sudden  increase  in  probe  current 
and  resulting  bridge  imbalance.  The  imbalance  is  amplified  and  fed  back 
to  the  bridge  top  to  restore  balance.  The  anemometer  output  voltage 
meanwhile  is  monitored  on  an  oscilloscope.  Any  instability  in  the 
system  will  result  in  oscillations  on  the  anemometer  response.  Using 
the  square  wave  generator,  the  frequency  response  of  the  CTA  bridge  is 
optimised  with  the  bridge  balanced. 

3.2.4  Probe  contamination 

Before  taking  any  measurements  with  the  hot-film  probe,  it  is 
lightly  brushed  with  a  dilutic  solution  of  acetic  acid  to  remove  any 
impurities  from  its  surface  which  may  alter  the  heat  transfer  charac¬ 
teristics  of  the  sensing  element. 

A  filtration  unit  is  incorporated  in  the  flow  system  to  remove 
any  impurities  which  may  enter  the  system  via  the  channel's  free 
surface  or  from  corrosion  of  the  connecting  pipes.  At  low  velocities, 
the  probe  is  prone  to  air  bubbles  collecting  on  its  surface  which  will 
lead  to  'hot  spots'  on  the  sensing  element  which  say  cause  it  to  burn 
out.  In  the  case  of  air  bubbles  collecting  on  the  probe  surface,  the 
probe  current  is  switched  off,  the  probe  raised  above  the  water,  re¬ 
turned  to  its  original  position,  and  the  measurements  resumed. 

3.2.5  Probe  traversing  technique 

The  hot-film  is  held  in  position  using  a  carriage  similar  to  that 
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used  for  the  pitot-static  tube  and  the  current  meter.  A  scale  alone  the 
length  of  the  channel  plus  scales  on  the  vertical  and  cross-wise  axis 
of  the  probe  make  it  possible  to  relocate  the  sensing  element  to  any 
point  in  the  flov?.  The  probe  guide  tube  is  earthed  so  as  to  prevent 
the  build  up  of  charge  on  the  film  sensor. 


An  off-set  probe  is  used  to  take  measurements  very  close  to  the 
channel  floor;  the  orientation  of  the  probe  is  such  that  the  probe 
measures  the  velocity  fluctuations  in  the  x  and  z  flov/  directions.  The 
hot-film  element  is  primarily  sensitive  to  the  component  of  flov;  velo¬ 
city  normal  to  the  axis  of  the  sensing  element. 


If  the  probe  is  yawed  about  an  axis  normal  to  the  sensing  element 
a  progressively  shorter  length  of  the  element  will  be  exposed  to  the 
flow  cooling  effect  and  the  sensitivity  will  fall.  At  a  yaw  angle  of 
90°,  that  is  when  the  flow  is  parallel  to  the  sensing  element,  there 
will  be  a  minimal  cooling  effect  and  therefore  flov;  components  along 
the  sensor  can  be  neglected.  The  sensitivity  cf  a  hot-film  probe  to 
variations  in  yaw  angle  is  given  by 


Ueff  =  U  COS  al 


where  U  =  actual  flov;  velocity,  =  effective  cooling  velocity  and 

=  yaw  angle. 


Considering  the  probe  in  a  typical  flov;  field,  there  will  be  a 
mean  velocity  U  along  the  iixis  of  the  probe  support  and  fluctuating 
velocity  components  u,  v,  v;  in  the  x,  s  and  y  directions.  As  has 
already  been  discussed,  the  effect  of  flow  along  the  sensing  element 
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in  diameter  and  4.6  in  long  (“  200  diameters)  connected  via  a  globe 
valve  to  a  constant  head  reservoir  which  is  supplied  with  water  pumped 
from  the  channel  system-  Hence,  the  hot-film  probe  is  calibrated  in 
a  situation  where  the  intensity  of  turbulence  is  at  a  minimum  i.e.  at 
the  centre  line  into  smooth  turbulent  pipe  flow. 


Hie  calibration  tube  is  slung  between  the  sides  of  the  flow 
channel  above  the  level  of  water,  so  that  the  probe  can  be  raised 
above  the  channel  water  level  and  inserted  inside  the  outlet  end  of  the 
calibration  tube  and  a  calibration  performed  without  roving  the  probe 
from  the  channel  system. 


Various  flow  rates  can  be  achieved  by  adjustment  of  the  globe 
valve  on  the  outlet  cf  the  constant-head  tank.  The  flow  rate  in  the 
calibration  tube,  and  hence  the  mean  velocity,  are  determined  from 
the  discharge  rate  measured  over  a  minimum  period  of  30  seconds,  the 
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water  collected  in 


a  polythene  bag  and  the  volume  measured.  From  a 
knowledge  of  the  tube  flow  Reynolds  number  assuming  the  probe  to  be 
positioned  on  the  calibration  tube  centre  line,  then  the  relation 
between  the  mean  velocity  and  the  centre  line  velocity  {Figures  3-4  end 
3-5)  can  be  used  to  calculate  the  centre  line  'velocity  (NASA  1963) . 

The  centre  line  velocity  is  then  related  to  the  CTA  bridge  voltage  and 
the  probe  calibrated.  Typical  calibration  curves  are  shown  in  Figures 
3-6  and  3-7.  Calibrations  are  taken  for  each  of  the  Linearieer  unit 
exponent  values  and  best  fit  lines  by  least  squares  fitted  through  each 
set  of  data.  Thus  the  most  appropriate  exponent  value  for  linear  cali¬ 
bration  can  be  determined.  This  value  is  then  used  for  all  the  probe 
measurements. 

3.2.7  Probe  operation 

After  calibration  the  probe  is  immersed  i.i  the  channel  flow,  and 
the  operating  resistance  calculated  and  set  for  the  existing  water 
temperature.  The  probe  is  then  ready  for  the  measurement  of  water  flow- 
velocities  at  the  same  temperature .  If  the  water  temperature  changes 
because  of  the  heating  effect  of  the  pump  or  changes  in  the  atmospheric 
temperature  then  the  operating  resistance  must  be  re-calculated  and 
reset  on  the  CTA  unit,  thus  providing  a  constant  overheat  ratio. 

The  probe  is  manually  positioned  for  the  first  reading,  the  CTA 
unit  and  then  the  data  logger  switched  on,  and  the  linearised  bridge 
voltage  sampled  at  0.33  second  intervals  over  a  period  of  10  seconds. 
The  sampled  voltage  is  recorded  or.  punched  paper  tape  for  a  subsequent 
computer  analysis.  After  the  sampling  period  the  data  logger  is 
switched  off,  the  probe  moved  to  the  next  measurina  ooint  and  the 
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apart)  in  the  floor  shear  layer  region  and  as,  points  further 
(1  cn)  in  the  free-stream  reaion.  This  is  to  enable  an  aeeui 


measurement  of  the  comparatively 


K'Sioc!  tv  ors 


the  shear  laver- 


For  cross-channel  constant  depth  flow  measurements  the  probe  is 
moved  in  equal  comparative iy  larger  (1  inch)  steps  across  the  channel 
width  to  gain  a  general  picture  of  the  side  wall  shear  layer  development 
and  the  central  core  flow.  The  resulting  values  of  mean  voltage  for 
each  of  the  measuring  points  at  a  station,  along  with  a  corresponding 
value  of  free  stream  velocity,  are  converted  into  mean  velocities  using 
the  appropriate  hot-film  calibration. 


The  free  stream  velocity  is  measured  using  the  Ott  current  mater 
which  gives  a  value  of  velocity  integrated  over  the  area  generated 
from  the  propeller  radius,  which  is  then  assigned  to  the  mean  value  of 
measured  anemometer  voltage  in  the  same  region  of  flow. 


3.3  The  Dali  flow  tube 


A  Dali  tube  measuring  device  is  incorporated  in  the  pipeline 
between  the  pumping  station  and  the  channel  inlet  to  give  a  measure  of 
the  water  flow  rate  through  the  channel  at  any  particular  set  of  flow 
conditions. 


The  Dali  tube  is  a  modified  Venturi  tube.  The  tube  inlet  incorp¬ 
orates  a  dam  followed  by  the  inlet  converging  cone  and  then  a  narrow 


cylindrical  section  on  either  side  of  the  throat  slot.  The  flow  then 
encounters  the  recovery  diverging  cone  with  an  included  angle  of  about 
15  degrees  after  which  the  flow  undergoes  sudden  enlargement  to  the 
pipe  diameter.  The  cones  are  steep  in  addition  to  being  truncated,  and 
the  whole  device  is  only  about  2  diameters  long.  There  are  no  gradual 
curves  as  in  Venturi  tubes,  and  other  devices  -which  are  intended  to 
accelerate  or  decelerate  a  fluid  with  minimum  loss. 


The  differential  pressure  produced  by  a  Dali  flow  tube  is  much 
higher  than  that  of  a  corresponding  nozzle  or  Venturi  tube.  The  flow 
rate,  Q,  is  determined  from  the  pressure  differences  created  by  the 
constriction  in  the  meter,  where  Q  is  given  by 


Q  -  0^2 


and  which  is  valid  for  one  dimensional  flow  as  given  by  Bernoulli. 
Vacations  from  the  one-dimensional  theory  are  all  included  in  the 
particular  value  used  for  the  coefficient  C^.  A  typical  calibration 
curve  of  values  is  shown  in  Figure  3-8  and  values  of  for  diff¬ 
erent  Dali  flow  tube,  throat  to  pipe  diameter  ratios,  ,  are  plotted 
in  Figure  3-9 


Pipe  Reynolds  numbers  of  the  order  3.5  x  ICr  and  above  need  to  be 

reached  before  the  value  of  C.  can  be  assumed  to  be  constant.  In  the 

d 

calibration  rig  available  the  maximum  obtainable  Reynolds  number  is 

c 

approximately  3.7  x  10  and  therefore  it  was  not  possible  to  calibrate 
the  Dali  flow  tube  in  its  linear  working  region.  The  value  of  for 
the  Dali  flow  tube  was  therefore  determined  from  Figure  3-9  for  a 
pipe  ratio,  Bi  ,  of  0.643  and  for  Reynolds  numbers  greater  than  3.5  x  10^. 
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4.  CONSIDERATIONS  OP  THE  FLOW  ANALYSIS 


The  flow  analysis  uses  the  momentum  equation  applied  to  a  control 
volume  bounded  by:  (1)  the  channel  floor;  (2)  two  vertical  parallel 
planes,  separated  by  a  distance  Zsy,  which  coincides  with  the  extremities 
of  the  side  wall  shear  layers;  (3)  a  plane  parallel  to  the  water  surface 
but  at  a  distance  below  the  surface  such  that  the  water  surface  tension 
and  interfacial  shear  effects  can  be  neglected;  (4)  two  vertical  para¬ 
llel  planes  across  the  flow  channel  separated  by  a  distance  d* 

Figure  4-1  gives  a  schematic  diagram  of  the  control  volume. 

The  structure  and  the  level  of  turbulence  are  assumed  such  that 
the  contribution  the  resultant  of  the  normal  turbulent  stress  terms 
make  to  the  momentum  equation  is  deemed  to  be  small  by  comparison  with 
the  other  terms.  For  constant  depth  flow  the  resultant  of  the  forces 
Fj  and  F2»  on  the  ends  of  the  control  volume,  which  are  equal  and 
opposite,  is  zero,  and  the  component  of  the  normal  forces  in  the  direc¬ 
tion  of  fluid  motion  is  likewise  negligibly  small. 

The  weight  W  of  the  liquid  in  the  control  volume  has  a  component 

W  sin  a2  in  the  direction  of  motion,  where  (*2  is  the  angle  of  slope 

of  the  channel  floor.  The  maximum  slope  encountered  in  the  present 

-5 

investigation  gives  the  maximum  value  of  sin  ct2  as  1  x  10  ,  hence 

the  maximum  value  of  the  component  of  weight  in  the  direction  of  motion 
is  1  x  10-5W. 

Consider  that  part  of  the  control  volume  bounded  by  the  outer  edge 
of  the  channel  floor  shear  layer  (not  a  stream  surface) ,  and  apply  the 
continuity  equation  and  momentum  theorem  to  a  stream-wise  slice,  of 
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unit  traverse  width,  of  the  volume .  The  mass  flow  through  the  left  hand 
end  of  the  control  volume  is 


and  through  the  right  hand  boundary  surface  it  is 


(4-1) 


ra  +  dm  dx  =  i 
ax 


pudz  +  d_  I  Ipudz  dx 
J  dx  ^  -» 

0  0 


(4-2) 


The  difference  between  the  two. 


dm  dx  =  d_  f  f 
d*  dx  _  j1 


pudz  dx 


(4-3) 


must,  by  the  continuity  equation,  be  equal  to  the  flow  through  the  top 
boundary  surface.  The  x  component  of  momentum  flux  through  the  left 
hand  side  is 


H*  ‘  f 


pu1  dz 


and  through  the  right  hand  side  it  is 


I  +  dKx  dx  =  jo u*  dz  +  d  jfpus  dz] 

x  aT"  J  U  J 


the  difference  being 

6 

dM  ,  ,  [7  s-1  - 

_ x  dx  =  d_  I jpu  uzj  ax 

dx  ux  h)  J 

0 

The  x  comoonent  over  the  tor  surface  is 


(4-4) 


(4-5) 


(4-6) 


M  )  =  U  d  ff  ,1 

x  top  e  —  Ipuoz 

Q  U  J 


(4-7) 


The  only  forces  acting  on  the-  control  volume  are  assumed  to  be 
shear  forces  at  the  p] ate  (-tq  dx  )  and  the  component  of  the  v;eight 
(V'/As)  sin  G.  The  moin.entu.vi  law  for  steady  flow  states  that  the  sun  cf 
all  the  forces  acting  on  a  volume  in  a  given  direction  is  equal  to  the 
next  efflux  rate  of  momentum  from  this  volume  in  the  same  direction. 
Applying  this  law  to  the  control  volume,  then  in  the  x  direction 


*  [W] 


dx  -  U 


e  fix 
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=  -To  d*  +  (VJ/uz)  Sin  G 


On  substituting  typical  experimental  values  into  the  momentum 
equation,  the  effect  of  gravity  on  the  control  volume  is  found  to  be 
several  orders  of  magnitude  less  than  the  orher  force  components  and  is 
therefore  neglected  in  the  analysis. 

For  steady  uniform  incompressible  turbulent  flow ,  the  stream-wise 
boundary  layer  equations  can  be  approximated  by 


2H  +  *£  =  o 
d*  a* 


continuity 


(4-9) 


u  d\j  v  du 
d*  dz 


„  dU  ,  d  r  du  — ; — -  \  .  d  r  — — .«>  ,, 

=  U  e  +  —  v  - - u’v'  i  -f  —  I  v'“  -  u,2j  (4-i0) 

e  - —  dz  ^  dz  '  ax  1  J 
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d  (  v':  -  u 1 2 ) 
x  -  momentum 

The  turbulent  normal  stresses,  given  bv  d  t  —  3 

dx^  7  "  u  >• 

although  commonly  neglected  can  contribute  as  much  as  10  per  cent  to 
the  magnitude  of  the  momentum  equation  at  flow  separation.  Substituting 
T  «  n(dvi/dz)  -  ou'v'  into  the  x  momentum  equation  and  neglect¬ 

ing  the  turbulent  normal  stress  yields 


-a 


—  du  ,  —  du 

u  dx  v  dz 
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e  dx 


1  dx 
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The  continuity  is  used  to  eliminate  v(x,z)  in  favour  of 

u(x,z)  and  then  the  momentum  equation  integrated  with  resc-uct  to 
across  the  entire  boundary  layer  to  give  the  ordinary  differential 
eauation 


33  .  G  dU 

_  +  (2  T  h)  —  ^ 
e  d* 


since  bv  definition  C  =  — pr 

f  *2plJ‘ 


PV  "  2- 


:  o’  (1  -  r  dz 

0  e  e 


u  u 

where  ©  =  momentum  thickness  -  (1  -  — )  dz 

e  e 

=  momentum  shape  factor  =  5i/0 
=  displacement  thickness  =  (1  -  — —  )  dz 

nJ  K 


H  =  momentum  shape  factor  =  &\/0 


(4-12) 


(4-13) 


(4-14) 


(4-15) 


(4-16) 


Thus,  for  the  case  of  parallel  side  wall  shear  layers,  the  velocity 
distribution  is  measured  perpendicular  to  the  channel  floor  and  at 
various  stations  along  the  channel  centre-line  in  the  control  volume 
ABCD  (Figure  4-1).  From  the  velocity  distributions  the  values  of  0, 

6 l  ,  Ue  and  H  are  evaluated  and  these  substituted  into  equation  4-4 
to  give  the  local  values  of  the  surface  friction  coefficient  at  each  of 
the  measuring  stations  in  the  control  volume. 


The  analysis  is  modified  to  deal  with  the  case  of  developing  side 
wall  shear  layers  by  the  addition  of  a  convergence  term,  ©/  (xc  -  >:)* 
to  the  right  hand  side  of  equation  4.  *c  is  the  location  of  the  point, 
as  measured  from  the  channel  leading  edge,  where  the  side  wall  shear 
layers  merge.  With  the  convergence  of  the  control  volume  there  will  be 
a  resultant  streamwise  hydrostatic  pressure  force?  this  is  assumed  to 
be  neglifible. 


*  ?roc.  Computation  of  Turb.  Bound  Lays* 
AF03R  -  IFP  -  Stanford  Cont.  Vcl  1. 
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The  results  from  the  momentum  integral  analyses  are  compared  with 

corresponding  fluid  parameters  computed  from  the  existence  of  a  loga- 

♦ 

rithinic  relationship  for  boundary  layer  flow.  The  logarithmic  region 
is  formed  by  the  overlap  of  the  inner  and  outer  laws  for  two  dimensional 
turbulent  boundary  layer  flow  over  smooth  surfaces.  The  overlap  between 
the  two  laws  is  given  by 


u  1  ,  zu 

—  =  —  log  T  +  B 

u  k  e  - 

t  v 


(4-17) 


The  validity  of  the  inner  law,  including  the  logarithmic  overlap 
can  be  established  by  plotting  velocity  profile  data  in  terms  of  the 


inner  variables  u/u  and  zu  /v  .A  single  olot  of  several  sets 

T  T 

of  data  produced  by  workers  using  widely  different  flow  regimes  (Figure 
4-2)  clearly  shows  the  data  with  the  exception  of  separating  flow 
collapsing  onto  a  single  logarithmic  relation  in  the  overlap  region 

ZU 

35  <  _ r  <  350  ,  corresponding  roughly  to  0.02  <  z/5  <  0.2, 

v 

after  which  the  individual  data  curves  either  veer  off  upwards  in  the 
outer  (wakelike)  layer  or  downwards  in  the  inner  viscous  shear  layer. 
The  inner  raw  is  extended  into  the  outer  region  of  the  boundary  layer 
by  the  addition  of  the  wake  function  (Coles  1956) ,  thus 


u  l  ,  zu  _  .  JI  W  rz\ 

u  k  e  -  <  'V 


(4-18) 


For  flow  over  smooth  surfaces  the  surface  roughness  is  entirely 
submerged  in  the  viscous  sublayer  and  the  flow  adjacent  to  the  surface 
follows  the  contour  of  the  surface  roughness  almost  without  separation. 
With  'increasing  surface  roughness'  the  roughness  elements  will  begin 
to  break  through  the  viscous  sublayer  producing  localised  flew  separa¬ 
tion  and  eddy  shedding.  The  viscous  sublayer  will  become  increasingly 
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less  stable  until  the  roughness  elements  are  such  that  the  sublayer 
will  merge  with  the  outer  layers. 

Thus,  the  effect  of  surface  roughness  on  the  shear  flow  is  to 

modify  the  mean  velocity  distribution.  The  effect  on  the  senilogarith- 

mic  plot  cf  u/u  versus  log„  zu  /v  is  to  disnlace  it  downwards  from 
x  ex 

the  smooth  wall  profile  by  an  amount  Au/u  .  This  downward  shift  of 

~  T 

the  logarithmic  profile  due  to  surface  roughness  had  been  noted  by 
Nikuradse  (1933)  for  the  case  of  pipe  flow. 

For  smooth  surfaces  a  family  of  straight  lines  can  be  produced, 
each  line  corresponding  to  a  given  value  of  Cf -  Hence,  experimental 
points  can  be  plotted  on  such  a  plot  (Figure  4-3)  and  the  line  upon 
which  they  fall  gives  the  appropriate  value  of  C-.  The  resulting  value 
of  Cf  can  be  checked  against  both  the  slope  and  position  of  the  line 
(Clauser  1954) . 

Extending  this  method  to  rough  surfaces  produces  several  diffi¬ 
culties,  these  being  that  the  origin  of  z  is  not  known  and  also,  since 
the  roughness  causes  a  shift  in  the  logarithmic  profile,  the  value  of 
Cf  is  confirmed  only  by  the  slope  of  the  logarithmic  line  and  not  its 
position.  By  defining  a  height  c ,  measured  from  the  crests  of  the 
roughness  elements,  where  the  logarithmic  asymptote  is  located.  Perry 
and  Joubert  (1963)  were  able  to  extend  the  method  to  the  case  of  rough 
surfaces. 

Using  this  method  of  shifting  the  origin  of  the  mean  velocity 
distribution  to  get  the  best  fit  line  the  the  sexilogarithsic  plot, 
the  relationship  is  given  by 
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where  the  roughness  function  g  is 


_  1  ,  ku  ,  _  Au 

g  -  —  log  t  +  B  -  — 

k  e -  u 

v  T 


(4-20) 


for  large  values  of  the  roughness  Reynolds  number  R^  (R^  =  ku  /v) 


In  the  search  for  a  universal  roughness  function,  the  roughness 

Reynolds  number  R,  has  been  expressed  in  terms  of  various  roughness 
K 

height  parameters.  A  characteristic  roughness  height  in  terms  of  the 
root  mean  squared  roughness  is  used  by  Karlsson  (1978),  the  mean 
apparent  amplitude  of  the  roughness  and  centre  line  average  have  been 
used  by  other  workers  (Musker  et  al  1976) . 


The  inability  of  these  parameters  to  provide  a  universal  relation¬ 
ship  between  the  roughness  Reynolds  number  and  the  roughness  function 

Au/u  led  to  the  use  of  various  simple  spatial  parameters  in  the 
x 

characterisation  of  the  surfaces  under  investigation  (Betterman  1966, 
Dvorak  1969,  Simpson  1973).  Several  empirical  relationships  were 
subsequently  produced,  all  applying  to  the  special  case  of  the  fully 
rough  flow  regime. 

In  a  study  of  flow  through  pipes  lined  with  replicated  ship  hull 
surfaces  Musker  1977  set  out  to  determine  a  roughness  function  using 
parameters  from  a  i  atisticai  analysis  of  the  surfaces.  The  roughness 
Reynolds  number  i:.  edified  to 

kux  auT  (4-21) 
— -  -  — ( 1  +  ASp)  ( 1  +  bS.  ku) 

V  V  K 

(Musker  &  Lewkowicz  1978) 
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where  a  is  the  standard  deviation  of  the  surface  heights,  Sp  is  the 
average  slope  based  on  a  sampling  interval  equal  to  the  correlation 
length  and  'a'  and  'b'  are  empirical  constants,  the  parameter  Sp  is 
assumed  to  represent  the  degree  of  bluff- iess  of  the  roughness  elements, 
and  SK  and  are  the  average  skewness  and  kurtosis  of  the  surface 
height  distributions. 

For  the  surface  investigated  an  approximation  to  a  universal  curve 
for  the  roughness  functions  could  be  obtained  if  the  short  wavelength 
cut-off  was  set  close  to  the  extrapolated  value  of  the  average  turbu¬ 
lence  macro-scale  in  the  vicinity  of  the  roughness  elements.  Hence  for 
a  2  mm  cut-off  wavelength  the  values  a  =  0.5  and  b  =  0.2  produce  a 
nearly  universal  curve  (Figure  4-4)  for  the  roughness  function, 
plotter  against  the  modified  roughness  Reynolds  number. 


5.1  Introduction  to  surface  measurement 


There  are  many  instruments  available  for  the  measurement  of  surface 
geometry.  These  are  capable  of  measuring  surface  features  from  meters 
in  length  as  in  the  case  of  a  surveyor's  theodolite  right  down  to  looking 
at  the  molecular  level  as  measured  by  the  Scanning  Electron  Microscope. 

A  state  of  the  art  review  and  bibliography  of  Surface  Topography  in 
Engineering  is  given  by  Thomas  and  King  (1977). 

For  a  successful  analysis  of  surface  roughness  it  is  important 
that  the  measuring  instrument  used  produces  readily  available  quanti¬ 
tative  information  about  the  surface,  it  clearly  being  an  advantage  if 
the  data  is  in  the  form  of  an  electrical  signal  which  lends  itself  to 
electronic  manipulation  and  computer  analysis. 

It  is,  of  course,  important  that  the  measuring  device  is  capable 
of  supplying  information  about  surface  parameters  which  together  are 
capable  of  uniquely  defining  a  particular  surface. 

5.2  Measuring  techniques  for  surface  microgeometry 


The  techniques  available  for  looking  at  surface  micro-geometry , 
for  example  (1)  Optical  Methods?  (2)  Taper-Sectioning?  (3)  Pneumatic 
Gauging  and  (4)  Stylus  Instruments* 
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5.2.1  Optical  methods 


Optical  methods  include: 

(a)  Microscopic  examination  of  the  surface  with  normally  incident 
illumination  which  at  low  magnifications  gives  a  general  picture  of  any 
trends  in  the  surface  features.  At  higher  magnifications  it  becomes 
increasingly  difficult  to  focus  on  the  peak  and  valley  features  simul¬ 
taneously  thus  severely  restricting  the  scope  of  the  microscope  in  the 
examination  of  surface  roughness. 

(b)  The  Light  Section  microscope  partially  overcomes  the  problem 
of  focussing  on  the  peak  and  valley  features  simultaneously  by  throwing 
light  from  a  slit  onto  the  surface  at  an  incident  angle  of  45°.  The 
reflected  image  of  the  slit  is  seen  in  effect  as  a  section  through  the 
surface  'with  the  vertical  scale  magnified  by  a  factor  equal  to  the 
cosecant  of  the  angle  of  incidence. 

(c)  The  Interference  microscope  can  be  used  for  relatively  smooth 
surfaces,  it  showing  the  surface  height  variations  as  light  and  dark 
interference  fringes.  Under  suitable  conditions  it  cap.  give  a  complete 
contour  map  of  the  surface.  None  of  these  optical  methods  are  capable 
of  providing  quantitative  information  without  laborious  processing. 

(d)  The  Transmission  Electron  microscope  has  a  much  greater  res¬ 
olution  and  depth  of  field  than  its  optical  counterparts,  but  it  re¬ 
quires  a  tedious,  delicate  two-stage  replication  technique,  with  a 
consequent  danger  of  poor  reproduction  of  the  initial  surface. 
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(e)  The  Scanning  Electron  microscope  is  a  considerable  improve¬ 
ment  of  the  Transmission  microscope.  It  is  a  reflection  microscope  in 
which  an  electron  beam  is  modulated  to  scan  the  surface  in  a  series  c~ 
parallel  traverses.  The  reflected  beam  is  detected,  amplified  and  d: 
played  on  a  TV  screen.  Specimens  can  be  examined  directly;  but  the 
size  of  specimen  is  limited  by  the  size  of  the  specimen  vacuum  chambe 


(f)  Various  designs  of  glossmetsr  have  been  produced  (Westburg 
1967/8)  to  measure  the  reflectance  of  a  surface  to  give  a  relative 
measure  of  the  surface  average  roughness,  h  development  of  this 
technique  is  the  laser  scanning  analyser  {Clarke  &  Thomas  1979) .  A 
laser  beam  is  reflected  from  a  rotating  polygonal  mirror  in  order  to 
scan  the  workpiece’  surface  at  high  speed,  where  it  is  again  reflected 
into  a  wide  aperture  photodetector  which  measures  the  reflected  inten¬ 
sity.  The  detector  output  is  amplified  and  applied  to  the  vertical 
deflection  coils  of  an  oscilloscope  whose  time  base  is  provided  by  the 
rotation  of  the  mirror.  Localised  variations  in  the  reflectance  of  the 
surface  thus  appear  as  changes  in  signal  strength  whose  position  on 
the  workpiece  can  be  established  from  the  time  elapsed  since  the  start 
of  the  current  scan.  The  laser  beam  spot  diameter  can  be  set  from 
200  pm  upwards  at  623  ram.  scan  width,  the  angular  traverse  is  +30°  and 
the  scanning  speed  is  5  KHz  maximum. 


5.2.2  Taper  sectioning 


In  taper  sectioning  the  specimen  is  sectioned  at  a  shallow  angle 
to  the  surface  and  the  section  polished  and  examined  microscopically. 
The  result  is  a  profile  with  its  vertical  scale  exaggerated  by  a  factor 
of  cotangent  a  ,  thus  emphasising  the  surface  features.  The  method  has 
the  disadvantages  of  being  destructive  and  time  consuming. 


. .  . . 


In  pneumatic  gauging  air  is  allowed  to  escape  through  the  gaps 
between  a  flat  nozzle  and  the  surface  it  is  in  contact  with.  The  back 
pressure  is  a  function  of  the  surface  roughness.  The  gauge  is  in¬ 
sensitive  to  smooth  surfaces  and  does  not  give  quantitative  information. 

5.2.4  Stylus  instruments 


The  basic  form  of  the  stylus  instrument  is  that  a  fine  stylus  is 
drawn  at  a  constant  speed  over  an  irregular  surface  and  the  variation 
of  its  vertical  displacement  with  time  will  represent  the  variation  in 
height  of  a  section  of  the  surface  with  horizontal  displacement,  hs 
the  average  slopes  of  most  surfaces  are  small,  the  vertical  movement  of 
the  stylus  is  generally  electronically  amplified. 

The  one  main  advantage  of  the  stylus  instrument  which  makes  it 
superior  to  the  other  surface  instruments  is  the  availability  of  an 
electrical  signal  which  can  be  electronically  manipulated  to  give  any 
desired  roughness  parameter,  displayed  directly  on  a  chart  recorder, 
or  analysed  by  computer. 

5.3  Measuring  techniques  for  surface  macrogeomatry 

The  previously  mentioned  instruments  all  measure  surface  features 
of  horizontal  dimensions  up  to  a  few  millimetres.  In  the  experimental 
set-up  of  this  investigation  the  work  surface  is  approximately  20  metres 
in  length  and  therefore  to  characterise  the  surface,  surface  features 
of  up  to  20  metres  must  be  measured.  The  techniques  available  to 


64 


measure  the  macrcgeometry  of  the  surfaces  car.  be  divided  into  (1}  Optica 
techniques  and  (2)  other  techniques. 

5.3.1  Optical  techniques 

Optical  instruments  for  the  measurement  of  surfaces  include  the 
auto-collimator  and  micro-alignment  telescope.  Both  of  these  use  a 
moveable  target  which  is  placed  on  the  work  surface  and  whose  displace¬ 
ment  relative  to  the  sighting  instrument,  with  distance  between  the  two 
gives  an  indication  of  the  surface  geometry.  The  instrument’s  sampling 
interval  and  sample  length  are  determined  by  the  separation  of  the 
target's  feet,  the  limit  of  angular  deflection  of  the  light  source  det¬ 
ectable  by  the  sighting  device,  and  the  optical  efficiency  of  the 
instrument. 

5*3.2  Other  techniques 

The  other  methods  include  using  instruments  such  as  dial  gauges 
which  can  be  used  to  measure  spot  heights  along  the  surface  relative  to 
a  datum,  such  as  an  engineer's  straight  edge.  An  electronic  level 
measures  the  inclination  of  the  surface  in  steps  equal  to  the  separation 
of  its  feet,  its  range  of  surface  feature  measured  being  frera  the  mean 
distance  apart  of  its  feet  to  as  far  as  it  is  physically  possible  to 
move  the  instrument  along  the  surface.  The  values  of  surface  inclin¬ 
ations  can  then  be  easily  converted  into  height  readings. 

Another  instrument  w-d  for  measuring  surface  'macrcgeometry'  is  a 
coarser  version  of  the  stylus  measuring  device  used  for  measuring  the 
microgeometry  of  a  surface.  3SRA  have,  for  many  years,  used  such  an 


instrument  (Canham  1956)  to  take  in-situ  ship  hull  roughness  records 
when  ships  undergoing  trials  are  in  dry  dock.  The  instrument  in  question 
is  called  the  BSRA  Wall  Roughness  Gauge,  a  diagram  of  which  is  shown  in 
Figure  5-1.  The  gauge  consists  of  a  measuring  probe  mechanically  linked 
to  a  sapphire  stylus  which  scratches  the  measured  profile  onto  a  graph¬ 
itic  coated  glass  slide.  The  measuring  probe  assembly  is  held  in  a 
carriage  mounted  cn  wheels  and  traversed  across  the  hull  surface  on  a 
track  which  is  held  in  position  on  the  surface  by  means  of  two  magnets, 
one  at  either  end  of  the  track.  The  measuring  probe  which  has  a  i.6  a a 
diameter  ball  point  is  drawn  over  the  surface  in  between  each  side  of 
the  track.  The  probe  is  constrained  to  move  in  a  vertical  direction 
and  hence  as  it  is  drawn  across  the  surface,  its  movement  relative  to 
the  track-way  is  recorded  by  the  sapphire  stylus  on  the  graphite  coated 
glass  slide.  The  vertical  measurements  on  the  recorded  profile  are 

equivalent  to  the  vertical  movement  of  the  measuring  probe,  whilst  the 

* 

traverse  distance  has  been  reduced  by  a  factor  of  approximately  xlO  to 
fit  the  profile  onto  the  glass  slide.  Thus,  the  maximum  length  that 
can  be  traversed  by  the  measuring  probe  is  approximately  63.5  cm  and 
the  corresponding  recorded  profile  is  approximately  6.4  cm  long. 

The  wall  gauge  records  can  be  projected  optically  on  to  a  screen 
and  the  Mean  Apparent  Amplitude  (MAA)  for  a  surface  ’wavelength*  of 
2  inches  (5.1  cm)  manually  determined  and  taken  as  a  typical  measure 
of  the  surface  roughness.  The  reason  for  using  a  2  inch  length  being 
that  a  statistical  analysis  of  ship-model  data  (Chaplin  1967)  seemed 
to  indicate  that  a  relationship  exists  between  frictional  resistance 
parameters  based  on  the  MAA  per  2  inch  length  of  hull  surface  and  the 
corresponding  computed  snip-model  correlation  factor.  The  wall  gauge 
surface  profiles  can  also  be  digitised  and  statistically  analysed  by 
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computer.  Prior  to  the  analysis  the  optically  enlarged  profiles  are 
projected  onto  a  digitising  table  connected  to  a  paper  tape  punch. 

BSRA  have  recently  produced  an  electronic  Hull  Roughness  Analyser 
to  replace  the  roughness  gauge  which  has  been  in  use  for  over  30  years. 
The  new  portable  instrument  provides  a  digital  record  on  a  strip-printer 
of  peak  to  valley  roughness  amplitudes  in  alternate  50  mm  lengths  as  the 
sensor  is  propelled  manually  over  the  work  surface.  The  measuring  probe 
has  a  vertical  measuring  range  of  between  15  pm  and  2  mm. 

5.4  Surface  reolication 


The  use  of  surface  replication  in  the  measurement  of  surface 
topography  has  been  recently  evaluated  by  George  (1979) .  The  investi¬ 
gation  makes  use  of  the  power  spectral  density  function  to  compare 
various  replica  materials  and  tecliniques.  The  requirement  of  the 
replica  is  that  'the  geometric  shape  of  the  contact  face  should  be  as 
similar  as  possible  to  the  parent  surface'.  Some  replica  materials, 
such  as  Acrulite  Microtech  type  A,  are  known  to  distort  as  they  set 
causing  bowing  of  the  replica.  For  many  applications  the  bowing  is 
relatively  unimportant  compared  with  the  roughness  replication.  3owing 
cam  however  cause  considerable  inaccuracy  when  investigating  the 
relatively  long  surface  wavelengths. 

It  was  found  that  for  all  the  replica  methods  investigated,  the 
replication  resolution  dropped  off  significantly  for  surface  wave¬ 
lengths  less  than  30  pm.  Bowing  of  the  replicas  was  considerably 
reduced  by  inserting  a  steel  plate  in  the  replica  as  it  curved. 
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5.5  Surface  characterisation 


Surface  topography  is  generally  divided  arbitarily  into  three 
classifications  according  to  its  scale  (Figure  5-2) .  The  largest 
wavelengths  are  classified  as  errors  of  form,  intermediate  wavelengths 
as  waviness,  and  else  short  wavelengths  form  the  surface  roughness  - 
Clearly  the  longest  wavelength  present  can  be  attributed  to  the  nominal 
size  of  the  work  surface.  The  shortest  wavelength  present  is  sore 
difficult  to  define;  theoretically  it  can  be  the  atomic  dimensions  of 
the  surface-  The  wavelength  cut-offs  of  the  measured  profile  must  be 
taken  into  account  in  the  subsequent  analysis.  The  measured  parameters 
of  a  surface  are  dependent  both  on  the  type  of  surface  under  investi¬ 
gation  and  on  the  wavelength  cut-o^s  of  the  measuring  instrument.  For 
a  stylus  instrument  the  short  wavelength  limit  is  determined  by  the 
finite  dimensions  of  the  stylus,  acting  as  a  low -pass  filter,  and  the 
long  wavelength  limit  by  the  sampling  length. 

The  surface  parameters  derived  from  a  computer  analysis  of  the 
surface  profile  are  not  intrinsic  properties  of  the  profile  since  all 
the  digital  computers  in  effect  sample  the  profile  at  discrete  inter¬ 
vals  to  obtain  a  discrete  set  of  height  readings  defining  the  profile. 
Hence  the  computed  surface  parameters  are  a  function  of  the  finite 
dimensions  of  the  measuring  head,  the  sampling  interval  and  tne  sampling 
length  (Thomas  &  Sayles  1978) . 

The  first  task  in  analysing  a  surface  pr'-'lle  is  tc  define  a  mean 
line  or  datum  from  which  parameters  can  be  measured-  This  will  make 
it  possible  to  define  the  surface  height  variation,  but  a  means  is  also 
required  of  defining  how  the  heights  vary  in  the  plane  of  the  surface. 


Thus*  to  descriije  a  surface  at  least  two  parameters  are  gsneraliv 
required,  one  defining  the  roughness  height  and  the  other,  a  spatial 
parameter,  defining  how  the  roughness  varies  along  the  surface.  By 
assuming  the  surface  profile  to  be  tending  to  a  straight  line  parallel 
to  the  measuring  instrument's  datum  then  a  number  of  elementary  para¬ 
meters  can  be  used  to  describe  the  roughness  amplitude.  For  example, 
the  maximum  peak  to  valley  height  of  the  roughness,  Rt,  is  defined  as 
the  vertical  distance  between  the  highest  peak  and  lowest  valley  on  the 
profile.  It  is  easy  to  define  but  has  the  disadvantage  of  being  sens¬ 
itive  to  a  typical  surface  event  such  as  scratches.  The  highest  peak 
and  lowest  valley  on  a  profile  as  measured  by  a  stylus  instrument  are 
not  necessarily  the  highest  and  lowest  points  on  the  surface.  It  is 
unlikely  that  any  peak  on  a  profile  will  coincide  with  any  summit  of 
a  surface  asperity;  the  measuring  instrument  will  probably  follow  a  line 
that  takes  it  over  the  shoulders  of  surface  asperities.  The  two  rest 
widely  used  mean-line  roughness  parameters  are  the  root  mean  square 
(RHS) ,  Hq,  and  centre  line  average  Ra,  each  of  which  car,  be  used  to 
define  the  mean  line  or  plane  through  a  profile.  The  rms  mean  line  is 
produced  by  minimising  the  sun  of  the  squares  of  the  distances  above 
and  below  the  mean  line;  The  rms  roughness  is  defined  by 


where  z  is  measured  from  the  mean  line  and  L  is  the  profile  length. 


The  centre  line  average,  ?.a,  is  the  most  universally  used  rough¬ 
ness  parameter  because  of  its  ease  of  measurement.  The  mean  line  or 
plane  is  defined  by  the  condition  of  equal  areas  or  volumes  of  surface 


and  void,  anovo  and  below  the  mean 


The  R_  parameter  is  defined  for  a  profile  length  X, 


Ra  ~  [  I  sidy 

y  J 

0 


(5-2) 


Since  the  rms  roughness  is  weighted  by  the  square  of  the  surface 
heights  it  tends  to  be  more  sensitive  than  Ra  to  large  deviations  from 
the  mean  line.  In  some  cases  this  can  represent  an  advantage  in  the 
rms  over  the  Ra  parameter. 


A  roughness  measure  used  to  specify  ship  hull  roughness  is  the 
mean  apparent  amplitude,  MAA  (Chaplin  1967) .  Using  the  BSRA  wall  gauge 
a  sample  length  of  approximately  30  inches  (76  cm)  is  recorded  and 
divided  into  2  inch  (5.1  cm'  lengths.  The  maximum  and  minimum  points 
in  the  middle  13  x  2  inch  lengths  are  then  used  to  plot  a  peak  to 
valley  envelope  over  the  total  sample  length.  The  MAA  is  the  difference 
in  height  between  the  average  of  the  13  highest  points  on  the  envelope 
and  the  13  lowest  points. 


Statistical  techniques  were  used  to  evaluate  such  parameters  as 
asperity  peaks,  curvatures,  surface  heights,  slopes  etc.  and  form 
probability  distributions  representing  the  frequency  of  occurence  of 
each  parameter  (Greenwood  &  Williamson  1966) .  The  parameters  so  ob¬ 
tained  do  not  contain  spatial  information  on  the  variation  of  surface 
properties  about  the  plane  of  the  surface.  The  representation  of 
surface  parameters  spatially  originated  from  random  process  theories 
developed  for  the  analysis  of  random  noise  (Rice  1944/5) .  This  tech¬ 
nique  was  applied  successfully  to  surface  profile  analysis  by  Whitehouse 
and  Archard  (1970)  who  express  the  spatial  properties  of  the  profile  { 
in  terms  of  the  autocorrelation  function. 
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The  use  of  statistical  and  random  process  theories  for  representing 
the  surface  as  opposed  to  the  profile  was  due  originally  to  the  classic 
work  on  sea  waves  by  Longuet  Higgins  (1957,  1962).  This  was  later 
adopted  to  solid  surface  topography  by  Nayak  (1971). 


The  probability  or  distribution  function  p(£)  represents  the 
probability  that  the  random  variable  Z(y)  takes  a  value  less  than  or 
equal  to  £  .  The  probability  function  is  the  integral  of  the  prob¬ 
ability  density  function  p(z), 


p(z)  dz 


where  £  is  the  measured  height  and  Z  is  the  corrected  profile  height. 


The  shape  of  the  probability  density  function  in  terms  of  the 
moments  of  the  function  can  be  used  to  provide  information  on  the  be¬ 
haviour  of  the  process  under  investigatioi  .  The  first  moment  M^,  where 


oo 


p(z)  dz 


i  jSP 

5  I 


is  the  mean  or  DC  level  of  the  profile.  This  is  generally  removed 


1 

J4§i 

m 


prior  to  data  processing  and  is  therefore  zero. 


The  second  moment  M2  is  the  variance  o2  which  is  the  square  of  the 
rras  value  or  standard  deviation.  It  is  a  measure  of  the  spread  of  the 
distribution  and  therefore  of  the  data. 


The  third  moment  M3  is  the  skewness,  which  is  useful  for  defining 
variables  with  an  asymmetric  spread.  It  is  conventional  to  normalise 


i  1 
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Kg  thus 


z3  p(z)  dz 


(5-5) 


where  a  is  the  standard  deviation  of  p(Z).  The  symmetric  Gaussian  dis¬ 
tribution,  which  is  a  good  approximation  for  many  surfaces'  structures 
(Greenwood  a  Williamson  1966)  has  a  skewness  of  zero. 


The  fourth  moment  is  called  the  kurtosis  and  represents  the 
peakedness  of  the  distribution.  Distributions  having  a  kurtosis  greater 
than  k  =  3  are  called  Leptokurtic  and  those  with  a  kurtosis  k  <  3  are 
termed  Platykurtic,  where  k  is  the  normalised  fourth  moment 


k  = 


z4  p(z)  dz 


(5-6) 


The  symmetric  Gaussian  distribution  has  a  kurtosis  of  k  =  3. 


The  use  of  these  higher  moments  to  characterise  a  height  distri¬ 
bution  is  however  limited  because  in  assessing  the  significance  of  their 
numerical  values  it  is  necessary  to  remember  that  they  are  themselves 
distributed  as  the  result  of  random  sampling  (Thomas  &  King  1978) . 


The  statistical  representation  of  the  surface  or  profile  heights 
in  terms  of  the  probability  density  function  p(Z)  is  known  as  the 
height  distribution.  In  practice  many  surfaces  have  symmetrical 
Gaussian  distributions  of  the  form 


p(Z) 


1  -  (z3  /2o3  ) 

- -  ® 

c/2-n 


(5-7) 
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The  distribution  is  constructed  by  plotting  the  number  or  propor¬ 
tion  of  surface  heights  lying  in  a  height  interval  as  a  function  of  the 
height.  A  convenient  method  of  testing  for  'normal'  or  Gaussian  height 
distributions  is  to  plot  the  cumulative  distribution  on  probability 
graph  paper;  this  is  scaled  such  that  a  Gaussian  distribution  will 
produce  a  straight  line.  The  knowledge  that  a  distribution  is  Gaussian 
is  useful  particularly  when  dealing  with  random  process  models.  The 
'asperities  on  asperities'  model  (Archard  1967)  which  accounts  for  both 
height  and  spatial  properties  of  surfaces,  combines  the  basic  distri¬ 
butions  in  a  random  process  theory  approach. 

The  outstanding  problem  in  surface  analysis  is  the  representation 
of  spatial  properties,  i.e.  how  heights,  slopes  etc.  vary  with  distance 
in  the  plane  of  the  surface.  The  spatial  functions,  such  as  the  auto¬ 
covariance  function  or  in  its  normalised  form  the  autocorrelation 
function,  the  structure  function  and  the  power  spectrum,  offer  a  means 
of  representing  the  properties  of  all  the  wavelengths  present. 

The  autocorrelation  function  (ACF)  is  simply  the  autocovariance 
function  (ACVF)  normalised  by  the  profile  sample  variance  a2  (i.e.  the 
square  of  the  rms  roughness) ,  Thus 

R(t)  =  C(t)/o7  (5-8) 

where  T  =  sampling  interval,  where  the  ACVF  is  given  by 

L 

C(t)  =  -  Jz(y)  Z(y  +  t)  dy  (5~9) 

0 

where  L  is  the  sample  length. 
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A  model  found  to  fit  the  autocorrelation  function  for  many  fandom 
surfaces  is 

R(t)  =  exp  (-t/3*)  ( 5—10) 

This  has  a  correlation  length  of  3  =  2.3  6*  where  1/3*  is  the  decay 
rate  of  the  function  at  the  origin.  The  correlation  length  3  was 
originally  defined  oy  Whitehcuse  and  Archard  as  the  distance  over  which 
the  autocorrelation  function  R(t)  takes  to  decay  to  a  value  of  0.1. 

The  spatial  variation  can  also  be  represented  in  the  form  of  the 
structure  or  variance  function  (SF)  (Whitehouse  1971  ,  Savles  &  Thomas 
1977)  defined  for  a  surface  as 

S(r,9)  =  E  { (z (x,y)  -  z(<x,y)  +  r(e)))2}  (5-11) 

where  E  {  }  denotes  expectation.  Z(x,y)  is  the  surface  height  at 
coordinates  (x,y)  in  the  plane  of  the  surface  and  Z((x,y)  +  r(6))  is 
the  surface  height  at  a  radial  distance  (r)  from  (x,y)  in  a  direction 
(©) .  For  a  profile  SF  is  simply 

S(t)  =  E  {(z(y)  -  z (y  +  t))2}  (5-12) 

which  in  an  integral  form  is 

L 

s(t)  =  r|(z(y)  ■ z(y  +  T))Z  dy  (5_13) 

0 

vjhere  (y)  is  the  distance  along  the  profile  and  t  the  ordinate  sep¬ 
aration  in  the  (y)  direction. 


The  function  represents 


mean  sauare 


difference  in  height 


expected  over  any  spatial  distance  t . 


The  power  spectral  density  function  (PSDP)  or  power  spectrum 
represents  the  same  information  as  the  ACF  or  SP  but  in  a  different 
form.  The  PSDF  is  generally  given  in  terms  of  the  ACVF: 


G  (w) 


CO 

tf 


C(t)  cos  (wt)  dr 


(5-14) 


where  to  is  a  radial  frequency  2ir/X  and  X  is  the  wavelength.  The  PSDF 
can  also  be  obtained  directly  from  the  surface  profile  data  through 
the  equation 


00 

G (<ij)  =  7  Jz(y)e~jWy  dy 

—CO 


(5-15). 
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6.  THE  FLOW  ANALYSIS 


6.1  Introduction 


The  flow  channel  is  30.5  cm  wide  with  a  maximum  flow  depth  of  30  cm 
and  is  supported  at  three  points  including  two  jacking  points.  The 
jacking  facility  can  be  used,  along  with  a  control  valve  on  the  channel 
inlet  and  a  variable  height  weir  on  the  channel  outlet,  to  produce 
flows  at  varying  mean  velocity  and  mean  depth  in  subcritical  or  super¬ 
critical  regimes. 

The  regime  of  flow  in  a  rectuangular  channel  is  defined  by  the  use 
of  the  non-dimensional  Froude  number,  defined  by  F  =  U/»/gh  .  where 
U  is  the  mean  flow  velocity,  g  is  the  gravitational  acceleration,  and 
h  is  the  depth  of  flow.  From  the  dynamic  considerations  the  Froude 
number  may  be  shown  to  be  proportional  to  the  ratios  of  a  representa¬ 
tive  inertia  force  to  a  representative  gravity  force.  Hence  for  a 
liquid  of  constant  density  with  a  free  surface  the  effects  of  gravity 
on  the  fluid  inertia  are  expressed  by  the  Froude  number  F  where  the 
fluid  regime  is  defined  as  subcritical  for  F  less  than  unity,  and 
supercritical  for  F  greater  than  'unity. 

For  the  purpose  of  this  investigation  the  flow  is  maintained 
below  its  critical  value  in  order  to  minimise  flow  instabilities  and 
secondary  flow  effects  which  are  produced  by: 

(1)  Misalignments  of  the  glass  panels  and  the  steel  sections 
making  up  the  sides  and  the  floor  of  the  channel  respectively. 


(2)  Non-uniform  or  disturbed  flow  at  the  channel  inlet. 


(3)  The  introduction  of  flow  measurement  probes  into  the  flow 
regime. 

The  joints  between  the  side  panels  and  the  floor  sections  are  made 
as  unobtrusive  as  possibla  to  the  flow  by  careful  alignment  and  sub¬ 
sequent  smoothing  out  the  joints  with  a  filler. 

6.2  The  preliminary  flow  study 

6.2.1  The  flow  at  the  channel  inlet 


The  flow  system  was  initially  operated  with  the  flow  issuing  into 
the  channel  directly  from  the  deturbulating  tank  (Figure  2-1) .  An 
investigation  of  the  flow  conditions  at  various  stations  along  the 
channel  showed  the  existence  of  large  scale  eddies  which  were  producing 
transverse  instabilities . 

These  flow  instabilities  were  thought  to  originate  from  swirl 
produced  in  the  pipe  upstream  of  the  deturbulating  tank  in  combination 
with  the  multi-orifice  discharge  apparatus  through  a  blanked-off  tee 
located  near  the  bottom  of  the  deturbulating  tank.  These  were  sub¬ 
stantially  reduced  with  the  introduction  of  two  honeycomb  flow  straight 
eners  made  from  aluminium  sheets,  across  the  entire  section  of  the 
deturbulating  tank.  The  honeycomb  flow  straighteners  are  5  cm  thick 
with  a  honeycomb  diameter  of  approximately  1  cm.  One  of  the  honeycomb 
flow  straighteners  is  positioned  horizontally  in  the  deturbulating 
tank  40  cm  above  the  inlet  'pepper  pot'  inverted  T-piece  whilst  the 
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other  honeycomb  is  positioned  transversely  at  the  channel  inlet.  The 
net  effect  of  these  flow  straighter.ers  is  to  break  up  any  large  scale 
eddies  present  in  the  deturbulating  tank  and  produce  pseudo  two  dimen¬ 
sional  uniform  flow  at  the  channel  inlet. 

6.2.2  The  side  wall  shear  layers 


The  development  of  shear  layers  from  the  glass  channel  sides  will 
significantly  modify  the  flow  in  a  narrow  channel  by  producing  a  con¬ 
vergence  effect  on  the  free  stream  flow.  Although  a  convergence  term  * 
can  be  introduced  into  the  flow  analysis  it  was  found  that  the  magnitude 
of  the  term,  which  is  dependendent  on  measured  transverse  velocity 
profiles  at  various  stations  along  the  channel  length,  had  an  unaccept¬ 
ably  high  possible  error.  The  convergence  term  is  given  a  distance  Xc 
measured  from  the  channel  inlet  which  is  estimated  to  be  the  flow 
length  over  which  the  side  wall  shear  layers  develop  ar.d  eventually 
merge.  Thus,  to  graphically  estimate  the  distance  Xc  from  the  appro¬ 
priate  transverse  channel  velocity  profiles  requires  the  flow  to  be 
uniform  potential  flow  at  the  measurement  datum  (i.e.  at  the  channel 
inlet)  and  also  that  the  transverse  profiles  give  a  clear  measurable 
representation  of  the  state  of  the  side  wall  shear  layer  development. 
Figure  6-1  shows  a  schematic  representation  of  the  procedure  used  to 
interpolate  the  distance  Xc  from  transverse  channel  profiles.  Typical 
measured  transverse  profiles  are  shown  in  Figures  6-2,  6-4  and  6-5. 

6.2.3  The  flow  regime 

An  examination  of  the  flow  in  uhe  channel  using  the  transverse 
velocity  profile  data  from  the  hot-film  anemometer  (Figures  6-2,  6-4 

*  Proc.  Comput.  of  Turb.  Bound.  Layers,  Stanford  Univ.  1568 
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and  6-5)  showed  there  to  be  no  appreciable  growth  of  the  side  Wail 
shear  layers  in  a  2.4  ra  length  of  the  channel  chosen  to  be  the  test 
section  for  all  subsequent  measurements.  Hence  the  flow  analysis  used 
assumes  the  growth  of  the  shear  layers  over  the  channel's  glass  sides 
to  be  negligible.  Hence  from  4.5  n  downstream  of  the  channel  inlet  to 
6.9  m  from  the  inlet  a  series  of  boundary  layer  traverses  were  taken 
to  determine  the  variation  of  along  the  channel  centre  line,  for  a 
particular  work  surface  and  a  particular  set  of  flow  conditions. 


Thus  from  the  initial  investigation  of  the  flow  regime  in  the 
channel  a  test  section  was  chosen  such  that  the  flow  in  a  volume  2.4  m 
long,  30.5  cm  wide  and  with  a  height  bounded  by  the  work  surface  and 
the  open  surface  of  the  flow  is  assumed  to  be  substantially  an  acceler¬ 
ating  two-dimensional  flow  with  a  developing  work  surface  shear  layer 
and  relatively  unobtrusive  parallel  side  wall  shear  layers. 


6.3  The  analysis 


The  flow  analysis  is  divided  into  four  sections,  one  section  for 
each  of  the  four  surfaces  investigated. 


6.3.1  Flow  over  a  painted  steel  surface 


It  was  decided  initially  to  carry  out  a  study  of  the  flow  over  a 


work  surface  consisting  of  three  6.1  m  lengths  of  0.5  inch  thick  shot- 
blasted  steel  plate  butted  up  end-to-end  on  the  channel  floor.  The 
surface  finish  consists  of  a  brush-applied  coat  of  marine  paint*  over  a 
spray-applied  coat  of  primer  paint**. 


*  A  thixotropic  anticorrosion  marine  paint.  Silver  Primocon. 


**  A  two-pack  reinforced  edge-type  primer  paint,  Metagard  G250. 
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The  larger  capacity  punip  B  is  connected  into  the  channel  system 
and  the  smaller  axial  flow  pump  isolated.  The  radial  flow  pump  B  which 
operates  under  conditions  of  a  negative  static  suction  head  is  primed 
using  a  vacuum  line  connected  to  a  compressed  air  driven  air  ejector. 
Prior  to  tailing  any  flow  measurements  the  system  is  purged  of  air 
pockets  and  the  not-film  sensor  calibrated  (see  Chapter  3.2).  By 
arrangement  of  the  channel  slope,  weir  height  and  pump  outlet  valve 
setting  a  number  of  flow  conditions  are  set  up  in  the  channel.  The 
floiJ  analysis  assumes  the  flow  to  be  sub-critical  at  a  constant  depth, 
and  substantially  two-dimensional. 

Three  sets  of  flow  conditions  were  analysed  in  the  present  section, 
namely  flow  conditions  1,  2  and  3.  The  flow  parameters  and  values  of 
determined  from  measurements  along  the  centre  line  of  the  channel 
flow  control  volume  for  each  of  these  conditions  are  summarised  in 
tables  6-1  to  6-3.  Transverse  flow  profiles  for  the  flow  conditions 
are  shown  in  Figures  6-2  and  6-3.  The  transverse  profiles  show  a 
pseudo  free  stream  region  across  the  middle  two- thirds  of  the  flew 
where  the  measured  mean  flow  velocities  may  fall  off  to  90  per  cent  of 
the  maximum  velocity  recorded. 

The  range  of  Reynolds  number  over  which  the  flow  measurements 
could  usefully  be  taken  was  limited  chiefly  by  the  constraining  effect 
of  the  channel  sides,  causing  significant  departures  from  the  assumed 
two-dimensional  flow .  Gravity  surface  waves  were  also  encountered 
immediately  downstream  of  the  honeycomb  sheet  at  the  channel  inlet. 

The  amplitude  of  the  surface  waves,  for  a  flow  Reynolds  number  of 
5.3  x  106  and  a  depth  of  flow  of  24.5  cm  was  negligible  over  the  test 
section,  which  corresponded  to  the  maximum  Reynolds  number  attained 
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in  flow  3  (Table  6-3}  with  a  Broude  number  of  0.53.  A  typical  value  of 
5  over  the  test  section  for  flow  3  was  7  cm  'which  leaves  the  remaining 
two-thirds  or  more  of  the  flow  to  attenuate  any  surface  wave  effects  on 
the  boundary  layer  development.  At  Reynolds  numbers  much  greater  than 
5.3  x  10  the  flow  regime  begins  to  significantly  break  away  from  the 
simplified  model  which  the  flow  analysis  requires.  The  calculated 
values  of  the  local  coefficient  of  surface  friction  are  consistent  'with 
those  based  on  the  momentum  approach  and  are  approximately  equal  in 
magnitude  to  the  published  values  for  hydrodynamics lly  smooth  surfaces 
(Schlicnting,  'Boundary  Layer  Theory'). 

The  scatter  in  the  deduced  values  of  the  flow  parameters  and  the 
low  values  of  B,  the  shape  factor,  are  thought  to  be  a  consequence  of 
calibrating  the  flow  measurement  instrumentation  in  a  flow  environment 
dissimilar  to  the  one  in  which  it  is  to  be  operated.  The  hot-film 
probe  is  calibrated  on  the  centre  line  in  smooth  turbulent  pipe  flew 
at  pipe  Reynolds  numbers  between  3.7  x  10’  and  16.2  X  1G3  corresponding 
to  a  range  of  centre  line  velocities  between  approximately  25  cm/sec 
and  100  cm/ sec.  the  hot- film  probe  is  operated  in  open  channel  flow 
at  varying  distances  from  the  work  surface  and  therefore  varying 
turbulent  intensity  at  channel  Reynolds  numbers  ranging  from  2.3  /.  10** 
to  5.3  x  10**.  The  degree  of  turbulence  will  depend  on  the  roughness 
of  the  wall  and  the  distance  from  the  wall. 

The  hot-film  anemometer  is  assumed  to  have  a  linear  response  to 
values  of  mean  streaerwise  flow  ve  lei  ties  u  measured  at  points  ranging 
from  adjacent  to  the  channel  floor  to  the  free-stream  region  of  the 
channel  flow.  The  flow  velocities,  so  determined,  range  from  a  minimum 
of  approximately  10  cm/sec  within  J  ram  of  the  channel  floor  to  a 
maximum  of  approximately  80  cm/sec  in  the  free  stream. 


I’MlK  Kim 


The  determination  of  Cf  using  the  carter* turn  integral  equation 
requires  the  evaluation  of  thv?  two  integrals  for  momentum  and  displace¬ 
ment  thickness  over  the  boundary  layer  thickness  5  -  The  integrals 
are  inherently  very  sensitive  to  measurement  inaccuracies  in  the  inner 
regions  of  the  boundary  layer.  With  the  instrumentation  used,  errors 
in  the  determination  of  z  and  u  are  unfortunately  at  a  maximum  in  this 
f  critical  region  near  the  solid-liquid  interface. 

i  -  6.3.2.  Flow  over  a  gravel  surface 

The  second  surface  flew  regime  to  be  examined  consisted  of  flow 

i 

§ 

over  the  channel  floor  covered  With  a  closely  packed  layer ‘of  gravel. 

i 

To  reduce  any  movement  of  the  gravel  to  a  minimum  the  free  stream 
velocity  of  the  flew  was  reduced  by  isolating  the  larger  capacity  pump 
B  and  connecting  into  the  channel  system  the  smaller  capacity  axial 
flow  pump  A.  Prior  to  any  flow  measurements  the  hot-film  sensor  is 
re-calibrated  and  the  system  purged  of  any  air  pockets  by  throttling 
back  the  flow  at  the  channel  inlet  valve  and  opening  vents  in  the  over¬ 
head  pipeline,  the  channel  slope,  weir  height  and  channel  inlet  valve 
were  readjusted  to  produce  the  desired  flow  conditions  and  then  the 
boundary  layer  velocity  profile  sampled  at  5  stations  along  the  channel 
test  section.  The  subsequent  flow  analysis  is  summarised  in  Table  6-4. 

Transverse  flew  profiles  are  shown  in  Figure  6-4.  These  show 
I  a  pseudo  free  stream  region  across  the  middle  two-thirds  of  the  flow 

where  the  measured  mean  flow  velocity  may  fall  off  to  80  per  cent  of 
the  maximum  velocity  recorded. 

The  flow  parameters  are  determined  for  channel  Reynolds  nunbers. 
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Rc,  ranging  from  1.5  x  lO^  to  2.4  x  10*^  corresponding  to  free  stream 

velocities  of  36.7  cm/sec  to  38.8  cm/sec  respectively.  The  flow  was 

well  below  the  critical  Froude  value  and  there  were  no  detectable  surface 

waves.  The  flow  depth  was  15.5  cm  and  the  measured  boundary  layer 

velocities  ranged  from  a  minimum  of  1  cm/sec  adjacent  to  the  gravel 

surface  to  38.8  cm/sec  maximum  in  the  free  stream.  The  determined 

shape  factors  H  show  the  flow  to  be  separated  and  the  values  of  Cf 

from  the  momentum  integral  are  consistant  with  flow  over  a  sand  roughened 

plate  for  xA  =  1500  (Schlichting)  where  k  is  the  grain  size, 
s  s 

Hence  for  x  =  570  cm  this  corresponds  to  a  grain  size  of  approximately 
4  mm  which  is  consistent  with  the  surface  analysis  of  the  gravel 
(Chapter  7) . 


6.3.3  Flow  over  a  ship  hull  replica 


After  removing  the  gravel  from  the  channel  floor  the  4.8  m  long 
replica  of  the  hull  surface,  described  in  the  surface  analysis  Chapter 
7,  was  lifted  into  the  channel  and  positioned  on  the  channel  floor  with 
its  upstream  end  3.0  m  from  the  channel  inlet.  The  replica  fitted 
tightly  between  the  channel  sides  and  any  small  gaps  were  filled  with 
waterproof  filler.  An  aluminium  tapered  leading  edge  similar  to  one 
positioned  against  the  upstream  end  of  the  three  steel  plates  was  posi¬ 
tioned  at  the  upstream  end  of  the  replica,  the  tapered  upstream  insert 
caused  a  gradual  change  in  the  invert  level,  i.e.  between  the  channel 
floor  and  the  superposed  replica.  An  abrupt  change  would  undoubtedly 
cause  a  local  disturbance  with  perhaps  a  separation  bubble.  A  tapered 
section  was  also  positioned  downstream  of  the  hull  replica,  again  to 
obviate  any  abrupt  changes  in  the  flow  which  would  engulf  part  of  the 
flow  over  the  replica. 


If ito  Mu?  will  tm 
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The  replica  was  positioned  on  the  channel  floor  such  that  the  test 
section  was  sufficiently  downstream  of  the  replica's  leading  edge  and 
upstream  of  the  trailing  edge  to  eliminate  any  streamwise  acceleration 
or  deceleration  effects.  Transverse  flow  profiles  we re  taken  at  stations 
along  the  replicas  length  (Figure  6-5)  and  these  show  a  pseudo  free 
stream  region  across  the  middle  two-thirds  of  the  flow  where  the  measured 
mean  velocities  may  fall  off  to  90  per  cent  of  the  maximum  velocity 
recorded. 


As  with  the  steel  plate  surface  the  larger  capacity  radial  flow 
pump  B  can  be  used  for  the  flow  over  the  replica  surface  and  the  smaller 
axial  flow  pump  A  isolated.  With  the  possibility  of  precipitating 
disturbances  at  the  leading  edge  of  the  replica's  surface,  it  was  deemed 
appropriate  to  investigate  the  f lev;  conditions  immediately  downstream 
of  the  replica's  leading  edge  (Table  6-5). 


The  preliminary  study  of  the  flow  pattern  on  the  upstream  part  of 
the  replica's  surface  showed  no  measurable  detrimental  effects.  Hence 
with  the  chosen  flow  conditions  a  series  of  9  boundary  layer  profiles 
were  recorded  along  the  working  surface  test  section.  A  summary  of  the 
subsequent  flow  analysis  is  given  in  Table  6-6.  The  determined  values  of 
the  local  coefficient  of  surface  friction,  using  the  momentum  integral, 
approximate  to  published  values  for  flow  over  hydrodynamically  smooth 
surfaces  (Schlichting  'Boudary  Layer  Theory').  The  mean  value  of  the 
shape  factor,  H,  of  1.4  shows  the  boundary  layer  to  be  well  below  flow 
separation  and  the  mean  value  of  the  Froude  number  of  0.4  agrees  with 
the  tranquil  visual  appearance  of  the  flow. 


6.3.4  Flow  over  a  coated  ship  hall  replica 


The  final  surface/ flow  combination  studied  in  this  investigation 
used  the  4.8  m  long  hull  replica  with  its  surface  covered  with  a  brush 
applied  coat  of  thixotropic  anti-corrosion  marine  paint*  (See  Chapter 
7.6.1).  The  flow  system  was  operated  under  identical  conditions  to  the 
flow  over  the  uncoated  replica's  surface  and  a  series  of  transverse 
flow  profiles  are  recorded  (Figure  6-6) .  Boundary  layer  traverses  were 
taken  at  the  nine  stations  along  the  test  section  and  the  resulting 
flow  analysis  is  summarised  in  Table  6-7.  The  magnitudes  of  the 
values  determined  from  the  momentum  integral  approximate  to  those  for 
the  uncoated  replica  surface.  An  inspection  of  the  surface  analyses 
(Chapter  7)  for  these  two  surfaces  shows  there  to  be  a  relatively  small 
difference  in  the  surface  geometry  with  respect  to  the  overall  geometry 
of  the  surfaces. 

6.3.5  The  determination  of 

To  determine  the  local  values  of  the  coefficient  of  surface  fric¬ 
tion  it  is  necessary  to  determine  the  free  stream  velocity  and  momentum 
thickness  gradients  along  the  test  section  for  each  flow/surface  in¬ 
vestigated.  Hence,  best  fit  lines  are  drawn  through  the  momentum  thick¬ 
ness  and  free  stream  velocity  data  for  each  flow  system,  and  the  re¬ 
spective  gradients  used  in  the  analysis.  Plots  of  free  stream  velocity 
versus  streamwise  distance,  x,  and  momentum  thickness  versus  x,  are 
shown  in  Figures  6-7  to  6-10. 

The  error  bars  on  the  free  stream  velocity  plots  show  a  total 
possible  error  of  2.5  cm/sec  on  the  calculation  of  a  free  stream 

*  Silver  Primocon  marine  paint 
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velocity  value,  taken  as  the  maximum  difference  between  five  current 


meter  readings.  The  error  bars  on  the  momentum  thickness  plots  are  made 
up  of  the  errors  derived  from  the  integration-by-strips  method  of 
determining  0  and  the  error  derived  from  the  location  of  the  y  origin 
on  the  velocity  distribution.  Hence  each  valut  of  momentum  thickness, 
so  determined,  has  a  possible  error  of  0.05  cm. 

The  plots  force  the  experimental  data  onto  best  fit  straight  lines 
and  hence  make  the  assumption  that  both  the  free  stream  velocity  and 
momentum  thickness  vary  linearly  along  the  channel  length. 

Variations  in  the  boundary  layer  thickness  and  the  possibility  of 
deviations  from  two-dimensional  flow  will  produce  significant  errors  in 
the  determination  of  the  flow  parameters.  A  simplified  picture  of  the 
turbulence  intensity  across  the  measureu  boundary  layers  can  be  pro¬ 
duced  by  displaying  the  alternating  anemometer  bridge  voltage  from  the 
hot-film  probe  on  an  oscilloscope.  The  degree  of  turbulence  will 
depend  on  the  roughness  of  the  wall  and  the  distance  from  the  wall. 


An  oscilloscope  with  triggering,  storage  and  camera  attachment 
facilities  was  used  along  with  a  polariod  camera  to  record  the  altern¬ 
ating  component  of  the  CTA  bridge  voltage.  Hence  a  series  of  pictures 
was  taken  of  the  alternating  component  of  the  CTA  bridge  voltage  for 
the  hot-film  probe  positioned  at  points  across  the  working  surface 
shear  layer.  Th«  photographs  were  taken  with  the  hot-film  probe 
positioned  6.9  m  from  the  channel  inlet,  above  the  coated  replica  work 
surface .  The  oscilloscope  time  base  was  set  on  5  m  sec  and  the  verti¬ 
cal  magnification  on  0.1  volt/division.  Hence  Figures  6-1 la  to  6-llh 
show  the  low  'turbulence'  level  in  the  region  very  near  to  the  solid 
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surface,  the  increased  turbulence  level  further  away  from  the  wall  in 
the  turbulent  shear  layer,  and  the  lower  turbulence  intensity  in  the 
free  stream  region.  In  contrast  to  the  varying  degree  of  turbulence 
in  the  channel  flow,  the  hot- film  anemometer  probe  is  calibrated  in  a 
situation  where  the  intensity  of  turbulence  is  at  a  minimum  at  the 
centre  line  in  smooth  turbulent  pipe  flow.  Hence  there  will  be  an 
element  of  uncertainty  in  the  interpretation  of  the  anemometer  bridge 
voltage  as  the  probe  is  traversed  through  the  channel  floor  boundary 
layer.  This  will  add  to  uncertainty  in  the  calculation  of  Cf. 

The  location  of  the  origin  of  the  z  axis  on  the  work  surface  is  of 
special  interest  in  this  and  other  similar  analyses.  For  set  conditions 
of  flow  over  a  smooth  flat  surface,  the  plane  of  the  surface  can  be  used 
as  a  datum  from  which  the  boundary  layer  velocity  profile  can  be 
measured.  Since  all  points  on  the  surface  are  on  the  datum,  then  the 
flow  velocity  at  any  point  in  the  flow  will  be  related  to  the  distance 
z  measured  from  the  plane. 

For  a  rough  surface  a  best  fit  plane  can  be  fitted  through  the 
surface,  but  not  all  the  points  on  the  surface  will  be  on  the  plane. 
Hence,  if  the  plane  is  used  as  a  datum,  the  origin  of  the  datum  will 
not  always  coincide  with  points  on  the  surface.  The  choice  of  reference 
plane  must  somehow  be  related  to  the  interaction  between  the  fluid  flow 
and  the  rough  surface.  For  example,  a  surface  with  a  relatively  large 
value  of  mean  slope  may  have  steep  sided  peaks  between  which  the  fluid 
will  be  relatively  unaffected  by  the  streamwise  fluid  motion.  Pools 
of  fluid  will  occupy  the  valleys  between  these  steep  sided  peaks,  with 
the  bulk  of  the  fluid  shearing  over  the  tops  of  the  surface  peaks  and 
these  pools.  Hence,  a  suitable  datura  would  seem  to  be  definable  by  a 


flexible  plane  or  line  which  is  forced  to  lie  across  the  tops  of  the 
surface  peaks.  The  problem  then  is  to  decide  when  the  geometry  of  the 
surface  is  such  that  the  fluid  flow  near  the  surface  cannot  be  adequately 
represented  by  the  viscous  shear  approach.  The  scale  and  intensity  of 
turbulence  near  to  the  solid  surface,  the  values  of  the  local  surface 
slopes  and  curvatures  and  the  distance  between  the  surface  peaks  may 
well  be  important  in  deciding  the  shape  of  the  surface  datum. 
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FIG  6-2  MEAN  CTA  BRIDGE  VOLTAGE, 
NON-DIMENSIONLIZED  WITH  MAXIMUM  VOLTAGE, 
VERSUS  DISTANCE  Y  ACROSS  CHANNEL 
FDR  FLOW  OVER  PAINTED  STEEL  SURFACE 

(FLOW  1) 

each  plotted  point  represents  the  mean  and 
standard  deviation  of  readings  taken  at  2  cm. 
intervals  from  z=2cm.  to  the  fluid  free  surface 
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FIG  6-3  MEAN  CTA  BRIDGE  VOLTAGE. 
NON-DIMENSIONLIZED  WITH  MAXIMUM  VOLT 
VERSUS  DISTANCE  Y  ACROSS  CHANNEL 
FOR  FLOW  OVER  PAINTED  STEEL  SURFACE 

( FLOW  2) 
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:  6-5  MEAN  CTA  BRIDGE  VOLTAGE. 
N-DIHENSIONLIZED  WITH  MAXIMUM  VOI 


RSUS  DISTANCE  Y  ACROSS  CHANNEL 
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FIG  6-6  MEAN 

NGM-DIMhNSIOi''Liicu  wi  i  n  ham 
VERSUS  DISTANCE  Y  ACROSS  CHANNEL, 

FLOW  OVER  COATED  REPLICA  OF  HULL  SURFACE 
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ANALYSIS1  OVER  PAINTED  STEEL  SURFACE 
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FIG  6-7  PLOT  OF  FffEE-STRWrfEM 
VERSUS  DISTANCE  ALONG  FLOW  CHANNE 
FOR  REPLICA  OF  HULL  SURFACE 
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FIG  6-9  PLOT  OF  FREE-STREAM  VELOC1 
VERSUS  DISTANCE  ALONG  FLOW  CHANNI 
FOR  COATED  REPLICA  OF  HULL  SURFACE 
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FIG  6-10  PLOT  OF  MOMENTUM  THICKNESS 
VERSUS  DISTANCE  ALONG  FLOW  CHANNEL 
FOR  COATED  REPLICA  OF  HULL  SURFACE 
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In  this  investigation  tests  have  been  carried  out  on  various 
surface  geometries  which  make  up  the  floor  of  a  17,4  m  long  glass-sided 
flow  channel  (Figure  2-1).  An  attempt  is  made  to  relate  a  section  or 
sections  of  the  spectrum  of  surface  wavelengths  present  to  changes  in 
the  flow  regime  resulting  in  relatively  large  increases  in  the  surface 
frictional  resistance.  The  study  is  divided  into  two  main  sections: 
the  measurement  and  analysis  of  the  channel  floor  surface  geometry  and 
the  measurement  and  analysis  of  the  fluid  flow  over  the  channel  floor. 
The  results  of  both  sections  are  then  brought  together  and  the  inter¬ 
relation  of  the  surface  geometry  and  frictional  resistance  deduced. 

7.2  Preliminary  roughness  measurements  on  ship  hull  surface 
replicas _ 

Surface  measurements  were  taken  on  two  surface  replicas  and  the 
subsequent  statistical  analyses  of  the  surface  were  used  to  help 
identify  significant  characteristics  of  the  type  of  surface  under 
investigation.  The  first  of  these  replicated  surfaces  was  a  positive 
acrylic  replica  taken  from  a  heavily  corroded  ship  hul?  by  the  British 
Ship  P.esearch  Association  (BSRA)  . 

A  network  of  50  mm  long  surface  profiles  was  recorded  using  the 
Talylin  stylus  instrument  connected  to  a  data  logging  system.  The 
surface  measured  approximately  15  cm  by  30  cm.  A  specially  designed 
stylus  is  used  for  these  relatively  rough  surfaces,  the  stylus  having 
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an  included  angle  of  30  degrees  and  a  tip  width  of  460  pm.  The  electri¬ 
cal  output  from  the  stylus  instrument  is  digitised  and  recorded  on 
punched  paper  tape  for  computer  analysis,  the  results  of  which  are 
given  in  Tab].es  7-1  to  7-8  and  Figures  7-1  to  7-3.  The  analyses  show 
the  surface  peaks  to  have  a  much  higher  radius  of  curvature  than  the 
valleys  and  the  height  distributions  to  be  negatively  skewed.  These 
results  are  typical  of  a  profile  which  has  lost  its  peaks  in  some  way. 

The  most  common  explanation  for  the  removal  of  peaks  is  given  in  terms 

<e. 

of  some  kind  of  wear  mechanism  (Thomas  1972) .  The  peaks  will  be  worn 
do\m  towards  their  broader  bases  increasing  the  distance  across  the 
peak,  whilst  the  valley  will  be  relatively  unaffected.  The  peaks' 
higher  radius  of  curvature  and  the  negatively  skewed  height  distribution 
can  also  be  explained  by  the  modification  of  the  surface  by  the  repli¬ 
cation  process.  The  physical  properties  of  the  replication  material 
may  be  such  that  it  cannot  exactly  follow  the  contours  of  the  specimen 
surface,  bridging  valleys  or  building  up  on  the  shoulders  of  peaks. 

The  second  replj  ca  surface  consisted  of  a  positive  resin  replica 
taken  from  the  hull  of  a  Shell  oil  tanker.  The  amplitude  of  the 
surface  features  was  such  that  the  Talylin  instrument  using  both  the 
standard  stylus  and  the  specially  designed  stylus  was  incapable  of 
recording  the  large  height  variations  encountered  on  the  surface. 

Hence  an  extension  stylus  with  a  ball  end  of  1.6  mm  was  constructed 
producing  a  2:1  lever  reduction  on  the  standard  stylus.  The  maximum 
stylus  traverse  length  was  reduced  by  half  to  50  mm  which  coincided 
with  the  reference  length  used  by  the  BSRA  Wall  Gauge  (Canham  1956) . 

The  stylus  was  used  to  record  ten  50  mm  profiles  from  the  surface. 

The  measurements  consisted  of  two  parallel  traverses  along  the  length 
of  the  replica  each  made  up  of  five  consecutive  50  mm  sample  lengths. 
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A  summary  of  the  analyses  of  the  ten  sample  lengths  is  given  in 
Tables  7-9  and  7-10.  The  results  again  show  the  mean  peak  radius  of 
curvatures  to  be  greater  than  the  mean  valley  radius  of  curvatures, 
hence  the  peaks  are  more  rounded  than  che  valleys. 

The  measurement  of  replicated  surfaces  is  often  found  to  be  the 
only  practical  method  of  carrying  out  a  comprehensive  analysis  of  a 
surface,  as  has  been  the  case  for  ship  hulls,  whose  inaccessibility 
and  surface  condition  is  such  that  the  necessary  delicate  electronic 
surface  measuring  equipment  is  generally  inoperative  in  such  conditions 
The  need  for  measurements  on  actual  ship  hull  surfaces  is  underlined 
here  by  the  inconclusive  analyses  of  the  replicated  surfaces.  Hence 
new  surface  measuring  instruments  are  required  which  are  capable  of 
working  in  various  orientations,  unaffected  by  possible  contact  with 
salt  water  and  marine  growths,  and  capable  of  providing  the  information 
necessary  for  a  comprehensive  surface  analysis. 

7.3  Surface  measurements  on  detachable  ship  hull  plates 

Surface  measurements  were  carried  out  on  four  specially  painted 
steel  plates  before  and  after  their  attachment  to  the  bilge  keel  of  a 
supertanker  for  the  duration  of  a  6  month  deep  sea  voyage.  These 
measurements  we re  primarily  to  investigate  the  surface  properties  of 
a  specific  marine  paint  under  the  action  of  fluid  flow.  The  paint 
used,  produced  by  Intematicn  Marine  Coatings,  is  specified  as  a  'self¬ 
polishing'  marine  paint.  Under  the  action  of  fluid  flow  over  the 
coated  surface,  the  tops  of  the  highest  surface  peaks  are  'polished' 
down  whilst  the  valleys,  protected  by  adjacent  peaks,  are  relatively 
unaffected,  resulting  in  an  overall  reduction  of  the  surface  roughness. 
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The  specimen  plates  used  we re  100  cm  long  by  70  cm  wide.  Three  of 
the  plates  were  roughly  painted  to  simulate  a  typical  in-service  re¬ 
painting  whilst  the  fourth  plate  was  smoothly  painted.  Two  traverses 
using  a  Talylin  stylus  instrument  were  taken  on  each  plate  as  shown  in 
Figure  7-4.  Each  traverse  is  made  up  of  between  9  and  11  consecutive 
50  mm  gauge  lengths  which  are  analysed  separately.  The  standard  Talylin 
•stylus,  a  steel  chisel  with  a  radius  of  curvature  of  0.6  -  1.2  mm  could 
be  used  on  the  smoothest  plate  whilst  for  the  other  three  roughly  painted, 
plates  the  coarser  stylus  with  a  ball  end  of  1.6  mm  diameter  was  used 
(see  previous  surface  measurements,  Section  7.2).  A  summary  of  the 
subsequent  surface  analysis  is  given  in  Tables  7-11  to  7-14.  In  each 
column  the  first  figure  is  the  average  for  all  the  500  mm  gauge  lengths 
in  that  traverse;  the  second  figure  (in  parentheses)  is  the  standard 
error  of  the  mean;  the  third  figure  is  the  number  of  gauge  lengths 
averaged  over.  The  results  from  the  two  parallel  traverses  on  each 
plate  are  reasonably  self-consistent  on  all  plates  except  MBl.  This 
is  true  both  for  the  initial  and  later  results,  and  is  a  good  indica¬ 
tion  that  the  surfaces  are  homogenous. 

The  smoothest  plate,  MB2,  does  not  appear  to  become  any  smoother, 
in  fact  its  CLA  roughness  has  increased  (Table  7-11),  while  the  peak 
and  valley  radii  of  curvature  increased  by  between  lb  and  3b  times. 

The  figures  for  the  slopes  are  not  very  informative,  since  one  has 
increased  while  the  other  has  decreased.  The  other  three  plates 
(Tables  7-12  to  7-14)  become  smoother  by  a  consistent  70  to  80  per 
cent.  Each  of  the  three  roughness  parameters  CLA,  F.MS  and  MAA  show 
this  effect  to  nearly  the  same  degree.  Their  mean  slopes  are  lower  by 
70  to  80  per  cent  and  the  peak  and  valley  radii  of  curvature  are  higher 
by  factors  of  between  lb  and  8. 
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In  Figure  7-5  cumulative  surface  height  distributions  are  plotted 
on  a  normal  probability  scale  and  it  can  be  seen  that  the  earlier  dis¬ 
tribution  is  much  rougher  over  the  upper  portion  of  the  surface  than 
over  the  lower.  This  result  agrees  with  its  visual  appearance  as  a 
relatively  smooth  plane  covered  with  raised  droplets  of  dried  paint. 

The  later  distribution,  though  overall  much  smoother,  shows  a  similar 
trend,  suggesting  that  the  loss  of  roughness  is  not  necessarily  all  at 
the  peaks.  Other  evidence  which  suggests  this  is  that  the  increases 
in  peak  radii  of  curvature  are  no  greater  than  the  increases  in  the 
valley  radii  of  curvature.  The  computed  values  for  the  power  spectral 
density  functions,  derived  from  the  autocovariance  function  are  shown 
in  Figure  7-6.  The  figure  shows  the  drop-off  in  power  of  the  wave¬ 
lengths  present  on  one  of  the  roughly  painted  plates,  1JT1,  after  the 
ocean  voyage  in  comparison  to  before  the  voyage. 

The  above  results  show  that  the  marine  coating  used  has  the  useful 
hydrodynamic  property  of  its  surface  becoming  smoother  with  the  action 
of  fluid  motion.  Thus,  a  hull  surface  may  be  rapidly  recoated  in  dry 
dock  and  the  resulting  uneven  surface  finish  will  become  smoother  whilst 
the  ship  is  in  service.  This  is  a  much  more  economic  method  of  pro¬ 
viding  a  smooth  surface  finish  than  by  producing  the  same  surface  in 
dry  dock. 

The  afore-mentioned  results  from  the  coated  detachable  ship  hull 
plates  can  be  explained  by  the  following  hypothesis:  7;hen  the  roughly 
coated  surfaces  are  placed  in  situ  on  the  ships  hull,  the  highest  peaks 
of  the  coating  may  protrude  through  the  quiescent  inner  region  of  the 
ship’s  boundary  layer.  The  summits  of  the  highest  peaks  will  subse¬ 
quently  be  exposed  to  the  relatively  violent  action  of  the  outer  turbu- 
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lent  flow.  The  smaller  peaks  and  the  surface  depressions,  enveloped  by 
the  inner  viscous  layer,  are  protected  from  the  wearing  action  of  the 
turbulent  outer  flow  and  so  will  remain  unchanged.  The  properties  of 
the  coating  are  such  that  the  fluid  action  will  produce  a  surface  with 
an  overall  reduction  in  the  mean  value  of  roughness. 

This  implies  that  the  viscous  layer  adjacent  to  the  specimen  sur¬ 
face  has  a  thickness  comparable  to  the  surface  roughness  amplitude. 

In  fact,  the  final  value  of  the  average  or  rras  roughness  amplitude  will, 
if  the  polishing  mechanism  is  correctly  explained,  give  an  indication 
of  the  thickness  of  the  viscous  layer  adjacent  to  the  specimen  surface. 
The  mean  value  of  the  rms  roughness  of  the  surfaces  after  the  ocean 
voyage  is  approximately  15  |im,  hence  if  all  the  surface  features  pro¬ 
truding  through  the  viscous  layer  have  been  removed  then  the  thickness 
of  the  viscous  layer  is  in  the  region  of  15  |im. 

7.4  Introduction  to  the  measurement  of  the  channel  test  surface 

geometry  _ . 

The  initial  channel  floor  working  surface  consisted  of  three 
lengths  of  half  inch  thick  shot-blasted  steel  plate  butted  up  end  to 
end  along  the  flow  channel.  A  selection  of  surface  measuring  instru¬ 
ments  were  used  to  record  as  wide  a  spectrum  of  surface  features  as 
possible.  In  this  way  as  much  information  as  possible  about  the  surface 
geometry  can  be  made  available  for  an  analysis  of  the  work  surface. 

7.5  Surface  measurement  instrumentation 


The  measurement  procedure  was  such  that  the  longer  surface  fea- 
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tures  were  investigated  initially  and  then  instruments  capable  of 
recording  surface  features  of  decreasing  dimensions  were  employed. 


(a)  The  Talyvel*  electronic  level  is  used  to  record  profiles 
along  the  entire  length  of  the  channel  floor.  The  Talyvel  has  a  samp¬ 
ling  interval  of  a  few  centimetres  and  a  sampling  length  of  up  to  the 
total  length  of  the  work  surface.  The  instrument  measures  the  inclin¬ 
ation  of  a  surface,  upon  which  its  feet  are  resting,  by  means  of  a  sus¬ 
pended  pendulum  whose  flat  top  rests  against  two  variable  inductance 
transducers.  The  transducers  provide  an  electrical  displacement  signal, 
related  to  the  tilt  of  the  pendulum,  which  is  amplified  and  applied  to 
a  centre-zero  meter.  The  instrument  is  traversed  along  the  channel 


centre  line  in  steps  equal  to  the  mean  separation  of  the  feet,  9  cm. 
The  inclination  of  the  surface  is  measured  in  units  of  inches  of  in¬ 


clination  per  inch  and  a  small  computer  program  is  used  to  convert 
these  into  height  readings. 


(b)  A  dial  gauge,  held  in  a  magnetic  clamp  with  its  measuring 
probe  perpendicular  to  the  work  surface,  is  used  to  record  spot  heights 

r 

on  the  surface  relative  to  a  122  cm  long  straight  edge  laid  lengthwise 
along  the  centre  line  of  a  section  of  the  work  surface.  The  dial  gauge 
is  used  to  record  surface  features  of  horizontal  dimensions  from  a  few 
millimetres  to  tens  of  centimetres  in  length. 


The  dial  gauge  has  a  ball  end  measuring  head  with  a  diameter  of 
2  mm.  One  small  division  on  the  dial  scale  represents  a  vertical  move¬ 
ment  of  2.54  f*m. 


*  Talyvel  is  a  Rank  Taylor  Hobson  instrument. 
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(c)  Measurements  can  also  be  taken  along  the  entire  .length  of  the 


work  surface  using  a  mechanical  stylus  instrument,  the  BSRA  wall  gauge. 
The  gauge  consists  of  a  metal  stylus,  in  contact  with  the  surface, 
connected  by  mechanical  linkage  to  a  smaller  stylus.  The  smaller  stylus 
duplicates  the  motion  of  the  larger  stylus  by  scratching  its  path  onto 
a  graphite  coated  glass  slide,  the  length  of  the  recorded  profile  is 
reduced  by  a  factor  of  10,  whilst  the  transverse  motion  remains  at 
unity.  The  motion  of  the  stylus  is  relative  to  a  track-way  magneti¬ 
cally  clamped  to  the  work  surface  upon  which  the  gauge,  mounted  on 
wheels,  is  manually  propelled. 

The  recorded  surface  profile  can  be  analysed  by  projecting  the 
magnified  image  of  scratched  profile  onto  a  digitising  table  connected 
to  a  paper  tape  punch.  The  sampled  height  readings  can  then  be 
computer  analysed  to  give  a  measure  of  the  surface  roughness. 

(d)  A  Ferranti  three-axis  measuring  table  is  used  to  measure 
spot  height  profiles  for  surfaces  which  have  a  variation  in  heights 
too  great  to  be  recorded  by  the  afore-mentioned  instruments.  The 
tapered  steel  measuring  probe  has  a  tip  radius  of  approximately  500  pm. 
The  probe  is  held  in  a  chuck  along  the  vertical  axis  of  the  measuring 
table.  The  probe  can  then  be  traversed  along  the  three  axes  manually 
to  any  position  on  the  work  surface  placed  on  the  measuring  table  below 
the  probe.  A  series  of  traverses  can  be  taken  across  the  work  surface, 
for  which  the  probe  is  moved  horizontally  in  convenient  step  lengths 
and  the  probe  lowered  to  come  into  contact  with  the  surface  at  each 
measuring  point.  The  orientation  of  the  probe  relative  to  the  throe 
axes  can  be  read  off  from  three  digital  displays  which  are  connected 
to  a  print-out  unit  which  automatically  records  the  relative  position 
of  the  probe  on  punched  paper  tape. 


(e)  Three  electronic  stylus  instruments*  were  available  to 
measure  the  relatively  small  surface  wavelengths  in  the  range  from  a 
few  micrometres  to  a  few  centimetres  in  length.  A  typical  stylus 
instrument  (Figure  7-7)  consists  of  a  pick-up  driven  by  a  gear  box 
which  draws  the  stylus  across  the  surface  at  a  constant  speed.  The 
stylus  is  connected  to  a  transducer  which  converts  its  vertical  movement, 
relative  to  some  datum,  into  an  oscillating  electrical  signal  which  is 
then  electronically  amplified.  The  output  can  then  be  either  fed  to  a 
chart  recorder  for  a  profile  trace  to  be  made,  or  it  can  be  digitised 
and  stored  on  punched  paper  tape  for  a  subsequent  computer  analysis. 

The  Talysurf  and  Surf com  stylus  instruments  are  such  that  for  this 
investigation  they  can  not  be  used  in  situ  on  the  channel  floor  and 
hence  a  representative  piece  of  the  channel  floor  surface  had  to  be 
brought  to  the  instruments  for  measurements  to  be  taken. 

7.6  The  measurement  and  analysis  of  the  surfaces  under  investi- 
gation _ _ 

7.6.1  A  preliminary  study  of  the  means  available  to  alter  the 
geometry  of  a  surface _ 

In  order  to  correlate  surface  parameters  with  related  flow  para¬ 
meters  it  is  essential  to  have  a  means  of  producing  a  variety  of  work 
surfaces  with  predictable  trends  in  their  overall  geometry.  Applying 
successive  coats  of  paint  to  a  surface  was  though  to  be  a  possible  way 
of  producing  a  number  of  surfaces  with  differing  surface  features. 


*  The  Talylin  1  and  the  Talysurf  3  stylus  instruments  produced  by  Rank 
Taylor  Hobson  and  the  Surf com  produced  by  Ferranti. 
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Hence  a  study  was  carried  out  to  examine  the  changes  in  the  stat¬ 
istical  microgeornetry  of  a  steel  surface  'with  successive  coats  of  paint, 
using  a  stylus  instrument.  Preliminary  measurements  on  the  coating 
material*  showed  the  stylus  instrument  to  be  capable  of  recording 
accurate  profiles  of  the  coating  surface.  The  subsequent  measurements 
showed  that  on  the  application  of  the  initial  coat  of  paint  the  surface 
roughness  decreased  substantially,  whereas  on  further  application  of 
paint  the  surface  roughness  does  not  change  by  anywhere  near  the  initial 
amount. 


The  results  are  given  in  full  in  Appendix  3,  as  published  in  the 
Journal  of  Coatings  Technology,  August  1978. 


7.6.2  Looking  into  the  possibility  of  damage  caused  by  stylus 
instruments  on  coated  surfaces 


The  diamond  stylus  of  a  Talysurf  3  instrument  was  drawn  across 
the  surface  of  a  coated**  glass  slide  and  a  Vickers  high  power  micro¬ 


scope  used  to  view  the  area  of  the  surface  traversed  by  the  stylus. 

A  mark  left  by  the  stylus  could  be  detected  and  an  approximate  value 
of  its  depth  was  obtained  using  the  microscopes  calibrated  micrometer 
scale.  By  focussing  on  the  top  and  bottom  of  the  stylus  mark  the 
depth  of  which  was  found  to  be  approximately  1  pm. 


A  light-section  microscope  set  on  x400  magnification  is  used  to 


view  the  section  of  surface  over  which  the  stylus  has  passed.  Any 


*  Silver  Primocon  marine  paint. 

**  The  coating  consisted  of  a  dip-applied  coating  of  Silver  Primocon 
paint,  applied  twenty  four  hours  previous  to  the  investigation. 


length  of  10.2  cm  is  user!  with  a  sampling  interval  of  190  pm  for  each 
of  the  recorded  profiles.  A  summary  of  the  statistical  analysis  for 
the  surface  is  given  in  Table  7-15. 

A  small  sample  of  the  painted  steel  surface,  approximately  15  cm:  , 
is  produced  so  that  it  can  be  fixed  to  the  bed  of  the  Talysurf  3  stylus 
instrument  for  subsequent  measurement.  A  single  sample  is  taken  from 
the  downstream  end  of  the  steel  plates  due  to  the  destructive  nature  of 
this  measurement .  The  measurement  on  the  sample  is  assumed  to  be  rep¬ 
resentative  of  the  surface  as  a  whole  and  this  is  to  some  extent  justi¬ 
fied  by  the  homogenous  appearance  of  the  painted  steel  surface.  The 
Talysurf  stylus,  in  the  form  of  a  truncated  diamond  pyramid  of  tip 
dimensions  8  pm  by  2.5  pm,  is  traversed  across  the  surface  for  a  sample 
length  of  8  mm.  The  analysis  of  the  measured  profile  is  summarised  in 
Table  7-15. 

Figure  7-10  shows  a  single  plot  of  the  power  spectral  density 
functions  computed  from  the  measurements  made  by  the  various  surface 
instruments  used  on  the  painted  steel  surface.  The  curve  is  a  reason¬ 
able  fit  to  the  inverse-square  law  for  a  surface  with  a  gaussian  height 
distribution. 

The  analyses  as  a  whole  (Table  7-15)  show  the  measured  peak  and 
valley  curvatures  and  surface  slopes  to  decrease  with  increasing  samp¬ 
ling  interval  and  sample  length.  This  is  because  unlike  true  random 
data  which  nas  limits  of  plus  and  minus  infinity,  surface  heights  are 
distributed  within  a  relatively  narrow  band  governed  by  the  mechanical 
properties  of  the  constituent  material  and  the  surface  function. 
Similarly  the  RMS  roughness  increases  with  increasing  sampling  interval 


and  sample  length. 


7.6.4  The  gravel  surface 


A  specimen  of  the  gravel  surface  was  mace  uo  by  bonding  a  closely 
packed  layer  of  gravel  to  a  rigid  backing  surface  producing  a  surface 
jO  cm  .  It  was  not  posssible  to  record  the  relatively  large  surface 
slopes  present  on  the  gravel  surface  with  the  conventional  surface 


measurement  instrumentation.  Hence,  the  specimen  surface  was  placed  on 
a  Ferranti  three-axis  measuring  table  ana  a  profile  Of  the  surface  re¬ 
corded  on  punched  paper  tape.  The  sample  length  used  is  27  cm  with  a 
sampling  interval  of  1  mm.  A  summary  of  the  statistical  analysis  of 
the  recorded  profile  is  given  in  Table  7-17.  The  recorded  profile  is 
snovm  in  Figure  7-11.  The  distribution  of  surface  slopes  is  shown  in 
Figure  7-13  and  the  distribution  of  curvatures  in  Figure  7-12. 


7.6.5  The  ship  hull  renlica  surface  and  the  coated  replica  surface 


The  replica  of  the  ship  hull  surface  was  produced  using  a  two- 
stage  replication  process  perfected  by  Shell  Research  Ltd. 


The  technique  involves  clamping  a  framework,  to  retain  the  repli¬ 
cation  material,  against  a  vertical  flat  section  of  the  dry  docked  ship 
hull.  Neoprene  rubber  solution  is  poured  into  the  space  between  the 
hull  and  the  framework  and  cures  whilst  constrained  to  take  up  the 
geometry  of  the  hull  surface.  The  surface  replicated  was  taken  from 
the  side  of  a  badly  corroded  hull. 


The  resulting  'negative*  replica  of  the  surface,  approximately 
3  cm  thick  by  5  m  long  by  0.4  m  wide,  is  peeled  away  from  the  mating 
hull  surface  after  removing  the  constraining  framework.  The  first 
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stage  of  the  replication  process  can  then  be  transported  from  the  dry 
dock  to  the  laboratory  containing  the  flow  analysis  equipment.  A  flat, 
levelled,  wooden  base  was  made  up  on  the  laboratory  floor  and  the  neo¬ 
prene  sheet,  replicated  surface  upwards,  placed  on  it.  A  closely 
fitting  framework  was  bolted  around  the  edges  of  the  neoprene  replica 
to  form  the  sides  of  an  open-topped  reservoir  into  which  a  layer  of 
resin*  is  poured.  The  resulting  resin  cast  (Figure  7-14),  approximately 
1.5  cm  thick  by  483  cm  long  by  30.5  cm  wide ,  is  allowed  to  harden  for 

7  days  and  then  carefully  pulled  away  from  the  neoprene  negative. 

The  underside  of  the  resin  cast  forms  the  positive  replica  of  the 
hull  surface,  hence  the  up- turned  cast  can  be  lowered  into  the  flow 
channel  to  form  the  work  surface.  The  mass  of  the  replica  plus  its 
flexibility  allowed  it  to  take  the  form  of  the  channel  floor,  producing 
a  rigid  surface. 

A  statistical  analysis  of  the  replica  surface  is  computed  for  each 
of  four  instruments  capable  of  recording  a  wide  range  of  the  surface 
features.  The  hull  replica  is  investigated  both  for  its  untreated 
surface  and  with  a  coat  of  paint  applied  to  the  surface  (see  Chapter 
7.6.1).  Hence  for  each  surface,  digitised  profiles  were  recorded  using: 

(1)  A  Talyvel  electronic  level  is  used  to  record  3  parallel  pro¬ 
files  along  the  entire  length  of  the  surface.  Ths  traverses  are  taken 

8  cm  and  17  cm  from  the  front  channel  side  (see  Figure  2-1)  and  8  cm 
from  the  rear  side  of  the  channel.  The  recorded  profiles  ar*5  shown 

in  Figure  7-15  for  the  virgin  replica  surface  and  the  analyses  tor  the 
virgin  replica  surface  and  the  coated  replica  surface  sr©  summarised 
in  Tables  7-18  and  7-19. 

*  Shell  Chemicals:  Epikote  816  with  Episurc  KTV. 
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(2)  The  dial  gauge  plus  datum  are  used  to  record  three  118  cm  long 
profiles  spaced  out  along  the  length  of  the  work  surface ,  moving  to¬ 
wards  the  channel  inlet.  The  first  profile  started  at  the  outlet  end 
and  15.9  cm  from  the  front  side  of  the  flow  channel.  The  second  profile 
started  192.7  cm  from  the  outlet  end  of  the  replica  and  11.1  cm  from 
the  front  side  of  the  channel.  The  third  profile  was  taken  360.7  cm 
from  the  outlet  end  of  the  replica  and  9.5  cm  from  the  front  side  of 
the  flow  channel. 

The  recorded  profiles  for  the  untreated  replica  surface  are  shown 
in  Figure  7-16  and  a  summary  of  the  analyses  for  the  untreated  replica 
surface  and  the  coated  replica  surface  is  given  in  Tables  7-18  and  7-19. 

(3)  The  Talylin  instrument  is  used  to  record  surface  features 
measuring  from  a  few  millimetres  in  length  to  a  few  centimetres.  The 
standard  Talylin  stylus  cannot  faithfully  reproduce  the  relatively 
rough  surface  Oj.  the  replica.  Although  the  standard  stylus  could  not 
be  used,  the  specially  constructed  extension  stylus  (Chapter  7.2)  could 
be  sue  issfully  traversed  across  a  specimen  piece  of  the  replica  approx¬ 
imately  15  cm2  .  Prior  to  the  measurements  the  specimen  surface  was 
securely  mounted  on  a  relocation  table  (Williamson  &  Hunt  1968) .  The 
relocation  table  allows  repeat  profiles  to  be  recorded  along  the  same 
section  through  a  surface  and  hence  record  the  changes  in  the  surface 
with  the  application  of  paint  (Figure  7-17) .  The  paint  is  applied  by 
brush,  hence  profiles  were  also  recorded  across  and  along  the  visible 
brush  marks  (Table  7-20) .  The  results  of  the  analysis  of  the  reloca¬ 
ted  profiles  are  summarised  in  Table  7-21. 
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(4)  The  Surfcom  stylus  instrument  is  used  with  the  relocation 
table  to  record  relocated  profiles  of  the  changes  in  geometry  of  a 
section  through  the  replica  surface  with  the  application  of  paint 
(Figure  7-18).  Two  profiles  are  also  records  1  along  and  across  the 
brush  marks  left  in  the  coated  surface.  The  results  of  the  analysis 
for  the  relocated  profiles  are  summarised  in  Table  7-22. 


The  measuring  head  of  the  Surfcom  consists  of  a  diamond  stylus 
with  a  3  Jim  tip  and  a  loading  of  0.5  gram.  The  instrument  is  used  with 
a  sampling  interval  of  100  Jim  and  a  sample  length  of  2.5  cm.  The 
vertical  magnification  is  on  its  lowest  setting  of  x200  and  the  horiz¬ 
ontal  magnification  is  xlO. 


The  combined  analyses  for  the  replica  surface  and  coated  replica 
surface  data  are  shown  in  Tables  7-18  and  7-19  respectively.  The 
results  show  the  RMS  roughness  of  the  coated  surface  to  have  increased 
by  approximately  5  per  cent  for  sample  lengths  of  2.5  mm  and  10.2  cm. 
Since  a  reliable  and  rapid  measure  of  the  change  in  roughness  of  a 
surface  can  only  be  gained  from  repeat  measurements  through  the  same 
section  of  the  surface,  the  measurement  results  for  sample  lengths  of 
1.18  m  and  4.83  m  can  be  neglected  in  a  comparison  of  the  surface 
roughness.  The  measurements  in  fact  show  a  RMS  roughness  increase  of 
10  per  cent  and  a  decrease  of  10  per  cent  fc  the  sample  lengths  of 
1.18  m  and  4.83  m  respectively. 


The  values  of  the  measured  slopes  and  curvatures  for  each  surface 
decrease  with  increasing  sampling  interval  and  sample  length,  with  the 
exception  of  the  2.5  mm  sample  measurements.  The  maximum  vertical 
movement  of  the  2.5  mm  cut-off  instrument  is  such  that  a  relatively 
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smooth  section  of  the  surface  had  to  be  chosen  for  a  complete  profile 
to  be  successfully  recorded,  hence  the  artificially  3ow  values  for  the 
surface  slopes. 

The  distributions  of  profile  slopes  and  curvatures  are  shown  in 
Figures  7~i9  and  7-20  respectively.  The  distributions  are  for  a  sample 
length  of  5  cm  and  the  solid  lines  show  the  uncoated  replica  surface 
whilst  the  broken  lines  show  the  coated  surface.  There  is  little  diff¬ 
erence  between  the  distributions  for  the  coated  and  uncoated  surface. 


TABLE  7-3  SURFACE  ANALYSIS  OF  REPLICA  FROM  CORRODED  HULL 
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TABLE  7-11  SURFACE  ANALYSIS  OF  DETACHABLE  HULL  PLATE  MB2 
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TABLE  7-16  ANALYSIS  Of  SURFACE  RECORDS  FROM 
TRAVERSING  THE  B.S.RA  WALL  GAUGE  OVER  THE 
PAINTED  STEEL  SURFACE 


record  number  Mean  Apparent  Amplitude 


A 

1 

MAAljjm) 

56 

2 

69 

17 

71 

18 

69 

19 

83 

20 

60 

21 

27 

22 

39 

23 

47 

24 

30 

55  mean 

19  standard  deviation 


sample  length  -50mm. 

*  — * 

ball-ended  probe diameter 
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SURFACE  ANALYSIS  FOR  REPLICA  WITH 
SUCCESSIVE  COATS  OF  MARINE  PAINT 


8.  THE  SURFACE  ANALYSIS /FLOW  ANALYSIS  CORRELATION 


The  experimental  flow  data  for  the  test  surfaces  is  plotted  using 
the  inner  lav;  for  the  boundary  layer  velocity  distribution  (Figure  8-1), 
The  plot  shows  the  existence  of  the  logarithmic  region  in  the  mean 
velocity  distributions.  The  gravel  surface  data  is  displaced  downwards 
by  an  amount  Au/ut ,  dependent  on  the  roughness  Reynolds  number  (Musker 
ewkowicz  1978) .  The  painted  steel  surface  and  replica  surfaces 
-■  t- ■  oth  lie  on  the  well  known  universal  curve  for  smooth  walls 
(Clauser  1954) .  The  values  of  determined  for  the  paint  steel  and 
replica  surfaces  agree  with  the  inner  law  plot,  but  the  surface  analyses 
show  the  replica  surfaces  to  have  R^  values  for  a  sample  length  of 
1.14  m,  of  approximately  10  times  that  of  the  relatively  smooth  painted 
steel  surface  (Table  7-15) .  The  curvatures  and  slopes  for  the  replica 
profile  are  approximately  20  times  the  magnitude  of  the  curvatures  and 
slopes  for  the  painted  steel  surface. 

The  gravel  surface  has  a  CLA  roughness  of  approximately  equal 

magnitude  to  the  replica  surface  for  comparable  sample  lengths.  The 

mean  peak  curvature  for  the  gravel  profile  measurement  is  70  times  that 

of  the  replica  profile  for  a  sample  length  of  1.14  m,  whilst  the  grand 

pro."  -le  mean  valiey  curvature  and  profile  slopes  are  40  times  and  10 

times  respectively  the  values  for  the  replica.  The  local  coefficients 

of  surface  friction,  0.,  determined  for  the  gravel  surface  are  approx- 

f 

imately  2  times  the  magnitude  of  the  values  for  rhe  replica  surface 
at  comparable  Reynolds  numbers.  Therefore  the  difference  in  the  values 
of  the  local  coefficient  of  s  sface  friction  for  the  gravel  and  ship 
hull  replica  surfaces  is  clearly  bound  wirh  the  magnitude  of  the  surface 


curvatures  and  slopes,  ir.  particular  to  the  peak  curvatures. 


The  profile-  slope  distribution  for  the  gravel  surface  (Figure  7-13) 
has  a  kurtosis  of  2.44  whereas  the  replica  surface  and  coated  replica 
surface  have  profile  slope  distributions  of  kurtoses  5.32  and  6.08 
respectively  (Figure  7-20) . 

The  profile  slope  distribution  obtained  for  the  gravel  surface 
shows  there  to  be  a  relatively  large  proportion  of  high  angled  slopes 
present  on  the  surface.  The  profile  curvature  distributions  for  the 
replica  surfaces  show  there  to  be  a  relatively  large  proportion  of  small 
curvature  slopes  present  on  the  surfaces.  Hence  the  gravel  surface  has 
a  relatively  large  proportion  of  high  a'.gled  slopes  and  high  surface 
curvatures  with  respect  to  the  replica's  surfaces. 

As  can  be  seen,  this  analysis  uses  a  very  simplified  approach  to 
the  determination  of  flow  parameters  in  a  not  so  simply  defined  flow 
regime.  Added  to  this  are  the  complex  specimen  surfaces,  such  as  the 
replica  surface  with  its  relatively  long  surface  wavelengths,  'which  will 
affect  the  boundary  layer  development  and  hence  the  flow  parameters  in 
a  random  manner. 

The  results  of  the  flow  and  surface  analysis  for  the  painted  steel 
surface  are  consistent  with  published  results  for  hydrodynamically 
smooth  surfaces.  The  gravel  surface  gives  results  comparable  to  pub¬ 
lished  sand-grain  surface  friction  results. 

The  replica  surfaces  have  values  of  the  local  coefficients  of 
surface  friction  which  are  comparable  with  those  for  a  smooth  surface. 
The  replica  surface  analyses  show  them  to  have  CLA  values  comparable 


with  that  .‘or  the  gravel  surface,  but  with  shallower  slopes  and  smaller 
surface  curvatures. 

Whilst  it  is  feasible  that  the  magnitude  of  the  local  coefficients 
of  surface  friction  for  the  replica's  surface  is  less  than  that  for  the 
gravel  surface,  it  is  not  possible  to  justify  the  result  that  the 
replica's  surface  has  C,  values  comparable  to  a  hydrodynamically  smooth 
surface. 

The  inability  of  the  flow  analysis  to  detect  the  effect  of  the 
smaller  changes  in  the  geometries  of  the  surfaces  on  the  values  is 
attributed  to  an  over-simplification  of  the  flow  in  the  open-channel. 
The  single  sensor  hot-film  anercometry  used  in  the  investigation  is 
found  to  give  reasonably  reproducable  values  for  the  local  coefficients 
of  surface  friction  along  the  centre  line  of  the  test  section  surfaces. 

The  instrumentation  and  techniques  available  for  the  measurement 
and  analysis  of  the  test  surfaces  are  such  that  the  surfaces  can  be 
analysed  to  a  much  greater  degree  of  accuracy  than  the  analysis  of  the 
corresponding  fluid  motion.  Hence,  to  obtain  a  realistic  correlation 
between  fluid  parameters  and  surface  parameters,  the  flow  analysis 
must  be  refined  to  obtain  a  fuller  picture  of  the  fluid/surface  inter¬ 
action.  In  or^er  to  achieve  this  a  more  sensitive  and  reliable  flow 
measuring  instrument,  in  the  form  of  a  laser  doppler-shift  instrument, 
is  to  be  incorporated  into  the  flow  system. 

Thus  the  present  work  has  shown  the  feasibility  of  using  open- 
flume  flow  in  an  investigation  of  the  surface  frictional  resistance  of 
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nominally  flat  surfaces.  The  results  are  plotted  using  the  semi-empiri¬ 
cal  formula  in  terms  of  the  centre-line  average  roughness  height,  CL A. 


XCLA  5-62  1O910 


7.752  ku_/v 

_  T _ 

“0.005  ku  /v  _  „  ,  . 

e  t  +  0.44  ku  /v 

r 


(Musker  et  al  1976) 


where  k  is  the  CLA  roughness  parameter.  Figure  8-2  shows  the  results 
from  the  painted  steel,  gravel  and  replica  surfaces  plotted  in  terms  of 
the  roughness  function,  xCL^»  versus  the  roughness  Reynolds  number,  R^. 
The  plot  shows  good  agreement  with  the  results  of  Musker  et  al  (1976). 
The  results  are  also  plotted  using  the  modified  roughness  Reynolds 
number, 


kuT  =  0UT  (l  +  aS  ) (1  +  bS.  K  )  (Musker  &  Lewkowiez  1S78) 
~  ~  P  k  u 


Figure  8-3  shows  the  results  plotted  along  with  data  from  Musker 
&  Lewkowiez  (1978) .  The  results  from  the  present  work,  for  the  painted 
steel  and  replica  surfaces,  correlate  to  some  extent  with  the  published 
data,  whilst  the  gravel  surface  data  lies  well  clear  of  a  mean  curve 
fitted  by  eye  through  the  published  data.  It  has  been  found  that  no 
reasonable  curve  fit  could  be  produced  for  a  surface  sample  length 
higher  than  10  mm  regardless  of  the  values  assigned  to  the  empirical 


constants  a  and  b  (Musker  &  Lewkowiez  1978).  Hence,  Figure  8-3  is 
included  merely  to  give  an  indication  of  any  trend  in  the  present  data. 
The  local  surface  friction  coefficient  was  also  calculated  using  the 
surface  friction  law  by  Ludwieg  and  Tillman  (1950) , 


=  0.246.10  °-678H  (Re)-0'268 
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FIG  8-1  INNER  LAW  PLOT  USING  THE  FLOW  ANALYSES 
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Introduction 


■me  resistance  or  a  snip  m  staaav  motion 


cairn  water  can 


resolved  into  the  following  principal  components :  (a)  Surface-friction¬ 


al  resistance;  (b)  Wavemaking  resistance;  {c} 


and  {d}  Hind  and  air  resistance. 


The  shin  is  generally  assumed  to  be  in  steady  notion  in  order  to 


simplify  the  resistance  analysis. 


The  surface  friction  is  made  up  of  the  summation  of  the  components 
in  the  direction  of  motion  of  the  tangential  stresses  which  the  water 
produces  on  the  hull  surface,  and  apart  from  the  viscosity  and  density 
of  the  water,  it  depends  primarily  on  the  area,  length  and  roughness  of 
the  wetted  surface  of  the  hull,  as  well  as  the  speed  of  the  vessel. 


The  wavemaking  resistance  consists  of  the  summation  of  the  compo¬ 
nents,  in  the  direction  of  motion,  of  the  normal  pressures  on  the  hull 
due  to  the  surface  wave  systems  set  up  by  the  action  of  the  hull  on  the 
water.  Wavemaking  resistance,  apart  from  the  density  of  the  water,  is 
primarily  effected  by  geometrical  features  such  as  proportions,  fine¬ 
ness  and  shape  of  hull,  as  well  as  the  strength  of  the  gravity  field. 


The  third  Deponent ,  that  is,  resistance  due  to  eddying  wakes,  is 
associated  with  energy  losses  due  to,  firstly,  the  formation  of  eddies 
behind  the  stern  of  the  ship;  more  especially  these  associated  with  the 
rudder  and  stem  post  and  the  propeller  shaft  supports.  Under  certain 


i 


. . . 1 . 1 . . 


form  of  normal 


as 

ll 


i 


s 

i 


The  fourth  component  is  the  resistance  due  to  the  relative  wind 
acting  on  the  above-water  part  of  the  ship.  Due  tc  the  low  density  of 
air  compared  with  that  of  water,  this  resistance  is  generally  quite 
small  compared  with  the  water  resistance. 

Therefore  it  can  be  said  that  there  are  two  principal  components 
of  resistance,  namely,  surface  friction  and  wave-making  (Figure  fti-1) . 


Of  the  two  principal  components ,  the  surface  friction  resistance 
is  predominant  at  low  Fronde  numbers,  that  is  fox  relatively  low  speeds 
such  as  for  tankers.  At  higher  speeds  the  wave-making  component  be¬ 
comes  predominant. 


1.2  Ship-Modal  correlation 


The  inability  to  calculate  these  resistances  directly  from  hydro- 
dynamic  theory  led  to  the  use  of  tests  on  models  geometrically  similar 
to  the  shic. 


Typically  a  ship  will  have  a  resistance  comprising  of  up  tc  85 
per  cent  frictional  (Schoenherr  1932) .  Methods  of  measuring  the  re¬ 
sistance  of  ships  were  originally  based  on  the  determination  of  the 
power  required  to  propel  a  ship  at  a  certain  speed  and  then  comparison 
with  past  results. 


o 


n  more  convenient  and  much  used  method  in  the  early  1900*6  was  the 
Froude  Method.  This  assumed  that  the  frictional  component  of  the  total 
resistance  could  be  determined  independently  and  subtracted  and  that 
the  remaining  residuary  resistance  is  free  from  viscous  effects  and  is 
subject  to  the  lav/s  of  dynamical  similitude,  The  residuary  resistance 
is  related  to  inertia  and  gravitational  forces  only,  and  determinable 
by  ship-model  experiments  using  Froude1 s  lav?  of  dynamical  similarity. 

i.2.1  Flat  plate,  analysis 

in  1872  Froude  carried  out  a  series  of  experiments  to  investigate 
the  influence  oh  the  surface  friction  of  the  velocity  and  the  dimensions 
of  a  flat  plate  towed  edgewise  through  still  water,  in  particular  the 
length  of  the  plate  in  the  direction  of  motion.  The  first  known  tests 
on  the  surface  friction  of  flat  surfaces  in  fluids  were  those  of  Beau- 
foy  in  1793.  In  1880  Froude  ahd  Tideman,  using  the  experimental  data- 
of  Froude 's  earlier  towing  experiments ,  proposed  a  formula  of  the  form 
F  =  CAUm  for  the  frictional  resistance  of  plates.  That  is,  the  surface 
friction  was  assumed  to  vary  jointly  as  the  wetted  surface  area  h,  and 
the  speed,  U,  raised  to  some  power  in. 

Thus,  in  the  analysis  of  the  surface  friction  of  a  ship,  it  was 
found  convenient  initially  to  consider  the  surface  friction  of  a  flat 
surface  moving  in  its  own  plane.  It  was  found  that  the  above  relation 
held  quite  well  for  a  given  plate  but  that  C  and  m  varied  from  one  plate 


to  another. 


1 ; 2 . 2  Friction  lines 

When  the  experimental  results  of  Froude  ware  re-analysed  on  the 
more  acceptable  basis  of  Reynolds  number  and  due  account  taken  of 
laminar  flow  on  the  smaller  planes  (Lackenby  1937),  the  results  could 
then  be  fitted  with  a  universal  curve,  or  friction  line  (Figure  Al-2) . 

With  the  theory  of  the  boundary  layer  as  proposed  by  Prandtl  in  1904  it 
was  now  possible  to  define  a  characteristic  length  for  the  flow  along 
a  plane  surface  which  corresponded  to  the  radius  in  the  case  of  de¬ 
veloped  flow  in  pipes.  This  characteristic  length  was  called  the 
"boundary  layer  thickness".  Thus  with  the  flow  in  cylindrical  pipes 
and  the  flow  along  plane  surfaces  being  reduced  to  a  common  base.  Von 
Karman  and  Prandtl  showed  independently  how  experimental  data  obtained 
by  measuring  pressure  losses  in  smooth  pipes  could  be  used  to  deduce  a 
formula  for  the  frictional  resistance  of  smooth  plane  surfaces.  (The 

i 

investigations  were  extended  to  rough  pipes  by  Nikuradse  who  success¬ 
fully  correlated  the  frictional  data  of  dense  sand-grain  roughnesses 
of  different  sizes  in  various  diameter  pipes.)  This  formula  and  other  ; 

I 

! 

subsequently  produced  formulae  were  experimentally  and  not  theoretically 

i 

derived  and  hence  they  were  no  more  accurate  than  formulae  based  on 
direct  measurements  of  the  resistance  of  plates  towed  through  still 
water  or  suspended  in  a  wind  tunnel.  This  led  back  to  various  experi¬ 
mental  work  on  towing  plates  through  water  as  Froude  had  done  in  1872. 

The  first  investigation  of  the  resistance  of  plane  surfaces  in  air  was 
carried  out  by  Zahm,  in  1903. 

Each  of  these  individual  experiments  covered  only  a  limited  range 
of  length  and  speed  and  hence  produced  a  series  of  equations  for  the 
frictional  resistance  of  plane  surfaces  deduced  from  mean  lines  drawn 
through  each  of  the  plots  of  results. 
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Schocnhorr  (1932)  collected  together  ali  the  existing  «ota  of 
oxporimontal  data  for  the  frictional  resistances  of  smooth/  flat 
surfaces,  each  sot  being  for  a  limited  range  of  length  and  spoed. 
Experiments  wore  carried  out  at  low  Reynolds  numbers  (2  x  105  to  2  x  106) 
to  provide  new  data  where  existing  data  was  widely  scattered.  The 
kinetic  viscosity  was  varied  in  addition  to  the  length  and  speed.  The 
whole  of  the  existing  data  was  then  used  to  produce  a  single  plot  of 
surface  friction  coefficient  versus  Reynolds  number. 

Bearing  in  mind  that  these  tests  were  made  at  different  times/  in 
different  media  and  with  widely  different  apparatus,  the  plot  showed 
that  on  the  whole  the  results  were  in  fairly  good  agreement,  and  that 
an  average  curve  would  be  representative  for  all  the  tests  (Figure  Al-3) . 
Accordingly  the  best  representative  mean  curve  was  fitted  to  the  data. 

It  was  found  that,  on  attempting  to  fit  an  equation  to  this  mean  curve, 
none  of. the  existing  formulae  fitted  over  the  entire  range.  From  an 
extension  of  the  work  done  by  Von  Karman,  an  equation  of  the  form 


K/i^cJ  =  log  (RC^)  +  M  Schoenherr  (1932) 

was  applied  to  the  mean  line  through  the  data. 

It  was  found  that  with  the  following  values  for  the  constants  K 
and  M,  M  =  0  and  K  -  0.558,  the  equation  coincided  almost  exactly  with 
the  mean  line  over  the  whole  range  of  Reynolds  numbers. 


5 


1.2.3  Shin  hull .resistance 

The  extension  of  the  flat  plate  analysis  to  the  frictional  resis¬ 
tance  of  the  uneven  and  Curved  surfaces  of  ships*  hulls  has  been  attemp¬ 
ted  by  various  workers.  Several  of  these  set  out  to  determine  the  in¬ 
fluence  on  drag  of  different  types  of  roughness  to  be  found  on  a  typical 
ships  hull,  and  also  to  produce  a  parameter  for  the  effect  of  curvature 
of  a  shj.ps  surface.  But,  in  order  to  be  of  value  these  attempts  had  to 
be  checked  against  actual  frictional  resistance  of  the  ship  itself. 
Hence,  it  was  necessary  to  have  a  method  of  measuring  the  frictional 
resistance  of  ships  to  a  fair  degree  of  accuracy.  For  this  purpose,  it 
was  proposed  to  tow  ships,  as  had  been  done  previously  by  Froude,  at 
low  speeds  where  the  wave-making  resistance  was  considered  small.  Such 
tests  however,  are  expensive  and  the  results  are  of  limited  value . 

Hence  from  Froude *s  method  of  calculating  ship  resistance  which  is 
denoted  by  the  equation 


r/au2  -  R/yAU2  =  f/au2  -  F/  AU2 

where  ship  units  are  denoted  by  capital  letters  and  ship-model  units  by 
small  letters;  total  resistance  =  r,R;  frictional  resistance  =  f,F; 
area  of  wetted  surface  =  a, A;  velocity  =  u.U;  specific  gravity  of  salt 
water  =  y 

Hence  the  frictional  resistance  of  the  ship  can  be  determined  when 
the  total  resistance  is  known. 


It  was  shown  that  the  relation  between  the  resistance  of  a  ship 
and  the  propeller  thrust  to  overcome  this  resistance  involves  the  so' 


called  thrust  deduction  coefficient,  which  could  be  approximately  ob¬ 
tained  from  the  self-propelled  tests  of  the  ship's  model.  Methods  had 
also  been  developed  for  determining  the  wind  resistance  of  ships’ 
superstructure.  Prom  these  considerations  the  total  water  resistance 
of  a  ship  could  be  obtained. 

It  was  usually  assumed  that  the  frictional  resistance  of  the  model 
is  equal  to  that  of  a  smooth  plane  surface  of  equal  length  and  area,  and 
the  frictional  resistance  of  the  ship  is  equal  to  that  of  a  smooth  plane 
surface  of  equal  length  and  area  increased  by  an  allowance  for  the  re¬ 
sistance  of  the  butt  straps,  laps  and  other  projections  on  the  ship's 
hull.  Although  the  correlation  of  the  resistance  of  ship  models  of 
different  sizes  and  the  reliable  prediction  of  actual  ship  resistance 
from  model  results  is  of  fundamental  importance,  it  was  not  possible  to 
accurately  predict  the  frictional  resistance  of  a  ship  model  from  the 
results  of  another  model  to  a  different  scale. 

Later  work  by  Froude  on  plates  coated  with  fine,  medium  and  coarse 
sand  was  re-analysed  by  Telfer  (1951)  giving 


where  ks  is  the  sand  grain  size  and  L  the  ship  length. 

1.2.4  Geometrical  similarity  between, ships  and  their  models 

The  chief  criticism  made  against  the  use  bf  the  Schoen'nerr  ( 1932) 
curve  fit  (friction  line)  was  not  so  much  against  the  line  itself  but 
that  the  data  producing  the  line  shows  such  a  large  spread  of  values 


so  invalidating  the  use  of  a  mean  line,  mating  prediction  by  extrapo¬ 
lation  to  ship  values  uncertain.  The  range  of  variation  of  value  in 
the  Schoenhefr  plot  was  of  the  order  of  20  per  cant  in  the  range  of 
Reynolds  numbers  from  on  million  to  twenty  million.  Hughes  (1952) 
decided  to  examine  critically  the  basis  of  the  existing  data  and  to 


determine  the  possible  reasons  for  these  large  variations.  Differences 
in  tlie  experimental  conditions  were  thought  to  be  the  most  likely  cause 
of  these  variations.  It  was  found  from  this  examination  that  there  was 
a  promise  of  a  high  degree  of  consistency  (and  hence  conformity  to 
Reynold's  Law)  if  proper  account  was  taken  of  variations  in  the  ratio 
of  length  to  breadth  of  the  plane  surfaces  used  by  the  various  experi¬ 
menters,  or  in  other  words  if  proper  respect  ’were  given  to  the  condition 
of  geometrical  similarity  implicit  in  any  dynamic  similarity  lav;.  Hence, 
the  lack  of  geometrical  similarity  appeared  to  be  one  of  the  major  causes 
of  'apparent'  inconsistencies  revealed  by  the  survey  of  the  data,  as 
verified  by  Figure  Al-4.  Other  major  causes  appeared  to  be  due  to  the 
additional  effect  of  wave-making  resistance,  pressure  resistance,  and 
transitional  flow  phenomena. 


Thus,  before  any  new  data  was  obtained  it  was  decided  that  firstly 
flat  surfaces  having  a  range  of  values  of  the  length  to  breadth  propor¬ 
tions  would  be  used.  Secondly,  the  plates  were  made  as  thin  as  possible 
so  that  both  form  resistance  and  wave-making  resistance  would  be  mini¬ 
mised.  This  had  to  be  balanced  against  preserving  the  flatness  of  the 
surface,  the  flatness  becoming  more  difficult  to  achieve  as  the  thick¬ 
ness  of  the  plate  is  decreased. 


It  was  also  realised  that  it  was  essential  that,  the  flow  ever  the 
test  surface  must  be  fully  turbulent  from  as  near  the  leading  edge  as 
possible  and  at  the  lower  speeds  (Figura  Al-5) . 


8 


With  the  increased  interest  in  the  prediction  of  the  resistance  of 
ship  forms  from  tests  on  models  a  review  of  the  subject  was  carried  out 
by  Granville  in  1956.  It  became  apparent  that  when  plotting  the  total- 
resistance  coefficients  of  a  family  of  geometrically  similar  vessels, 
it  is  necessary  to  specify  non-dimensional  ship  lengths  to  satisfy 
dynamical  similarity,  h  graphical  method  of  analysing  the  resistance 
of  geometrically  similar  forms  was  developed  in  terms  of  flat-plate 
resistance. 

On  examination  of  the  existing  experimental  data  for  families  of 
geometrically  similar  models  it  v?as  found  that  there  was  insufficient 
accuracy  to  establish  definitely  the  superiority  of  any  one  prediction 
method.  Furthermore,  a  lack  of  tests  on  hydrodynamicaily  smooth  full- 
scale  vessels  prevents  verification  of  full-scale  predictions. 

1.2i5  Ship-Model  flow  regime 


For  drag  measurements  in  general  to  be  ’useful  it  is  necessary  to 
have  a  means  of  determining  the  type  of  flow  present,  whether  it  be 
laminar,  transitional,  or  fully  turbulent.  In  the  case  of  ship  models 
it  had  been  established  that  some  of  the  models  have  a  laminar  region 
of  flow  during  testing  (Figure  Ai-6) .  It  was  decided  that  ships  models 
for  reasons  of  comparing  results,  should  be  operated  with  fully  turb¬ 
ulent  boundary  layers.  The  techniques  that  became  available  for  exam¬ 
ining  the  state  of  the  boundary'  layer  flow  were  in  general  cumbersome, 
requiring  special  and  often  very  extensive  tests,  in  addition  to  the 


resistance  measurements. 


In  i960  Tovmsin  used  a  technique  for  flow  detection  which  did  not 
obstruct  the  flow  in  any  way,  the  sensing  elements  forming  part  of  the 
model  surface;  hence  simultaneous  flow  observations  and  resistance 
measurements  could  be  made.  The  method  used  was  similar  to  hot  wire 
anemometry,  but  instead  of  vising  wires,  which  are  easily  broken  in  water 
and  need  a  probe  to  hold  them  in  position,  it  used  platignum  hot  film 
sensors. 


Flow  measurements  were  made  over  the  forward  end  of  three  models 
of  different  lengths,  two  of  which  with  turbulence  stimulators  in  the 
form  of  studs.  The  third  one  was  a  9  foot  model  of  the  ship  Lucy  Ashton 
tested  both  with  a  trip  wire  as  a  stimulator  and  with  no  stimulator  at 
all.  Studs  are  the  generally  accepted  form  of  turbulence  stimulator. 


It  was  found  that,  even  at  the  top  end  of  the  speed  range,  a  lam¬ 
inar  region  persisted  behind  the  studs  on  the  largest  of  the  models 
(15.4  feet  in  length)  covering  half  to  one  per  cent  of  the  wetted 
surface.  The  smaller  the  length  to  breadth  ratio  of  the  vessels  tested, 
the  more  prone  they  were  to  laminar  flow. 


The  model  of  the  Lucy  Ashton  (1953) ,  which  has  a  relatively  large 
length  to  breadth  ratio,  was  found  to  be  free  of  laminar  flow  even  at 
the  lowest  speeds  when  the  trip  wire  was  used.  Without  the  trip  wire, 
a  laminar  region  covering  some  5  per- cent  of  the  model  surfaces  per¬ 
sisted  up  to  the  higher  speeds  used.  The  results  showed  conclusively 
the  need  to  study  the  flow  regime  during  the  testing  of  ship  models. 
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■scale  snip  trial: 


Ship  trials  at  sea  are  used  to  evaluate  the  accuracy  with  which 
prediction  of  the  full-scale  resistance  of  vessels  can  be  Bade  from 


tests  or 


sis.  Thrust  or  toraue  measurements  on  the  propeller  shaft; 


are  related  to  ship  resistance  on  the  basis  of  additional  tests  on 
models.  The  usual  result  from  sea  trials  was  that  the  resistance  so 
determined  is  greater  than  the  predicted  resistance,  the  difference 
between  them,  the  correlation  coefficient,  being  made  up  of  the  hull 
roughness  and  errors  due  to  scale  effects  in  the  propeller  allowances 
and  errors  in  the  prediction  method. 


The  most  recent  extrapolation  methods  for  ship-model  tests  (ITTC, 
September  1978)  are  good  for  most  classes  of  vessel  ana  the  correlation 
coefficient  is  now  almost  entirely  due  to  hull  roughness.  Differences 
between  the  prediction  method  and  the  ship  trial  results  give  rise  to 
the  following  roughness  allowance 


ACf  x  10*  105  Oc  /L) 


,,i/3 


where  k„_ _  is  the  BSEA  mean  apparent  amplitude  (HAA) . 


To  eliminate  the  uncertainties  due  to  the  action  of  the  propellers, 
more  difficult  and  costlier  tests  were  performed  on  hulls  without  pro¬ 
pellers,  the  propulsion  being  achieved  by  a  towing  vessel  or  more 
recently  by  aircraft  engines.  _ The  classical  test  of  this  kind  is  that 
of  the  Greyhound  by  William  Froude.  More  recently  3SRA  carried  out  an 
investigation  in  which  a  former  Clyde  steamer,  the  ’Lucy  Ashton’,  {1953; 
was  propelled  by  aircraft  jet  engines.  Tne  combined  measured  thrusts 
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ana  estIsati.cE  surface  friction  between  using  a  red  oxide  paint  and  sharp 
sears,  and  a  slightly  smoother  aluminium  paint  with  faired  seams 
(Figure  AI-8) .  On  allowing  the  hull  to  foul,  the  surface  resistance 
increased  by  about  50  per  cent.  The  fouling  consisted  of  scattered 
small  barnacles  about  1/8  inch  high  on  the  hull  bottom,  together  with  a 
band  of  fine  grasses  on  the  sides.  Thus,  the  trials  showed  the  surface 
friction  to  be  very'  sensitive  to  small  roughnesses,  and  more  so  to 
fouling  (Figure  Al-9). 


1.2.7  Analysis  of  the  ship  hull  surface 

In  an  attempt  to  threw  more  light  c  the  'Lucy  Ashton*  results, 
records  of  the  hull  surfaces  were  taken  with  a  wall  gauge  on  the  two 
painted  surfaces  used.  Analysi  c  these  showed  that  the  nature  of  the 
roughness  was  a  complex  function  ..  superimposed  undulations,  the  mean 
amplitude  depending  on  the  wavelength,  and  there  seemed  to  be  no  simple 
parameter  which  could  be  used  to  define  the  surface.  A  form  of  harmonic 
analysis  was  used  which  showed  that  the  amplitude  decreased  with  de¬ 
creasing  wavelength  and  that  over  most  of  the  range  the  aluminium  paint 
surface  was  smoother,  the  amplitudes  being  about  two-thirds  those  of  the 
red  oxide  paint  'which  agreed  with  the  measured  resistances.  From  this 
work  B5RA  realised  the  importance  of  the  effect  of  hull  roughnesses  on 
the  performance  of  ships  and  decided  to  make  hull  roughness  surveys  on 
all  subsequent  ship  trials  in  the  hope  of  explaining  differences  in 
performance  between  sister  ships:  there  could  be  as  much  as  23  per  cent 


variation  in  the  performance  f: 
analysis  was  used  to  look  at  tJ 


css  tor  SisteJ*  snxos.  «  scata. 


separata  factors  affecting  trial 


performance  and  this  snowed,  as  expected#  that  plate  roughness  is  by 
far  the  mast  significant  factor^  The  greater  part  of  the  variation  it 
the  hull  surface  roughness  seemed  to  be  due  to  the  uneven  application 
of  oaint. 


In  order  to  look  at  the  effect  of  deterioration  of  null  surface 
due  to  breakdown  of  the  paint  system,  formation  of  corrosion  products 
and  organic  fouling,  3SRA  in  conjunction  with  Shell  'Tankers  carried  out 
a  series  of  repeat  trials  on  two  sister  ships  after  various  periods  in 
service  (Figure  Al-10) .  The  dry  dock  inspections  showed  that  there  was 
little  or  no  organic  fouling  ana  that  the  deterioration  of  the  hull 
surface  was  primarily  due  to  corrosion.  The  plot  of  measured  shaft 
horse  powers  against  corresponding  hull  roughness  (Figure  11-11)  shows 
considerable  scatter  but  has  a  marked  trend  of  shaft  horse  power  in¬ 
creasing  with  hull  surface  roughness. 


1.3  Ship  hull  roughness/ fluid  flow  considerations 


1.3.1  Tests  on  the  effects  of  different  types  oi  structural  hull 
roughness,  painting,  corrosion  and  fouling,  on  fluid  flows 


In  an  attempt  to  determine  the  influence  of  the  various  types  of 
roughness  likely  to  occur  on  a  ship's  hull,  Kespf  (1937)  ran  three 
comparative  sets  of  tests  on  pontoons  whose  flat  bottoms  were  covered 
with  the  types  of  surface  under  examination .  The  first  set  of  tests 
comprised  two  different  conditions  of  paint- roughness  ana  investigated 
the  influence  of  butt  weld  joints.  The  second  set  of  tests  was  to 
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Kenpf  then  went  on  to  measure  the  additional  roughness-resistance 
of  ships  using  two  new  approaches.  In  the  first  base  F rouse ' s  method 
of  measuring  the  total  resistance  of  a  towed  ship  and  of  her  model  was 
adopted  by  using  the  Prandtl-Schlichting  formula  for  minimum  frictional 
resistance  for  the  model  to  calculate  the  coefficient  of  frictional 
resistance  for  the  ship.  The  second  method  was  based  on  the  measure¬ 
ment  of  the  local  resistance  of  a  certain  small  piece  of  the  hull 
surface  from  the  paralle  middle  body  part  of  the  ship  (Ksmpf  1529) , 
and  provided  a  direct  measure  of  local  frictional  resistance. 


The  results  produced  C^,  values  for  the  ships  which  were  too  high 
at  the  lower  Reynolds  number,  i.e.  at  lower  speeds.  This  discrepancy 
was  explained  as  the  result  of  there  being  laminar  flow  round  the 
models  which  would  make  the  measured  model  resistance  at  lower  speeds 
too  small,  and  as  the  turbulent  frictional  resistance  is  subtracted, 
the  residual  resistance  for  the  ship  will  be  calculated  too  small. 


Kerar-f  then  went  on  to  look  at  the  effect  of  fouling  on  the  re¬ 


sistance  of  shies.  Included  In  this  category:  (1)  rusting,  (2)  grdwt 
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of  shell,  (3)  growth  of  seaweed.  He  formulated  the  following  conclu¬ 
sions  regarding  the  above  aspects: 

(1)  Rusting  is  a  gradual  process  over  the  years,  affecting  the 
whole  bottom  of  the  ship.  The  equivalent  sand  roughness  may  be  doubled 
or  more  by  rusting  producing  10  per  cent  to  20  per  cent  more  frictional 
resistance. 

(2)  The  resistance  of  ship  hull  surfaces  fouled  with  shells  can 
be  evaluated  using  an  equivalent  sand  grain  roughness,  where  the  shells 
are  assumed  to  be  similar  in  shape  to  the  sand  grains.  A  comparison 
was  made  of  the  experimental  data  obtained  for  a  clean  ship  hull  bottom 
and  for  one  covered  with  a  growth  of  shells.  The  shells  were  of 
'medium'  height  (approximately  4.5  mm)  with  a  density  on  the  surface 
producing  a  maximum  resistance  to  the  flow.  The  results  showed  that 
frictional  resistance  of  the  shell  coated  surface  could  be  as  much  as 
85  per  cent  more  than  the  clean  surface. 

(3)  The  resistance  of  seaweed  which  mostly  covers  a  small  strip 
under  the  water-line  was  not  tested,  but  from  ship  trials  it  seemed 
that  smooth  seaweed  did  not  affect  the  resistance  very  much  in  com¬ 
parison  to  the  sharp-edged  roughness  of  shells. 

Telfer  (1969)  considered  the  problem  of  fouling  in  the  case  of 
tankers  and  found  that  owing  to  their  quick  turn-round  in  port  they  do 
not  usually  suffer  badly  from  fouling,  but  in  the  event  of  fouling 
being  present  the  chip's  frictional  resistance  could  be  substantially 
increased. 
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At  the  present  time  toxic  anti-fouling  is  applied  to  the  hull  in 
the  form  of  a  paint.  The  anti-fouling  system  must  adhere  well  to  the 
underlying  anti-corrosive  coatings  or  else  tie  coats  and  barrier  coats 
may  be  required.  Removal  of  the  anti-fouling  coat  due  to  it  becoming 
ineffective  can  seriously  increase  the  surface  roughness,  resulting  in 
an  expensive  blasting  procedure  to  renovate  the  surface.  In  practice 
anti-fouling  coatings  have  to  be  renewed  frequently,  at  say,  one  to 
two  year  intervals. 

Allan  and  Cutland  (1956)  looked  at  the  problem  of  the  effect  of 
hull  roughness  on  the  ship  resistance,  by  studying  separately  the 
various  structural  roughnesses  occurring  on  a  ship's  hull.  Full  scale 
experiments  were  carried  out  on  pontoons  with  the  various  roughnesses 
on  plates  fixed  to  their  undersides.  The  results  from  these  tests  were 
then  applied  to  a  number  of  representative  ships  in  the  large  cargo 
ship  range  (450  feet  to  620  feet  long) . 

It  was  found  that  the  local  velocities  r>t  the  edge  of  the  boundary 
layer  calculated  in  the  fore-body  of  the  ship  were  above  that  for  the 
flat  surface,  whereas  those  in  the  after-body  were  below  that  for  the 
flat  surface.  Thus,  it  was  thought  probable  that  the  effect  of  rough¬ 
ness  in  the  front  quarter  length  of  the  ship  was  substantially  under¬ 
estimated,  whereas  the  opposite  was  probably  happening  for  the  last 
quarter  length,  so  the  two  effects  tended  to  balance  out.  Contrary  to 
belief,  the  middle  quarters  of  the  ship's  length  were  found  found  to 
be  as  important  as  the  leading  quarter  in  their  contribution  to  the 


overall  resistance. 


In  conclusion  it  was  thought  that  even  when  care  is  taken  to  ensure 
that  the  local  small-scale  roughness  does  hot  exceed  a  few  thousandths 
of  an  inch,  the  effect  of  resistance  will  still  be  considerable  at  ship 
Reynolds  number;  The  effect  on  resistance  of  roughness  arising  from 
bad  paintwork  was  shown  to  increase  the  surface  friction  by  as  much  as 
100  per  cent,  thus  again  showing  the  need  to  ensure  a  smooth  hull 
finish. 

Results  derived  from  measurements  of  ship  hull  roughness  and  per- 
formance  on  ship  trials  collected  over  a  period  of  years  by  the  BSRA 
(1977)  showed  that  a  one  per  cent  increase  in  power  is  required  to  main¬ 
tain  speed  for  every  10  pm  (0.0004  inch)  increase  in  the  mean  apparent 
amplitude  (Chaplin  1965)  of  the  hull  roughness.  Hence,  the  bad  surface 
in  Figure  Al-12  would  be  expected  to  require  17  per  cent  more  power 
than  the  average  surface,  whilst  the  good  surface  'would  require  8  per 
cent  less  power  than  the  average  hull  surface. 

A  typical  new  ship  has  a  hull  roughness  equivalent  to  130  pm 
(0.005  in)  mean  apparent  amplitude,  whereas  an  old  ship  in  poor  condi¬ 
tion  can  exceed  1,000  pm  (0.04  in).  Whilst  it  is  thought  possible  that 
improvements  of  100  pm  (0.004  in)  and  more  could  be  made  for  older 
ships,  the  finish  on  a  newly  built  ship  is  unlikely  to  be  better  than 
80  pm  (0.003  in). 

A  study  by  Townsin  et  al.  (1976)  found  that  the  condition  of  a 
ship's  bottom  and  the  effect  it  has  on  the  ship's  performance  could  be 
particularly  important  in  the  economic  operation  of  the  ships  studied. 
Weld  beads  could  contribute  to  the  increase  in  power  by  1  per  cent  or 
more.  Hence,  buffing  off  the  butt  welds  seemed  in  some  cases  to  be 


economically  attractive.  On  several  ships  a  J-rge  increase  in  the  hull 
roughness  was  observed  after  recoating  in  dry  dock.  The  popular  method 
of  recoating  by  rapid  application  of  the  paint  using  a  high  pressure 
spray  was  found  to  give  an  uneven  'over-sprayed'  surface  which  could  be 
rougher  that  the  uncoated  surface. 


Frequent  underwater  scrubbing  of  the  ship  hull  to  remove  marine 
growths  such  as  weeds  gave  'worthwhile  short  term  improvements  in  the 
ship  performance,  but  fouling  should  be  eliminated  completely  if  large 
increases  in  fuel  consumption  are  to  be  avoided  (Townsin  1978) . 

1.3.2  The  initial  preparation  of  the  hull  surface 

Hudson  et  al.  (1959)  studied  the  effect  of  surface  preparation 
on  the  performance  of  the  subsequently  painted  surface  of  ships’  bottoms 
in  an  attempt  to  provide  a  common  procedure  on  surface  preparation  for 
the  shipbuilding  and  shipping  industries  as  a  whole.  The  experimental 
work  consisted  of  producing  specimens  of  ship  plate,  in  both  the  as- 
rolled  and  the  descaled  conditions,  exposing  them  to  'weathering  for 
various  periods  of  up  to  12  months,  the  weathered  plates  were  then 
treated  by  several  accepted  practical  methods  of  surface  preparation, 
before  being  painted  with  three  coat  bottom  painting  schemes  and 
immersed  in  the  sea  for  one  year. 

It  was  found  that  the  paint  did  net  adhere  well  to  freshly  de¬ 
scaled,  unrusted  steel.  Adhesion  was  greatly  improved  by  a  chemical 
surface  treatment  and  allowing  the  plate  to  rust  a  little  after  pick¬ 
ling  or  grit  blasting,  provided  that  the  plate  is  wire-brushed  immed¬ 


iately  before  painting. 


The  surface  preparation  of  steel  by  abrasive  blast  cleaning  methods 
is  a  well  established  procedure  with  both  technical  and  economic  ad¬ 
vantages.  The  technology  of  abrasive  blast  cleaning  had  been  investi¬ 
gated  by  many  workers  whilst  the  effect  of  the  resulting  surface  rough¬ 
ness  on  the  thickness  of  paint  coatings  applied  to  such  surfaces  had 
received  very  little  attention.  Daniel  (1969)  investigated  the  impor¬ 
tance  of  film  thickness  in  relation  to  the  anti-corrosive  properties  of 
paint  coatings.  On  a  rough  surface,  the  minimum  film  thickness  will  be 
present  over  the  largest  'peaks'  of  the  surface,  hence  this  minimum 
thickness  will  largely  determine  the  initial  protective  properties  of 
the  coating.  In  other  words,  an  abrasive  blasted  surface  will  require 
more  paint  to  keep  it  free  from  corrosion  than  a  smooth  surface.  The 
basic  problem  is,  how  much  more  paint  is  required  and  can  this  extra 
paint  be  estimated  by  a  consideration  of  the  properties  of  the  paint 
and  of  relevant  measurements  carried  out  on  the  substrate.  In  order 
to  help  solve  the  problems,  methods  for  predicting  and  measuring 
minimum  paint  thickness  v?ere  required.  The  investigations  were  restric¬ 
ted  to  a  consideration  of  surfaces  produced  by  the  abrasive  blasting  of 
fairly  smooth  steel  surfaces. 


The  characterisation  of  blast  cleaned  surfaces  in  terms  of  the 
volume  per  unit  area  bounded  by  the  metal  surface  and  a  plane  locatea 
by  the  summits  of  peaks  (i.e.  surface  volume)  was  examined  and  a  method 
of  measurement  was  developed.  The  effects  of  surface  roughness  upon 
paint  thickness  was  investigated,  the  most  significant  parameter  of  a 
blast  cleaned  surface  being  the  'surface  volume'.  Variations  in  paint 
thickness  shown  by  paint/substrate  systems  were  experimentally  measured 
on  cross-sections  of  the  test  surfaces.  A  method  of  estimating  the 
minimum  paint  thickness  over  peaks  was  proposed  for  systems  when  the 
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volume  of  dry  Paint  present  per  unit  area  is  Each  greater  than  the 
surface  volume.  Instruments  were  used  to  measure  paint  film  thickness 
and  it  was  shown  that  an  approximate  measure  of  the  volume  of  dry  paint 
present  per  unit  area,  and  a  very  approximate  indication  of  minimum 
thickness  could  be  achieved,  with  certain  limitations. 


1.4  General  surface  rouqhness/fluid  flow  considerations 


1.4.1  Practical  aspects  of  surface  frictional  resistance 

In  the  field  of  aerodynamics.  Young  (1950),  for  instance,  looked 
at  the  particular  problem  of  the  effect  of  camouflage  paint  on  profile 
drag  of  aeroplanes,  and  thus,  assessed  the  effects  on  their  performance. 
This  problem  was  of  some  importance  during  the  1939-45  war  and  the 
results  were  of  general  interest  for  workers  in  the  field  of  surface 
frictional  drag.  The  profile  drag  is  assumed  to  be  the  sum  of  the 
surface  frictional  drag  and  the  form  drag. 

Measurements  were  made  in  a  high  speed  wind  tunnel  of  the  drag  of 
a  wing  with  various  grades  of  surface  roughness  in  the  form  of  camou¬ 
flage  paint.  The  drag  effect  of  each  roughness  tested  was  such  that 
an  equivalent  size  of  sand- roughness  could  be  associated  with  it.  An 
analysis  of  roughness  records  taken  with  a  roughness  gauge,  developed 
by  Tomlinson  of  NPL,  showed  that  the  equivalent  sand  roughness  heights 
for  each  of  the  surfaces  were  about  1.6  times  the  arithmetic  mean  of 
the  roughness  heights  k,  and  in  general  had  similar  amplitudes  to  the 
largest  roughness  elements  occurring  with  moderate  frequency.  Estimates 
were  then  made  of  the  effects,  of  the  various  surface  finishes  tested, 
on  the  aircraft  drag  for  different  values  of  wing  chord  or  body  length. 


20 


forward  speeds  ana  mughts,  T 


He  results  were  tnen  used  to  assess  the 


effects  on  the  performance  of  representative  aircraft. 


More  recently,  Re  da  (197^1)  carried  out  experimental  work  in 
connectic  .  with  a  missile  development  program.  During  flight  test,  the 
missiles  were  suffering  from  a  reduction  in  range  which  could  not  be 
accounted  for  in  the  original  drag  calculations.  It  was  thought  that 
the  surface  of  the  missile  motor  casing  was  in  some  way  amplifying  the 
surface-frictional  drag  levels  beyond  the  estimated  values.  On  inspect¬ 
ing  the  motor  casing  surfaces,  'which  were  constructed  from  fibre  glass- 
wound  material,  they  were  found  to  have  a  large  percentage  of  their 
surface  covered  with  ridges  and  grooves  formed  during  fabrication . 

These  surface  undulations  were  in  a  transverse  direction  to  the  local 
flow  direction,  and  it  was  thought  that  the  resulting  random  rcugh/wavy 
surface  v?as,  in  some  way,  increasing  the  surface  frictional  resistance 
levels  beyond  the  estimated  levels  which  were  based  on  the  smallest 
roughness  dimension,  the  fibre  material  filament  radius. 


In  the  resulting  investigation,  surface  contour  traces,  using  a 
stylus  instrument,  were  taken  from  an  actual  motor  case  section,  and 
these  showed  the  existence  of  three  scales  of  roughness  i.e.  (1)  a 
roughness  scale,  (2)  a  short  wavelength  scale  and  (3)  a  long  wavelength 
scale. 


Several  rough/wavy  patterns  were  fabricated  in  an  attempt  to 
simulate  the  different  types  of  surface  features.  The  most  complex 
was  a  roughness  scale  superimposed  on  a  short  wavelength  waveform,  both 
of  which  Were  superimposed  on  a  long  wavelength  waveform. 
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Three  types  of  models  were  tested:  (1)  send  grain  models  for  com¬ 
parison  with  previous  results,  (2)  machined  models  formed  by  impressing 
the  desired  surface  patterns  c»  aluminium  plates  and  (3)  a  mould  of  the 
roughness  pattern  taken  directly  from  a  full-scale  motor  casing.  The 
results  of  the  investigation  showed  that: 

(a)  Superposition  of  a  given  roughness  dimension  on  a  surface  which 
already  possesses  long  and  short  periodic  waveforms  increases  the  sur¬ 
face  shear  above  the  levels  measured  for  identically  rough  surfaces 
possessing  only  one,  or  none,  of  these  same  periodic  waveforms. 

(b)  Equivalent  sand-grain  roughnesses  for  the  random  rough/wavy 
wall  (motor  case  material)  were  of  the  order  of  four  to  eight  times  the 
physical  roughness  scale  attributed  to  it,  the  physical  roughness  scale 
used  being  peak  to  valley  height. 

(c)  Velocity  profiles  measured  over  rough/wavy  surfaces,  including 
random  rough/wavy  and  rough/multiple  periodic  waveform  surfaces,  were 
found  to  possess  logarithmic  regions  for  which  the  law-of-the-wall  is 
valid. 


(d)  When  plotted  in  terras  of  velocity  defect  coordinates  (White 
1974),  the  smooth,  rough  and  rough/wavy  profile  data  collapsed  onto  a 
near  universal  curve  in  outer  portions  of  the  layer,  in  agreement  with 
previously  published  results,  verifying  that,  even  for  complex  surface 
patterns,  roughness  effects  are  localised  deep  within  the  boundary 
layer. 
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There  have  been  a  nunibsr  of  other  expariments  on  various  types  of 
roughness  carries  out  in  various  flow  configurations.  Open-channel 
flow,  for  instance,  has  been  used  in  several  investigations  of  the 
frictional  resistance  of  rough  surfaces.  Powell  (1946)  carried  out 
tests  in  a  rectangular  flume  on  various  types  of  roughness  made  up  of 
.i. >,nt  arrangements  ox  sguare  s  teel  strips  layed  at  right  angles  to 
the  flow  and  extending  down  the  sides  and  across  the  bottom  of  the  flume. 

3 

Measurements  were  also  carried  out  in  a  smooth  flume.  The  analysis  of 
the  results  v?ere  based  on  the  fundamental  work  on  open  channels  by 
Keulegan  in  1938.  The  roughness  effect  was  found  to  depend  both  on  the 
height  of  the  roughness  elements  and  their  spacing.  Einstein  and  Banks 
(1950)  carried  out  a  series  of  tests  to  determine  the  total  resistance 
of  different  types  of  roughness  opposing  the  flow  of  water  through  an 
open  channel.  Previous  work  had  mainly  been  confined  to  the  effect  of 
uniform  roughness  elements,  whereas  here  the  roughness  was  made  up  of 
various  densities,  patterns  and  combinations  of  different  shaped  rough¬ 
ness  elements.  The  aim  of  the  investigation  was  to  find  if  the  total 
roughness  frictional  resistance  could  be  expressed  as  the  sum  of  the 
individual  resistances.  On  a  much  larger  scale,  this  type  of  argument 
had  been  suggested  in  a  study  of  idle  Salinas  river  where  the  total  force 
opposing  the  flow  was  thought  to  consist  of  the  superimposed  resistance 
exerted  by  the  river  bed  and  sides,  the  sand  bars,  and  the  vegetation. 

An  analysis  of  the  experimental  results  of  the  investigation  showed 
that  it  was  possible  to  express  the  total  resistances  as  the  sum  of  the 
individual  resistances  as  long  as  the  separation  between  the  roughness 
elements  was  sufficiently  large  so  that  individual  roughnesses  did  not 


exert  any  mutual  interference  on  one  another. 


were  derived  for  uniformly  rough  surfaces  m  conduits  were  rouna  to  db 
valid  for  surfaces  with  variable  roughness  if  the  roughness  dimensions 
were  averaged.  Separate  friction  factor  equations  were  derived  for 
three  distinct  types  of  flow  regime  occurring  in  rough  pipes  and  channel 
Thus,  if  for  a  particular  rough  surface,  more  than  one  of  the  flow 
regimes  were  to  be  produced,  then  the  friction  factors  calculated  for 
each  regime  present  could  be  added  to  give  the  apparent  friction  factor 
for  the  surface  as  a  whole. 


Perry  and  Joubert  (1963)  showed  that  the  effect  of  roughness  in  the 
outer  flow  can  be  accounted  for  using  an  equivalent  viscosity  which  is 
assumed  to  be  dependent  only  on  the  variables  at  the  wall,  sucn  as  shear 
stress,  fluid  density,  viscosity  and  the  roughness  size  and  geometry. 
Using  a  wind  tunnel,  flow  measurements  were  taken  over  roughness  ele¬ 
ments  consisting  of  timber  strips.  The  work  surface  could  be  adjusted 
in  the  tunnel  to  produce  various  pressure  gradients.  The  resulting 
velocity  distributions  were  used,  along  with  a  modified  graphical  method 
to  determine  the  local  boundary  layer  parameters  for  each  of  the  flow 
regimes .  The  graphical  method  used  in  the  analysis  was  primarily  tor 
hydrodynamically  smooth  surfaces  in  flows  with  and  without  pressure 
gradients.  For  rough  surfaces  the  origin  of  z  is  not  known  and  there 
is  a  shift  in  the  logarithmic  profile,  hence  in  using  the  logarithmic 
line  to  accurately  calculate  C^,  only  the  slope  of  the  line  and  not  its 
position  might  be  used.  The  logarithmic  asymptote  is  assumed  to  be 
located  at  a  distance  e  below  the  roughness  crests.  Hence  to  use  the 
position  of  the  line  in  calculating  the  value  of  e  must  be  estimated 
an  error  in  the  value  of  the  z  origin  distorting  the  logarithmic  line 


into  a  curve. 


An  attempt  war  made  to  treat  the  surface  roughness  as  a  random 
process  {Kayak  19715  by  Singh  and  Lutiley  {1871}  who  examined  the  problem 
of  steady  incompressible  flow  past  a  semi-infinite  rough  wall  whose 
roughness  was  randomly  distributed  in  a  way  which  can  be  described. 
Assumptions  were  made  in  the  analysis  which  linked  the  average  size  and 
extent  of  the  roughness.  The  average  slope  of  the  roughness  was  assumed 


to  be  very  much  less  than  1  and  it  ‘was  therefore  expected  that  the 
additional  mean  velocity  and  pressure  due  to  the  roughness  would  be 


small  compared  with  those  corresponding  to  a  smooth  wall.  The  problem 
was  then  to  calculate  the  contribution  made  by  wall  roughness  in  modi¬ 
fying  the  mean  velocity  profile  in  the  laminar  boundary  layeri  Calcu¬ 
lations  were  made  to  find  the  roughness-induced  velocity  in  terms  of 
the  spectral  density  of  the  roughness  and  an  influence  factor.  It  was 
concluded  that  flows  over  rough  surfaces  invariably  become  less  stable, 
probably  because  in  a  wide  spectrum  of  roughness  there  will  always  be 
soma  roughness  wavelengths  present  which  will  make  the  shear  layer 
unstable  and  increase  the  surface  friction. 


Monzavi  (1872)  investigated  the  relation  between  natural  roughnesses 
and  their  frictional  resistance.  It  had  already  been  shown  that  charac¬ 
terising  the  surface  roughness  by  the  height  of  the  roughness  elements 
alone  is  insufficient  to  uniquely  define  a  surface,  therefore  in  this 
investigation. nine  roug  surfaces  were  tested,  flow  measurements  taken 
and  the  surfaces  statistically  evaluated.  Experimentally  determined 
frictional  resistance  values  were  correlated  with  possible  hydromechani¬ 
cal  significant  statistical  parameters  of  the  surface  geometry. 


Prom 
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Ives  a  measure  of  the  roughness  intensity. 
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(c)  The  absolute  mean  value  of  the  surface  slooe,  S  ,  . 
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In  t.v;o  papers  presented  by  Leukowicz  and  Musker  at  the  1978  Inter¬ 
national  Symposium  on  ship  viscous  resistance  at  Gotehorg  the  results 
from  study  of  the  effect  of  ship  hull  surfaces  on  fully  developed  pipe 
flows  were  presented-  A  replication  technique  was  used  which  allowed 
cylindrical  replicas  of  ship  hull  roughnesses  to  be  tested  in  a  pipe 
line  with  compressible  fully  developed  turbulent  air  flow.  The  replicas 
form  the  internal  surface  of  the  pipe  line.  Dynamic  similarity  is  used 
to  achieve  a  comparison  between  the  fully  developed  pipe  flow  and  the 
inner  regions  of  the  turbulent  boundary  layer  over  a  flat  plate  which 
is  assumed  to  be  equivalent  to  the  flow  over  a  ship  surface.  The  flow 
analyses  were  correlated  to  statistical  analyses  of  the  respective 
surface  roughness  (Thomas  1975) .  It  was  notec  that  marine  rough  surface: 
like  engineering  surfaces  are  unlikely  to  have  a  surface  texture  similar 
to  Nikuradse’s  tightly  packed  ‘equivalent’  sand  grain  roughnesses 
(Nikuradse  1933) .  Engineering  surfaces  are  generally  random  and  have 
average  slopes  of  a  few  degrees.  Marine  surfaces  are  more  likely  to 
have  randomly  positioned  areas  of  different  roughnesses  with  abrupt 
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Reynolds  number  and  frictional  resistance. 

Manning  (1979)  investigated  the  influence  of  surface  roughness  on 
the  gas  flow  through  cracks.  Gas  flow  rates  were  measured  through  a 
number  of  idealised  cracks  whose  surfaces  'were  roughened  by  grit- 
blasting.  Crack  surfaces  were  used  with  roughnesses  between  0,7  and 
13  pm  R_,  and  with  a  number  of  different  widths.  The  results  were 
analysed  in  terms  of  friction-factor /Reynolds  number  plots.  The  power 
spectral  density  function  was  used  in  the  analysis  of  the  surface  top¬ 
ography.  The  computer  power  spectrums  showed  the  surfaces  could  be 
approximated  to  bandwidth-limited  white  noise. 


1.4.4  Flow  over  curved  surfaces 

The  work  so  far  reviewed  has  been,  in  the  main,  limited  to  analy¬ 
sing  the  friction  effect  of  the  relatively  short  surface  wavelengths. 
Since  a  considerable  amount  of  work  correlates  the  flow  measurements 
over  flat  surfaces  with  the  flow  over  the  curved  hull  of  a  ship  it  is 
of  considerable  interest  to  investigate  the  effect  of  curved  surfaces 
on  fluid  flow. 


Hsroney  and  Bradshaw  (1975)  investigated  turbulent  boundary  layers 
over  convex  and  concave  surfaces.  The  experimental  set-up  consisted 
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Characteris ing  surfaces  comprisin' 


£ 


In  order  to  determine  the  roughness  allowance  which  must  be  applied 
to  data  from  model  tests  irt  order  to  Sake  it  possible  to  accurately 
predict  the  resistance  of  the  ship,  it  is  essential  to  have  a  means  of 


characterising  a  rough  surface. 


Prom  the  classic  work  of  Nikufadse  (1933) ,  whose  rough  surfaces 
were  made  of  closely  packed  sand  grains,  a  great  deal  of  subsequent 
work  was  concerned  with  establishing  an  equivalent  sand  rougnness  for 
each  surface  considered,  so  that  the  resistance  could  be  obtained  from 
Nikuradse’s  results.  Simple  averages  and  root-mean- square  (rhS)  values 
of  the  roughness  amplitude  were  used  and  interpreted  in  terms  of  equiv- 


28 


ct ,1  Gii t  Sanci  roughness .  As  the  worn  progressed  it  became  clear  that  many 

forms  of  roughness  could  hot  be  defined  satisfactorily1  by  an  equivalent 
sand  roughness. 


1-5.2  Characterising  surfaces  of  non-uniform  dimensions 


In  most  practical  cases,  the  roughness  elements  are  of  non-uniform 
site,  so  that  the  'completely  rough*  regime  obtained  by  Nikuradse  in 
which  the  surface  friction  is  independent  of  Reynolds  number  and  depend¬ 
ent  only  on  the  relative  roughness,  is  hot  valid.  For  example,  Todd 
(1951)  lias  shown  that  ship  hulls  are  hardly  ever  'completely  rough'  in 
this  sense. 


In  attempting  to  characterise  practical  rough  surfaces  it  became 
clear  that  a  statistical  description  of  the  surface  by  centre  line 
average  (CL&)  or  RMS  roughness  amplitudes  and  other  similar  parameters 
was  not  adequate.  For  example,  two  surfaces  could  have  the  same  CLZl  ' 
or  RMS  roughness  amplitudes  but  have  completely  different  textures  and 
probably  produce  different  amounts  of  frictional  resistance  to  fluid 
flow  over  them  (Eisenberg  1950) .  This  need  to  have  a  more  comprehensive 
method  of  surface  characterisation  was  also  true  in  the  case  of  engin¬ 
eering  surfaces  (Peklenik  1968)  where  the  interest  is  in  the  study  of 
a  number  of  interfacial  phenomena  such  as  friction  and  wear,  and 
electrical  and  thermal  contact  resistance. 


Returning  io  the  problem  of  the  frictional  resistance  of  rough 
surface  (Eisenberg  1950) ,  on  looking  at  simple  cases  of  surfaces  with 
the  same  average  amplitude  parameters. (Figure  Al-13)  it  is  clear  that 


an  additional  parameter  to  cor.ro  Late  the  surface  With  its  frictional 
resistance  is  the  slope  of  slope  distribution . 

Thus,  an  adequate  description  for  the  uniform  two-dimensional  rough¬ 
ness  could  be  built  up  using  the  averages  or  frequency  distributions  of 
the  amplitudes,  slopes  and  curvatures  (Posey  1946) .  In  addition  to  the 
above  geometrical  properties,  for  rough  surfaces  in  practice,  it  was 
thought  necessary  to  specify  the  roughness  distribution  in  the  direction 
of  motion  and  the  traverse  direction,  and  the  roughness  density  in  the 
direction  of  motion.  It  was  also  thought  prudent ,  in  these  early  stages 
of  characterising  surfaces,  to  look  at  the  spectral  distributions  and 
correlations  (as  between  the  longitudinal  and  traverse  distribution):  of 
the  above  parameters. 


To  help  in  the  characterisation  of  practical  surfaces,  a  surface 
measuring  instrument  was  developed  in  order  to  make  it  possible  to 
produce  chart  recordings  of  surface  profiles  from  the  movement  of  a 
probe  which  follows  the  surface .  The  movement  of  the  probe  is  fed 
through  a  linkage  system  to  strain  gauges.  The  resulting  electrical 
signal  is  suitably  processed  and  used  to  drive  a  pen  in  the  recorder. 

The  instrument  was  designed  in  order  to  be  usable  in  situ  on  ships  hulls 
(T6dd  1951,  Hama  1954). 

A  surface  roughness  measuring  device  was  developed  by  the  British 
Ship  Research  Association  (BSRA)  (Canhsm  1956)  in  order  to  enable  them 
to  take  comprehensive  roughness  records  of  ship  hulls;  These  records 


could  then  he  used  ai0ng3i.de  the  ship  trial  data.  The  roughness  igeasijre 
ments  are  made  in  dry  dock  shortly  before,  or  immediately  after  trials 
and  involve  the  recording  of  roughness  profiles  at  a  number  of  positions 
on  the  hull  (Deer  1977) . 


A  statistical  analysis  was  used  to  examine  the  various  factors 
affecting  the  trial  performance  and  it  was  found,  as  expected,  that  the 
hull  plate  roughness  Was  by  far  the  most  significant  factor.  The  nature 
of  the  roughness  varies  in  different  parts  of  a  hull  (for  example,-  paint 
drops  tend  to  accumulate  on  the  flat  of  the  bottom)  and  alternatively., 
the  same  degree  of  roughness  may  well  have  different  effects  on  resist 
tance  according  to  the  position  on  the  ship's  hull. 

A  technical  memorandum  by  Canhaa  (1961)  set  out  to  give  some 
guidance  on  the  selection  or  recording  positions  and  the  number  of 
records  required  to  be  taken  with  the  BSRA  wall  roughness  gauge  in 
order  to  make  a  satisfactory  null  roughness  survey.  A  hull  roughness 
survey  is  required  to  establish  whether  the  surface  of  a  hull  has  a 
roughness  within  predetermined  limits,  or  to  determine  the  average 
roughness  of  the  hull  surface,  or  to  make  an  accurate  assessment  of  the 
frictional  resistance  of  the  hull.  The  hull  surface  is  divided  into 
definite  areas  in  each  of  which  there  is  a  significantly  different  type 
of  surface  roughness .  It  is  therefore  essential  to  survey  as  much  as 
possible  of  the  hull  before  deciding  on  the  number  of  records  and  the 
location  of  the  recording  positions  required.  The  number  of  records  to 
be  taken  should  increase  With  increasing  average  roughness  and  with 
increasing  non-uniformity  of  the  roughness,  and  should  not  be  less  than 


fifty  for  a  'reasonably  concise’  hull  survey 


Roughness  records  were  taken  on  a-  sample  of  68  new  Ships  (Lackenby 
1962}  and  the  resulting  distribution  of  average  toughness  was  found  to 
approximate  to  a  normal  gaussian  distribution  with  a  moan  roughness  valu 
of  0.0074  inch  (16.8  pm)  (Figure  Al-14).  Typical  roughness  tracers  for 
these  surfaces  are  shown  in  Figure  Al-12. 


In  a  subsequent  BSRA  report  dealing  with  the  various  methods  of 
analysing  hull  surface  roughness  records,  Chaplin  (1965)  proposed  that 
spectral  analysis  could  be  successfully  used  as  a  method  of  obtaining 
alternative  roughness  parameters  to  those  currently  in  use.  It  was  also 
shown  that  commonly  used  parameters  such  as  the  mean  apparent  amplitude 
and  the  centre  line  average  could  be  deduced  from  the  spectral  analysis. 
The  previously  accepted  method  of  analysis  of  hull  surface  roughness 
records  had  been  to  determine  the  mean  apparent  amplitude  for  a  sample 
length  of  2  inches.  Statistical  analysis  of  numerous  ship  trial  aha 
corresponding  model  test  results  indicated  the  existence  of  a  signifir 
cant  relationship  between  frictional  Resistance  parameters  based  on  the 
mean  apparent  amplitude  per  2  inch  length  of  hull  surface  and  the 
corresponding  derived  ship-model  correlation  factor.  It  was  clear  that 
alternative  methods  of  analysing  roughness  records  had  to  be  considered 
in  the  hope  of  throwing  more  light  on  the  relationship  between  surface 
roughness  and  ship  resistance  and  in  improving  the  ship-model  corre¬ 
lation  factors.  The  records  obtained  showed  the  profile  to  be  nearly 
always  random  in  nature  (Nayak  1971)  making  it  impossible  to  represent 
the  profile  by  any  simple  mathematical  function.  Much  the  same  sort  of 


problem  existed  in  the  study  o.f  the  sea  surface  and  electrical  'noise', 
and  this  resulted  in  the  development  of  the  technique  of  spectral 
analysis  to  a  high  degree.  Thus,  Chaplin  went  on  to  represent  the 


profile  by  the  sum  of  a  largo  number  of  harmonic  components  in  which 
the  amplitudes  and  phase  angles  were  random  variables.  These  quantities 
could  be  estimated  by  suitable  numerical,  techniques  and  an  amplitude 
spectrum  constructed. 

Just  as  the  energy  spectrum  can  be  used  to  resolve  white  light  or 
a  complex  sound  into  a  sum  of  pure  colours  or  tones,  so  a  roughness 
profile  can  be  resolved  into  a  series  of  component  sine  waves  each  with 
a  definite  frequency  and  contributing  a  definite  proportion  of  the  total 
variance  in  the  profile. 

Apart  from  providing  a  distribution  of  the  amplitudes  and  fre¬ 
quencies  of  the  component  waves  in  the  profile,  several  useful  measures 
of  roughness,  including  the  mean  apparent  amplitude  and  centre  line 
average,  and  other  statistical  properties  can  be  derived  from  the 
spectrum. 
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FIGA1-1  SHIP  RESISTANCE -RESULTS  FROM 
WIND-TUNNEL  TESTS  AND  TANK  TESTS  ON  A 
c  20  FT.  LONG  MODEL  (Lackenby  1962) 
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FIGA1-2  FROUDES  RESULTS,  FROM  TOWING  WOODEN 
BOARDS  IN  A  TANK, ANALYSED  ON  A  REYNOLDS 
NUMBER  BASE  (Lackenby  1962) 
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FIG  A1-B  COMPILED  EXPERIMENTAL  DATA 
ON  THE  FRICTIONAL  RESISTANCE  OF 
PLANE  SURFACES  ( Schoenherr  1932 ) 


FIG  A1-3  COMPILED  EXPERIMENTAL  DATA 
ON  THE  FRICTIONAL  RESISTANCE  OF 
PLANE  SURFACES (Schoenherr1932) 
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FIG  A1-4  PLOT  °F  NEW  AND  OLD  FRICTIONAL 
RESISTANCE  DATA, PLOTTED  TO  ABASE  OF 
LENGTH /BREADTH  RATIO{Hughes  1952) 
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FICA1-7  CORRELATION  OF  SHIP  AND  MODEL 
RESULTS  FOR  LUCY  ASHTON ( Lackenby  1962) 


FIG  A1-8  COMPARISON  OF  FULL-SCALE  RESISTANCE 
WITH  SCHOENHERR  SMOOTH  SHIP  PREDICTIONS 
FROM  MODELS  (Conn  etal.  1953) 
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FIG  A1-9  LOSS  IN  PERFORMANCE  WITH  TIME  IN 
SERVICE  FOR  THE  CROSS-CHANNEL  SHIP 
KONINGEN  ELISABETH  ( Lackenby  1962) 


FIGA1-10  PERCENTAGE  INCREASE  IN  TRIAL  SHP 
OVER  THE  ORIGINAL  TRIALtLsckenby  1962) 
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FIG  A1-11  PERFORMANCE  OF  18  000 -TON  DEADWEIGHT 
TANKERS.  VARIATION  IN  TRIAL  SHP  WITH  HULL 
ROUGHNESS(Lackenby1962) 
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FIG  A1-12  SURFACES  WITH  GEOMETRIES  WHICH 
CANNOT  BE  DEFINED  BY  AMPLITUDE  PARAMETERS  ALONE 
(Eisenberq  1950) 
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In  analysing  the  effect  of  surface  roughness  on  the  surface 
frictional  resistance  coefficient  C^,  a  number  of  fundamental  relations 
are  available  to  predict  the  motion  of  a  fluid  flowing  over  a  solid 
surface.  In  fact  the  fundamental  equations  of  viscous  flow  have  been 
known  for  more  than  a  hundred  years.  At  low  Reynolds  numbers  (laminar 
flow)  the  equations  are  very  difficult  to  solve,  even  on  digital 
computers,  and  at  high  Reynolds  numbers  (turbulent  flow)  the  equations 
are,  at  present,  unsolvable  because  of  the  random  nature  of  flow  fluc¬ 
tuations,  in  particular  the  boundary  conditions. 

Although  the  exact  solutions  to  the  fundamental  equations  of  com¬ 
pressible  viscous  flow  cannot  be  obtained,  a  great  deal  of  useful 
information  can  be  released  from  simplified  solvable  forms  of  these 
exact  equations.  A  pioneer  in  this  field  was  Prandtl,  who  using  the 
boundary  layer  concept,  which  stated  that  the  viscous  effects  are  con¬ 
fined  to  a  very  thin  layer  near  the  solid  surface,  gave  approximate 
solutions  of  the  three  laws  of  conservation  for  the  case  of  two-dimen¬ 
sional  flow  over  a  plane  surface. 

2.2  The  three  laws  of  conservation 

The  number  of  viscous  flow  equations  used  is  generally  reduced  to 
the  three  laws  of  conservation: 

(1)  Conservation  of  mass  (continuity). 


(2)  Conservation  of  momentum  (Newton's  second  law) . 

(3)  Conservation  of  energy  (first  lav/  of  thermodynamics) . 

The  exact  number  of  basic  flow  equations  depends  on  personal 
preference,  some  relations  being  more  basic  than  others.  Truitt,  for 
example,  used  nine  basic  equations  relating  to  nine  unknowns  (White  3S74 
In  the  present  investigation  the  flow  is  assumed  to  be  incompressible 
and  isothermal,  hence  the  three  laws  of  conservation,  for  fluid  flow, 
reduce  to: 

(1)  The  equation  of  continuity,  which  in  its  commonest  form  is 
given  by: 

“■  +  pdiv  V  =  0  or  dp  (2-1) 

Dt  —  +  div  pV  =  0 

dt 

where  v  is  the  three-dimensional  velocity  vector,  v  =  iu  +  jv  +  kw 

„  du  ,  dv  ,  dv/ 
and  div  V  =  —  +  ——  +  — 
dx  dy  az 

If  the  density,  p ,  is  constant  (incompressible  flow)  the  equation  of 
continuity  reduces  to 

div  V  =  0  (2-2) 

(2)  The  conservation  of  Momentum  equation,  com.,  lly  known  as 
Newton's  second  law,  is  given  by 


DV 

P^r 


-  pg  +  V.t. 


(2-3) 


where  o  -  gravitational  acceleration  and  V . t  -  dt  /< 3x. 

‘ij  -3'  J 

V  is  -the  gradient  operator,  -  is  the  stress  tensor,  which  for  a 


Newtonian  viscous 
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sneral  deformation  law 
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+  6  .  .C,  div  V 
1  j  bv 


(2-4) 


where  1  is  called  the  coefficient  of  bulk  viscosity.  Thus,  the  momentum 
equation  for  viscous  Newtonian  flow  is  obtained  by  substituting  the 
general  deformation  lav/  (equation  2-4)  into  Newton's  second  lav?  (equation 
2-3) .  The  resulting  equation  of  motion  originally  given  by  Navier  and 
Stokes  can  be  written  as 


DV  „  ±  a 

Dt  ^  dx. 


dv. 

1  + 

dv. 

_ i_ 

+  5  .  .Ct  div  V 
xi  bv 

■dx. 

dx.  3 

(2-5) 


For  incompressible  flow,  the  density  is  constant  and  div  V  vanishes 
from  the  Continuity  equation  (2-2)  and  the  coefficient  Cfav  vanishes 
from  Newton's  lav/.  Also  assuming  that  viscosity  is  constant,  the  Navier- 
Stokes  equations  can  be  considerably  simplified  to 


pDt  =  pg  "  Vp  +  |iVJv 


(2-6) 


Summarising,  the  two  laws  of  conservation  of  mass  and  momentum,  as 
used  for  fluid  motion  are 


dp 

dt 


+  div  pV  =  0 


(2-1) 


ny 

p—  =  pg  +  V.t.  .  -  Vp 

MDt  y  13 


(2-3) 


where  for  Newtonian  fluid,  the  viscous  stresses  are 
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ij  bv 


(2-4) 


Equations  2-1  and  2-3  involve  four  variables,  of  which  o  and  v 
are  assumed  primary,  and  the  remaining  two  variables  are  assumed  known 
from  auxiliary  relations  and  data  of  the  form: 


p  =  p(p,  T>  ;  |i  =  |X(p#  T) 


2.3  Approximate  solutions  of  the  laws  of  conservation  using  the 
boundary  layer  hypothesis  _ _ _  _ _ 


As  can  be  seen,  the  equations  to  be  solved  are  a  formidable  set  of 
non-linear  partial  differentials  equations.  No  method  of  solving  these 
equations  for  an  arbitary  viscous  flow  problem  exists,  although  digital 
computer  numerical  techniques  using  finite-difference  approximations 
show  considerable  promise.  Thus,  approximate  solutions  are  obtained 
using  generalised  forms  of  the  exact  equations. 

Several  exact  solutions  have  been  obtained  for  specific  types  of 
Newtonian  viscous  flow  (White  1974) .  The  high  Reynolds  number  flows 
considered  are  explainable  using  a  boundary  layer  hypothesis.  The 
boundary  layer  theory  states  that  the  viscous  effects  are  confined  to 
relatively  thin  layers  adjacent  to  solid  boundaries.  In  the  remaining 
region  av?ay  from  the  wall  the  influence  of  viscosity  is  assumed  to  be 
unimportant  and  the  flow  is  frictionless  and  potential.  Hence  Prandtl 
(1904)  simplified  the  Navier-Stokes  equations  to  provide  approximate 
solutions  for  two  dimensional  incompressible  boundary  layer  flow  on  a 
plate. 
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in  dimensional  form,  for  steady 


n'r.nq 

isothermal 


cha  boundary  layer  equations, 
two-dimensional  incompressible 


flow  are 


dn  dv 

dx  da 


0 


(2—7) 


du  udu  vdu  f  du  dU  yd2  u 

‘dt  dx  T  dz  dt  dx  '  dz2 


(2-8) 


These  are  numerically  solvable,  for  practically  any  distribution  of 
free  stream  velocity.  The  solution  of  the  boundary  layer  equations 
involves  a  great  deal  of  labour,  and  so  the  introduction  of  digital 
computer  techniques  lightened  the  burden  in  this  respect.  A  number  of 
finite-difference  solutions  are  available  using  both  explicit  ana  im¬ 
plicit  methods.  The  methods  involve  examining  the  flow  at  the  different 
stations  in  the  x-direction.  The  explicit  methods  (Gortlsr  1939)  have 
the  advantage  of  simplicity  but  give  inconsistent  results  unless  small 
steps  Ax  are  used,  usually  meaning  excessive  computation  time.  The 
implicit  methods  (Smith  6  Clutter  1963)  are  relatively  complicated  but 
give  more  stable  results. 


2.4  Integral  forms  of  the  basic  boundary  layer  equations  for 
laminar  isothermal  two-dimensional  incompressible  flow 


Turning  our  attention  to  more  approximate,  rapid  and  simpler 
techniques  we  consider  the  integral  equations  of  the  boundary  layer. 
The  basic  boundary  layer  equations: 


Continuity  - 


du  dv 

dx  dz 


0 


(2-9). 


Momentum  - 


du  udu  vdu 

dt  dx  dz 


r  dU  .  UdlJ  ’ 
dt  dx 
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1  dt 
p  dz 


(2-10) 
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for  the  case  of  steady  flow  ever  an  Impermeable  wall. 


The  H  or  shape  factor  is  always  greater  than  unity  (Figure  212-1) , 
and  for  laminar  flows  varies  from  a  value  of  about  2  at  the  stagnation 
point  to  about  3.5  at  separation.  In  turbulent  flow,  the  variation  is 
between  about  1.3  to  2.5. 


Continuing  with  the  derivation  of  the  integral  relations,  the 
mechanical-energy  integral  formed  by  multiplying  continuity  by  u5  -  U2 
and  momentum  by  2u,  subtract,  and  integrate  as  before  to  get 
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P  J  dz 
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(2-15) 


+  ib  f  u 


(U2  -  u*  )dz  -  U2  v 


which  can  be  reduced  to  J 


C  =  — 

D  pU 


_  ,  20_  du  i_  £l_ 

i*  U  dt  0  t  V  U2  dt  U1  dx  (~  3}  U  {2  16) 


where  D  is  the  shear  dissipation  term 


w 

m,  D  =  Tdu 
0J  dz 


63  is  the  kinetic  energy  thickness  or  dissipation  thickness, 


u  ,  J  u‘ .  , 

U  **  “  u2'aZ 


(2-17) 


(2-18)  and  is  the  energy  dissipation 


coefficient. 


-11  and  2-16  are  ths 


isothermal  laminar  two-dimensional  incosoressible  flow. 


2.5  Turbulent  boundary  layer  equsticnE 


Laminar  frow  is  inherently  unstable  and  above  a  critical  .  .olds 
number  there  is  an  abrupt  change  to  turbulent  flow .  Turbulence  is  an 
experimentally  observed  phenomenon  and  there  is  no  existing  theory  which 
states  that  turbulent  flow  should  occur  at  high  Reynolds  numbers.  The 
transition  process  can  be  made  to  occur  earlier  by  the  introduction  Of 
an  obstruction  or  obstructions  at  the  wall  such  as  roughness  elements, 
causing  additional  disturbances  in  the  boundary  layer,  which  distort 
the  velocity  profile  and  lower  the  critical  Reynolds  number.  The  trans¬ 
ition  can  also  be  affected  by  pressure  gradients,  wall  cooling,  or  free 
stream  turbulence.  Thus,  for  the  majority  of  cases  of  low-viscosity 
fluids  ('water,  mercury,  gases  ....)  the  flow  is  normally  turbulent,  hot 
laminar. 


Using  the  assumption  that  fluid  is  in  a  randomly  unsteady  turbulent 
state  the  time-averaged  or  mean  equations  of  motion  can  be  used  to 
approximate  the  How.  Thus 


du  dv  da  „ 

—  +  —  +  hr.  =  q 

dx  .  dy  dz 


(2-i5) 


(u'  .  u’  .) 
i  3 


=  pg  -  ?p  +  ji7J  V 


{2-20} 


where  the  turbulent  inertia  term  u'^  u'..  is  dependent  on  not  only  the 
physical  fluid  properties  but  also  on  local  flow  conditions  (velocity. 


geometry,  surface  roughness,  and  upstream  history) . 


which  there  is  no  avaiiabl 
that  the  turbulent  inertia 


seans  of  defining  their  values.  On  assumitig 
.eras  can  be  looked  upon  as  stress  terms  the 


mean  momentum 
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(White  1974)  where  t 
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The  turbulent  shear  -pu'v'  is  the  dominant  term  in  any  attempt  to 
resolve  the  momentum  equation  in  the  boundary  layer.  In  attempting  to 
resolve  the  turbulent  flow  momentum  equation  it  is  far  simpler  initially 
to  look  at  the  case  of  two-dimensional  boundary  layers.  One  method  of 
resolving  the  term  -pu'v'  is  to  assume  the  terra  to  be  a  function  of 
local  geometry  and  flow  conditions,  whilst  another  is  to  look  at  the 
local  velocity  profile  itself. 


Using  the  first  of  these  approaches,  PrandtX  (1926)  related  the 
turbulent  shear  to  the  local  mean  fluctuations  using  an  eddy- transport 
theorem  based  on  the  much  smaller  scale  kinetic  theory1  molecular  trans¬ 
port  model.  The  local  RMS  fluctuations  were  related  to  the  local 
velocity  gradient  through  a  turbulent  length  scale  £, 


^  (2-22) 
1  dz 

/V*2  =  fi  (2-23) 

i  dz 


where  £.  and  £«,  called  the  mixing  lengths,  are  analogous  to  the  mol- 


ocular  mean  free  path.  Prandtl's  mixing-length 


;1  can  then  be 


written  as 


H  H 


where  |£,  =  p4b  Hr! 


(2-24) 


Other  approaches  included  the  introduction  of  an  eddy  viscositv  e 
by  setting  -u' v~  =  e,  du/dz  or  by  utilising  the  process  by  which  the 
boundary  layer  acquires  additional  turbulent  fluid,  entrainment  {Prcc, 
coraput.  of  Turb.  Boundary  Lavers  1358) . 


2.6  The  general  form  6f  the  boundary  layer  velocity  distribution 


The  velocity  profile  across  the  boundary  layer  can  be  split  into 
three  regions:  The  region  next  to  the  wall  where  viscous  shear  dominates, 
and  an  overlap  region  where  both  types  of  shear  are  present.  Relation¬ 
ships  have  been  put  forward  to  predict  the  shape  of  the  profile  in  these 
regions. 


Prandtl  in  1333  deduced  that  in  the  inner  region  the  mean  velocity 
depends  upon  the  wall  sheaf  stress,  the  fluid  properties,  and  the  dis¬ 


tance  z  from  the  wall;  Thus  the  law  of  the  'wall  is  of  the  form 


u  .f  zu  I 

-  -  f  — 

T  v  V  3 


(2-25) 


For  the  outer  region,  Karman  in  1330  deduced  that  the  wall  acts 
only  to  retard  the  flew,  reducing  the  local  velocity  u  below  the  free 
stream  value  Uss  independently  of  viscosity  but  dependent  upon  the  wall 
shear  stress  and  the  distance  z  over  which  the  flow  is  affected. 
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(U  -  u)/u  =  0(2/0) 


(2-26) 


For  the  overlap  region,  the  two  functions  are  merged  together  over 
seme  finite  region  between  the  inner  and  outer  (Millikan  1936) .  Overlap 


[«r 

it  v  J 

UJj 

U 

—  CO  — 

U 


(2-27) 


For  the  inner  part  of  the  overlap  region  the  relation  is  given  b 


u  1  /•zu  ^ 

—  -  —  log  t 

u  r  3e  - 

T  {  v 


+  B 


and  for  the  outer  part 


U  -  u 

CO 


1 


W  „ 

—  +  A 

1°; 


where  sc  is  Von  Kansan's  constant  and  T\  and  B  are  constants. 


For  the  inner  region,  the  form  of  the  law  of  the  wall  can  be 
solved  by  assuming  that  the  total  shear  is  constant  very  near  the  wall 
and  that  the  eddy  viscosity  follows  a  cubic  relation  near  the  trail  and 
an  exponential  relation  further  out  (Spalding  1361)  to  produce  a  for¬ 
mula  for  the  inner  law  which  was  found  to  be  equally  valid  in  the 
logarithmic  overlap  region  (Figute  A2-?.) . 


55 


Although  all  the  profiles  for  the  inner  law  collapse  onto  ?  ijii- 
versal  curvs  (Fwiirs  P* 2—3)  xt  can  bs  sesi^  that  at  too  end  of  the 
curve  i.e.  in  the  overlap  region,  there  are  increasing  deviations  frca 
the  predicted  curve.  Thus  the  outer  law  is  such  more  sensitive  to 
external  influences,  such  as  pressure  gradients.  The  addition  of  a 
par arete r  involving  pressure  gradient  to  the  outer  law  was  first  attemp¬ 
ted  by  Clauser  (1S54,  1956). 

The  work  previously  done  on  the  problem  of  the  turbulent  boundary 
layer  had,  in  the  Bain,  been  directed  towards  looking  at  the  simplest 
of  cases,  such  as  flow  over  e  smooth  surface  with  zero  pressure  gradient. 
Only  a  few  empirical  methods  existed  for  predicting  the  effect  of 
pressure  gradients  on  turbulent  boundary  layers  and  these  were  unreli¬ 
able,  doing  little  more  than  to  correlate  the  set  of  data  on  which  they 
were  based  originally. 

For  some  time  it  had  been  realised  that  the  state  of  the  turbulent 
boundary  layer  depended  not  only  upon  the  local  environment  such  as  the 
local  values  of  pressure  gradient,  the  wall  shear,  the  surface  roughness, 
the  boundary  layer  thickness,  etc,  but  also  upon  the  previous  history 
of  the  layer.  Thus,  when  Clauser  (1954)  carried  out  his  research  on  the 
turbulent  boundary  layer  he  began  by  attempting  to  study  in  detail  a 
few  cases  that  had  single  and  well  defined  pressure  histories.  The 
work  was  carried  out  on  boundary  layers  having  relatively  constant 
histories,  which  provided  a  particularly  simple  and  well  defined  case 
to  study.  The  resulting  set  of  profiles  for  these  so  defined  pressure 
histories  were  celled  “equilibrium"  profiles. 
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Thus,  trial  piessuj 
profiles  at  a  number  of 
coefficients  were  de  ten, 


e  distributions  would  be  set  up  and  the  velocity 
stations  would  be  measured.  The  surface- friction 
lined  and  (u  -  U,5)/u  plotted  against  z/o  for 


each  station.  If  the  velocity  profiles  were  not  similar,  the  pressure 
distribution  was  altered  and  the  procedure  repeated  until  similar  pro¬ 


files  were  obtained. 


Clauser  found  that  the  plots  of  (u  -  fJj/'u  aaainst  z/6  produced 

e  T 

a  universal  curve  which  showed  an  expected  difference  from  the  universal 
curve  at  constant  pressure  (Figure  A2-4) .  Next  it  was  decided  to  ex¬ 
periment  with  a  steeper  dimensionless  pressure  distribution.  When 
analysing  the  results  it  was  found  that  many  of  the  previously  used 
boundary  la} nr  parameters  were  unsatisfactory.  The  boundary  layer 
thickness  has  the  difficulty  that  the  outer  edge  of  the  turbulent  layer 
is  ill-defined.  The  displacement  thickness 


CO 


0 


and  the  momentum  thickness 


(2-29) 


could  not  be  obtained  directly  from  the  universal  plot  that  formed  the 
basis  of  the  correlation  of  turbulent  profiles. 


Thus,  parameters  were  required  whicn  did  not  need  additional  data 
for  their  specification.  The  simplest  universal  thickness  so  defined 


I; 


|! 

!  • 
£ 


was 


CO  03 


0  1  0 


(2-30)  - 


For  constant  pressure  profiles  the  integral  has  the  value  3.6,  and  for 

V 

each  set  of  equilibrium  profiles  there  was  a  different  value  for  the 
integral. 


In  addition  to  the  defect  thickness  A,  a  suitable  parameter  was 
required  to  c) aracterise  the  various  members  of  the  universal  family  of 
profiles.  A  ratio  of  moments  on  the  universal  plot  was  used, 


CO  00  CO 


The  new  thickness  A,  and  the  integral  family  parameter  G,  were  related 
to  6  and  6  as  follows 


(2-32) 


(2-33) 


.  The  assumption  that  the  entire  profile  xs  represented  by  the 

universal  function  of  (u  -  U) /u  versus  z/6  was  known  to  be  not  true 

T 

near  the  wall.  The  calculated  values  of  the  integrals  for  A  and  G  for 

a  given  distribution  of  (u  -  U)  /u  versus  z/6  in  the  outer  portion  of 

T 

the  layer  were  applicable  only  as  long  as  the  surface  roughness  elements 
are  small  compared  with  the  thickness  of  the  viscous  sublayer.  Pressure 
gradients  were  found  to  have  an  unexpectedly  large  effect  on  surface 
friction.  By  combining  his  work  with  the  results  of  previous  researchers 
(Ludwieg  &  Tillman  1950, .Schubaner  &  Klebanoff  1950)  a  more  complete 
picture  of  turbulent  surface  friction  was  produced. 


58 


The  univar: 


vej.ocj.t 


..ion  ior 


cbuient  velocity  pro 1 


Les 
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layer  in  which  the  v 


:ure  2-3)  showed  a  region  outside  the  viscous  sub- 
ocitv  distribution  could  be  represented  by 


u 

u 

T 


5.6  log1Q 


+  4.9 


(2-34) 


even  when  pressure  gradients  are  present. 


For  the  case  of  a  rough  wall,  it  had  already  been  shown  by  Prandtl 
that  the  semilogarithmie  part  of  the  curve  was  displaced  downward  para¬ 
llel  to  itself  by  an  amount  Au./u  ,  which  depends  on  the  roughness 

i  T 

Reynolds  number  ku  /v  •  Plots  of  experimental  data  on  the  relations 

between  these  two  parameters  for  various  types  of  roughness  showed  that 

when  ku  /v  is  greater  than  approximately  100,  the  asymptotic  form 
t  * 

is 

Au/u  =  5.6  log(ku  /v  )  +  constant  (2-35) 

x 

Assuming  that  zu^/v  at  the  outer  edge  of  the  viscous  sublayer 
is  approximately  50,  then,  when  the  roughness  is  greater  than  about 
twice  the  thickness  of  the  viscous  sublayer,  the  universal  lav;  for  the 
turbulent  region  near  the  wall  becomes 

u/u  =  5.6  ,log(z/k)  +  constant  (2-36) 

T 

The  above  result  held  for  smooth,  moderately  rough,  and  extremely 
rough  walls  without  pressure  gradient  and  for  smooth  walls  with  pressure 
gradients.  Using  dimensional  reasoning  showed  that  the  effects  of 
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pressure  gradients  on  conditions  near  the  wall  should  be  small,  it  then 
followed  that  the  results  we  re  true  for  rough  walls  with  pressure 
gradients.  For  all  these  cases,  there  existed  a  turbulent  region  next 
to  the  wall  for  which  u/u  is  linearly  related  to  log  z. 

Looking  at  the  ’universal  clot  of  (u  -  U)  /u  versus  z/A  it  was  seen 

"  x 

that  for  the  inner  portion  of  the  plot,  (u  -  U)  /'u  should  be  linearly 

proportional  to  log  z/A.  Further  the  slope  of  this  semi logarithmic 

curve  should  be  5.6,  the  same  as  that  on  the  plot  of  u/u  versus  loa 

T 

zu^/v  or  log  z/k  (Figure  A2-5) .  Thus,  there  is  a  region  of  overlap 
between  the  two  plots  for  which  data  for  the  inner  and  outer  portions 
of  the  layer  may  be  correlated  and  used  to  establish  surface  friction 
laws  for  turbulent  layers . 

An  analysis  along  these  lines  produced  a  universal  surface  friction 
law  for  all  turbulent  boundary  layers  having  constant  pressure  gradients. 
It  includes  the  effects  of  Reynolds  number,  pressure  gradient,  and 
surface  roughness.  Thus, 

/!“  =  5.6  log  R0  -  ^1_  \f*i  +  fe]  G  +  4.3  (2-37) 

f  6  u  '  "  2  ^  u  •' 

T  T 

The  effects  of  pressure  gradient  on  surface  friction  were  found  o 
be  large,  both  for  laminar  and  turbulent  layers.  Gradients  considerably 
smaller  than  those  required  to  produce  separation  being  effective  in 
reducing  the  surface  friction. 

On  comparing  these  results  'with  those  of  previous  authors  it  was 
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found  that  they  showed  little  agreement. 


Since  the so  Kethods  did  not 


agree  with  the  relatively  simple  experimental  results  obtained  in  this 
analysis  it  was  concluded  that  a  reliable  method  of  predicting  the 
behaviour  of  turbulent  layers  under  the  influence  of  pressure  gradients 
was  still  required. 


2.8  The  law  of  the  'wake 


Carrying  on  from  the  work  of  Clauser  (1954),  Coles  (1956)  repre¬ 
sented  the  turbulent  boundary  layer  velocity  profile  by  a  linear  com¬ 
bination  of  two  universal  functions,  one  being  the  law  of  the  wall  and 
the  other  the  lav/  of  the  wake. 


The  lav;  of  the  wake  is  characterised  by  the  profile  at  a  point  of 
separation  or  reattachment.  The  development  of  the  turbulent  boundary 
layer  was  interpreted  in  terms  of  an  equivalent  wake  profile,  which 
supposedly  represents  the  large  eddy  structure  and  is  a  consequence  of 
the  constraint  provided  by  inertia.  This  equivalent  wake  profile  is 
modified  by  the  presence  of  a  wall,  at  which  a  further  constraint  is 
provided  by  viscosity.  The  wall  constraint,  although  it  affects  the 
entire  boundary  layer,  is  effective  mainly  in  the  sublayer  and  in  the 
logarithmic  profile  near  the  wall. 


Thus,  the  main  part  of  Coles'  work  was  the  study  of  the  function 
h(x,y)  in  the  general  mean-velocity  formula 


u/Uj  =  ftzu^/v)  +  h(x,y) 


(2-38) 


An  examination  of  existing  e: 
that  the  function  h(x,z)  could  be 
larity  lav/.  That  is,  the  general 
in  the  form 


fperimental  data  led  to  the  conclusion 
reduced  to  a  second  universal  siini- 
mean- velocity  formula  could  be  written 


u/u.  «= 


•T 


fUu  /v)  +  - 

T  K 


s/6) 


(2-39) 


where  H  is  a  profile  parameter  which  independent  of  x  and  y,  and  the 
function  W( z/o)  is  supposedly  common  to  all  two-dimensional  turbulent 
boundary  layer  flows. 


Coles  noticed  that  when  the  deviations  of  the  velocity  above  the 
overlap  region  on  the  plot  of  u/u  versus  zu^/v  (Figure  2-6)  were 
normalised  by  the  maximum  deviation  (at  z  =  6)  the  result  would  be  a 
single  wakelike  function  of  z/o .  The  deviations  were  correlated  using 


u+  -  2.5  log  z+  -  5.5 
_ _ e 

U  +  -  2.5  log  6*  -  5.5 
e  e 


0.5  W 


(2-40) 


The  function  W  is  called  the  wake  function  and  is  normalised  to  be  zero 
at  z  =  0  and  have  a  value  W  =  2  at  z  =  <5 .  From  plots  of  W  against  z/5 
(Figure  A2-6) ,  W  has  a  point  of  inflection  about  the  region  of  z  3  5/2 , 
W  =  1.  Thus,  the  following  curve  fit  approximation  for  the  wake 
function  was  proposed 


(2-4 i) 


The  law  of  the  wake  can  be  used  in  both  the  overlap  and  outer 
layers  and  is  given  by 
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(2-42) 


where  II ,  the  wake  parameter,  is  directly  related  to  the 
overlap  constant  A.  Thus,  n  ~  A/2  where  A  -  2.35  for  J 
flow  (flat  plate)  and  A  -  0.65  for  pipe  or  channel  flow. 
H  is  assumed  constant  in  an  equilibrium  boundary  layer; 


outer  variable 
loundarv  layer 
The  parameter 
e.g.  for  flat 


plate  n  -  0.5.  For  other  equilibrium  layers,  it 
IT  =  n  (3)  where  Clauser's  equilibrium  parameter, 
(Clauser  1954,  1956) . 


is  assumed  that 


3  =  6l.dpe  (2-43) 


Data  from  near  equilibrium  flows  (Coles  &  Hurst  1968)  was  plotted 
along  with  some  typical  non-equilibrium  data  as  shown  in  Figure  A2-7 
(White  1974)  and  a  reasonable  fit  to  the  equilibrium  points  is  given  by 


i 

i 

i 

t 

i 


n  =  0.8(8  +  0.5) °* 75 


(2-44) 


which  is  not  too  bad  a  fit  even  for  the  non  equilibrium  flows.  Coles 
and  Hurst  (1968)  recommended  that  the  fit  of  the  law  of  the  wake  to 
velocity  profiles  should  be  confined  to  the  range  of  0.1  <  2/6  <  0.9, 
there  being  discrepancies  near  the  wall  and  near  the  edge  of  the  bound¬ 
ary  layer. 

It  was  also  suggested  that  three-dimensional  flows  could  be  use¬ 
fully  represented  by  the  same  two  universal  functions,  considered  as 
vector  rather  than  scalar  quantities.  When  the  wall  component  was 
defined  to  be  in  the  direction  of  the  surface  shearing  stress,  then  the 
wake  component,  in  the  cases  studied,  was  found  to  be  very  nearly 
parallel  to  the  gradient  of  the  pressures. 
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that  previous  integral  methods  viere  developed  without  any  knowledge  of 
the  shear  stress  behaviour.  This  meant  that  the  moment  of  momentum  and 
the  mean  kinetic  energy  integral  equations  could  not  be  used  as  sources 
of  information.  Thus,  the  analysis  proposed  a  new  wake  hypothesis 
which  had  some  physical  basis  and  led  to  a  relationship  between  the 
mean  velocity  and  the  shear  stress  distribution  through  the  layer. 


The  model  fitted  the  experimental  data  satisfactorily  in  the  regioi 
adjoining  the  logarithmic  lav;  of  the  wall  and  implied  a  half  power 
distribution  in  the  mean  velocities  in  this  region  for  the  experimental 
data.  The  model  also  gave  a  universal  description  of  the  data  at  the 
outer  edge  of  the  layer.  Thus,  although  the  proposed  outer  flow  defect 
lane  was  not  valid  over  as  wide  a  range  of  pressure  gradients  as  the 
Coles  wake  hypothesis,  the  model  afforded  a  more  accurate  description 
of  the  experimental  data  and  showed  a  connection  between  the  mean  flow 
and  Reynolds  stress  fields. 


2.9  The  linear  momentum  equation  applied  to  an  inertial  control 
volume 


A  closer  look  is  now  taken  at  the  momentum  equation  as  developed 
for  a  control  volume. 

The  net  force  F  on  a  particle  or  system  or  particles  of  fixed  mass 
is  given  by  the  law  of  conservation  of  momentum  (Newton's  second  law), 


F  =  dM/at 


(2-45) 


is  tiie  total  linear  momentum  of  the  system.  If  the  force  is 


assumed  to  he  constant  over  a  tine  At.  then  F.  At  =  AM 


(2-45) 


Referring  to  Figure  A2-8  we  consider  a  certain  quantity  of  matter 
at  some  time  t  enclosed  by  the  solid  line.  At  some  time  t  +  At,  the 


■  f 


boundary  of  the  system  has  a  new  physical  location  as  represented  by  the 
dotted  line. 


Considering  the  ..regions  denoted  by  A,  B  and  C,  the  system  occupies 
region  A  at  time  t  and  regions  B  +  A  -  C  at  time  t  +  At.  If  m  repre¬ 
sents  the  mass  contained  in  the  different  regions  and  at  different  times. 


m  (t)  =  m  (t  +  At)  -  m  (t  +  At)  +  m  (t  +  At)  (2-47) 

/i  A  L  B 


Hence  the  change  in  the  total  momentum  of  the  system  is  given  by 


AM  »  M  (t  +  At)  -  M  (t  +  At)  +  M  (t  +  At)  -  M  (t) 

«  C  D  A 


Rearranging  and  dividing  by  At  gives. 


_M  =  Mft(t  +  At)  -  +  Mg (t  +  At)  -  Mc<t  4-  t) 

t  _  ___ 


(2-48) 


Taking  the  limit  as  At  ->  0,  the  first  term  on  the  right  of  equation  2-48 
becomes 


M&(t  +  At)  -  M  (t) 
limit  A  _ A 

At  0  At 


dt  {M)c.v. 
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=  ] VpV.dA 

c.s.' 

c.  v.  designates  the  control  volume  fixed  in  space  and  bounded  by  the 

control  surface  {c.s. )  where  EAM{t  +  At)  is  the  momentum  associated 

B 

with  the  mass  that  has  crossed  the  boundary  into  region  B  in  time  At. 

2 Am  is  the  time  rate  at  which  the  momentum  is  crossing  the  surf a c 

E 

into  region  B  at  time  t.  Thus,  equation  2-46  becomes 


„  a  r 

F=  jVpdv  +  jVpV  dA 


(2-49) 


The  total  force  F  is  made  up  of  the  total  surface  force  ?s  (pressure 
and  shear)  and  a  body  force  B  which  ±s  a  force  par  unit  volume .  Hence 
the  momentum  equation  for  a  control  volume  becomes. 


jsdv  = 


a£  Jvpu 


Vpdv  +  VpV  dA 


(2-50) 


This  equation  is  valid  only  for  an  inertial  system,  such  that 


Newton’s  laws  of  motion  apply. 
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FIG  A2-2  COMPARISON  OF  SPALDINGS 
LAW  WITH  THE  PIPE  FLOW  DATA  OF 
(White  1974) 
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APPENDIX  3:  STYLUS  MEASUREMENT  OP  THE  MICRQGEOMETRY  OF  A  COATE 
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,M.J.  King  and  T.R.  Thomas 
Teessids  Polytechnic* 


The  feasibility  of  using  a  stylus  instrument  for  looking  at  a 
section  through  a  coatee  surface  :s  discussed.  Experi¬ 
mental  measurements  of  the  deformation  of  paint  by  a 
stylus  are  presentee,  and  it  is  srsown  that  deformation  has 
no  significant  effect  on  roughness  measurements. 
Changes  in  the  statistical  microgecmetry  of  a  profile  after 
successive  coats  of  paint  are  examined  by  relocation  pro- 
filometry.  Most  of  ihe  changes  are  found  to  occur  after  the 
first  coat,  ihe  surface  oecomes  smoother.  its  slopes  be¬ 
come  gentler,  and  its  paaxs  more  rounded.  Howeve..  tne 
shortest  surface  v.'-veler.gtns  are  to  to  be  relatively 
unaffected  by  painting. 


INTRODUCTION 

It  is  of  some  interest,  both  practical  and  theoretical,  in 
the  coating  industry  to  establish  the  relationship  be¬ 
tween  the  microgeometry  of  the  coated  component  and 
that  of  its  coaling.  Such  a  measurement  would  show, 
for  example,  the  extent  to  which  the  coating  conformed 
to  its  substrate,  which  in  turn  might  shed  some  light  on 
the  physical  properties  ur'  the  coating  and  on  the 
mechanics  of  its  adhesion. 

For  some  tirm*  ihe  measurement  of  the  way  in  which  a 
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converts  its  vertical  movement  relative  to  some  datum 
into  an  oscillating  electrical  signal  which  is  electroni¬ 
cally  amplified  to  produce  a  trace  on  a  chart  recorder. 
This  system  has  been  modified  to  record  the  surface 
profile  digitally  on  punched  paper  tape  for  subsequ-r 
computer  analysis.  The  system  has  been  desc^’ 
detail  elsewhere.3 

The  stylus  instrument  has  previously  been 
Hansen  to  measure  profiles  of  painted  surfaces.4  Ha» 
sen  used  a  Talvsurf  4  instrument  (Rank  Taylor  Hooson, 
Leicester)  to  measure  a  number  of  substrates  and  sur¬ 
face  coatings.  These  measurements  showed  that  identi¬ 
fiable  differences  in  roughness  could  be  detected  with 
the  stylus  instrument  for  high  gloss,  semi-gioss,  and  flat 
coatings.  It  was  also  shown  possible  to  detect  the  de¬ 
gree  to  which  substrate  irregularities  are  filled,  and  that 
substrate  roughness  can  have  a  major  influence  on  the 
appearance  of  the  finished  surface.  Also  vertical  and 
horizontal  distances  could  be  measured  directly  from 
the  chart  recordings  knowing  the  relevant  magnifica¬ 
tions.  Thus,  it  was  shown  that  the  stylus  instrument  is  a 
potentially  applicable  tool  in  the  field  of  surface  coat¬ 
ings. 

In  this  paper  Hansen’s  work  is  extended  to  show  that 
the  same  stylus  instrument  will  faithfully  record  the 
profile  of  even  a  soft  coating.  Some  measurements  of  a 
painted  surface  made  with  a  stylus  instrument  and  their 
statistical  analysis  are  described  and  discussed.  All  the 
coating  measurements  described  were  carried  out  on  a 
shotbiasted  steel  surface  previously  spray-coated  with 
a  two-pack  reinforced  edged-type  primer  paint 
(Mctagard  G250)  and  given  subsequent  brush-applied 
coats  of  a  thixotropic  anticorrosion  marine  paint  (Silver 
Primocon4’).  Mctagard  250  is  for  use  with  shotbiasted 
surfaces.  Its  dry  film  thickness  is  about  20  ,um.  It  con¬ 
tains  14%  solids,  weighs  0.95  kg/I  and  surface  dries  in 
ihiee  minutes.  Silver  Primocon  has  a  recommended  dry 
film  thickness  of  75-125  |xm.  It  contains  39%  solids, 
weighs  |.I7  kg/I  and  surface  dries  in  6-8  hr.  The  micro¬ 
hardness  of  its  dried  film  was  measured  by  us  as  50 
VPN;  hardness  measurements  on  the  primer  were  not 
feasible. 
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STYLUS  MEASUREMENT  OF  MICROGEOMETRY  OF  COATED  SURFACES 


Figure  1— Oscillogram  of  the  transducer  output  as  the  stylus 
falls  towards  and  mokes  contact  with  the  painted  surface 


SURFACE  DAMAGE 

Users  of  stylus  instruments  are  sometimes  con¬ 
cerned  by  the  marks  it  occasionally  leaves  on  the  test 
surfaces,  particularly  if  the  material  is  soft,  and  the 
possibility  that  this  damage  can  affect  the  accuracy  of 
measurement.  In  this  investigation,  preliminary  meas¬ 
urements  were  made  to  see  if  a  stylus  instrument  could 
be  used  to  measure  a  painted  surface  without  the  stylus 
ploughing  through  the  surface  geometry  and  thus  mis¬ 
representing  the  surface. 

First  a  square  of  side  equal  to  the  traverse  length  of 
the  stylus  (8  mm)  was  scribed  on  a  glass  microscope 
slide  painted  three  hours  previously  with  a  coat  of  the 
marine  paint.  A  series  of  equally  spaced  traverses  were 
taken  within  the  scribed  square  parallel  to  one  of  its 
sides  at  intervals  of  lime  corresponding  to  paint  drying 
times  from  three  hours  to  eight  hours.  The  paint  was 
then  allowed  to  dry  for  a  total  drying  time  of  24  hr.  The 
Film  thickness  was  measured  to  be  13  ^.m. 

Another  series  of  parallel  and  evenly  spaced 
traverses  were  then  taken  within  the  scribed  squares 
and  at  right-angles  to  the  initial  traverses.  Thus,  the 
depth  of  ploughing  of  the  stylus  into  the  surface  could 
be  seen  from  the  relative  depth  of  the  troughs  on  the 
chart  recordings  corresponding  to  the  crossings  of  the 
initial  traverses  at  various  stages  along  their  length  and 


at  different  drying  times,  it  was  found  from  the  chart 
recorder  cross  traces  of  the  stylus  that  it  does  not  pro¬ 
duce  any  sign  of  surface  damage  after  the  paint  has 
been  allowed  its  manufacturer's  specified  drying  time. 

The  possibility  of  the  stylus  deforming  the  painted 
surface  elastically  and  thus  misrepresenting  the  surface 
geometry  was  investigated  further  by  pneumatically 
raising  the  stylus  from  the  work  surface  and  then  letting 
it  fall  freely  and  observing  the  path  of  the  stylus  as  it 
falls  towards  the  surface,  comes  in  contact,  and  any 
subsequent  motion  after  initial  contact.  This  was  made 
possible  with  the  connection  of  a  storage  oscilloscope 
to  the  amplified  output  of  the  stylus  and  by  arranging 
the  oscilloscope's  triggering  along  with  the  storage 
facility  to  produce  a  photographic  record  of  the  reaction 
of  the  painted  surface  to  the  contact  of  the  falling 
stylus. '■  From  the  results  {Figure  I)  there  is  no  evidence 
of  significant  elastic  deflection.  It  is  concluded,  on  the 
basis  of  these  preliminary  experiments,  that  accurate 
measurement  of  the  roughness  of  this  particular  coating 
material  is  feasible  with  a  stylus  instrument. 
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RELOCATION 

There  are  often  situations,  particularly  in  research, 
where  it  would  he  \er>  useful  to  look  at  a  particular 
section  through  a  surface  before  and  after  an  experi¬ 
ment,  such  as  an  investigation  of  wear,  to  sec  what 
changes  have  oceutretl  to  the  surface  geometry.  It  is 
cleariy  desirable  that  exactly  the  same  section  is 
traversed  each  time,  otherwise  the  changes  observed 
could  be  attributed  to  the  displacement  of  the  profile, 
and  small  but  significant  changes  might  not  be  observed 
at  all.  In  general,  however,  the  stylus,  which  is  only  8 
//.ni  wide,  will  not  retrace  its  path  exactly  unless  con¬ 
strained  to  do  so.  A  jig  for  providing  this  constraint  is 
termed  a  relocation  table.'1 

The  relocation  table  is  bolted  tc  the  bed  of  the  stylus 
instrument  and  the  specimen  stage  is  located  kinemati¬ 
cally  and  held  in  position  at  three  points  by  three 
pneumatically  operated  pistons.  The  stage  can  be  low¬ 
ered  and  removed,  an  experiment  of  sonic  kind  per¬ 
formed  oil  the  specimen  secured  to  the  stage,  and  the 
stage  with  specimen  replaced  on  the  table.  Relocation 
of  the  stylus  then  occurs  to  within  the  width  of  the 
original  profile.7  To  preclude  the  possibility  of  lateral 
displacement  it  is  necessary  to  raise  the  stylus  both 
during  the  return  stroke  of  the  pickup  and  when  the 
specimen  stage  is  lowered  and  raised.  'I  he  stylus  is 
raised  by  directing  two  converging  jets  of  compressed 
air.  controlled  by  an  eiectricaiiy  operated  valve,  onto 
the  underside  of  the  stylus  arm. 

The  relocation  table  up  to  the  present  time  had  been 
used  successfully  over  a  maximum  period  of  24  hroniy. 
Due  to  the  drying  time  of  the  paint  and  the  several  coats 
required,  it  is  necessary  in  the  present  investigation  to 
rely  on  the  relocation  tcchnqiue  over  a  period  of  one 
week.  This  enhanced  performance  was  achieved  by 


Table  1 — Changes  In  the  Microgeomotry  of  a  Shotblasted 
Steel  Surface 
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increasing  in  steps  the  air  pressure  to  the  three  location 
pistons  until  the  distance  between  the  specimen  stage 
and  table  as  measured  with  slip  gauges  stops  decreas¬ 
ing.  The  pressure  required  to  fully  relocate  the  stage 
plus  specimen  was  thus  established  and  was  then  used 
to  check  the  relocation  up  to  seven  days. 

hi  the  second  pail  of  these  measurements,  chart 
traces  were  taken  of  traverses  at  24  hr  intervals  over  a 
period  of  seven  days  on  a  section  of  the  primer-coated 
specimen  surfaces.  The  specimen  and  stage  were  low¬ 
ered  after  each  traverse  and  removed  from  the  reloca¬ 
tion  tabic  as  would  be  the  case  in  the  actual  investiga¬ 
tion  Prior  to  each  traverse  the  stage  with  specimen  was 
replaced  and  relocated  with  the  previously  determined 
relocation  pressure;  the  resulting  traces  were  compared 
and  were  found  to  show  acceptable  relocation  of  the 
specimen  over  the  period  of  one  week  (Figure  2). 


Figure  4 — Cumulative  height  distributions:  solid  line,  primer 
oniy,  broken  !ms,  one  top  coat;  chained  fine,  two  top  coats 
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COATiNG  MEASUREMENTS 

A  piece  of  Vi  in.  thick  steel  plate  approximately  7  in. 
square  was  securely  bolted  to  the  top  of  the  relocation 
specimen  stage.  The  relocation  housing  was  itself 
bolted  to  the  bed  of  the  stylus  instrument.  The  initial 
specimen  working  surface  was  coated  with  a  layer  of 
primer  paint  to  prevent  corrosion.  The  surface  was 
positioned  so  that  traverses  8  mm  in  length  could  be 
taken  approximately  in  the  middle  of  the  plate. 

The  raised  stylus  was  brought  to  the  start  of  its 
traverse,  'lowered,  and  the  chart  recorder  was  started. 
The  traversing  mechanism  was  put  into  gear  and  a 
traverse  was  taken.  At  the  end  of  the  traverse  the  gear¬ 
box  was  put  into  neutral,  the  stylus  raised,  and  the  arm 
was  returned  to  its  starting  position.  The  specimen 
stage  pius  piate  could  now  be  lowered,  removed  from 
the  rig.  and  then  part  of  the  working  surface  given  an 
evenly  brush-applied  coat  of  marine  paint. 

After  allowing  the  paint  to  dry  for  24  hr,  the  specimen 
stage  was  replaced  and  relocated  on  the  relocation  ta¬ 
ble.  The  stylus  was  lowered  and  a  traverse  taken.  This 
procedure  was  repeated  to  produce  chart  recordings  of 
profiles  of  the  initial  primered  surface  and  of  the  subse¬ 
quently  applied  five  coats  of  the  marine  paint.  Any 
change  in  the  surface  profile  produced  by  painting 
could  then  be  seen  by  aiigning  the  traces  by  their  pre¬ 
dominant  surface  features  (Figure  3).  The  output  from 
the  stylus  transducer  on  passing  through  the  instru¬ 
ment's  amplifier  circuits  represents  the  measured  sur¬ 
face  profile  but  with  its  horizontal  and  vertical  scales 
distorted  in  relation  to  the  actual  surface.  Thus,  when 
looking  at  the  traces  it  must  be  remembered  that  the 
vertical  magnification  is  exaggerated  relative  to  the  hor¬ 
izontal  by  a  factor  of  x50. 

The  whole  process  of  taking  traverses  and  painting 
was  repeated  on  another  part  of  the  working  surface 
and  this  time  the  stylus  transducer  output  was  digitized 


Degrees 

Figure  5— Slope  distributions:  iegond  as  for  Figure  4 


and  recorded  on  paper  tape  for  computer  processing. 
For  this  the  sampling  rate  of  the  system  had  to  be 
substantially  increased  to  sample  at  a  sufficiently  large 
number  of  points  on  the  profile  for  a  statistical  analysis 
to  be  made  possible.  The  data  logging  system  in  use3 
was  already  in  its  maximum  sampling  mode  and  to 
achieve  this  desired  greater  sampling  rate  the  sampling 
per  unit  length  was  increased  by  decreasing  the  stylus 
speed  of  traverse.  To  do  this  the  gearbox  was  switched 
to  neutral  and  a  special  gearbox  connected  to  the  stylus 
arm  drive  mechanism.  Thus  the  x20  horizontal  mag¬ 
nification  setting  on  the  stylus  instrument  for  the  chart 
recorder  traces  was  increased  to  x  KSOO  to  produce  the 
increased  number  of  samples  required  for  computer 
analysis.  The  vertical  magnification  was  unaltered. 

From  inspection  of  Figure  3  it  appeared  that  few 
significant  changes  took  place  in  the  profile  after  the 
second  applied  coat.  On  the  strength  of  these  observa¬ 
tions  it  was  decided  to  take  the  repeat  traverses,  using 
the  data  logging  system,  on  only  two  applied  overcoats. 
This  was  thought  adequate,  in  the  time  available,  for  a 
statistical  analysis  of  the  painted  surface. 

First,  a  general-purpose  program  written  in  Fortran 
IV  for  an  ICL  1905H  computer  was  used  to  analyze  the 
profile  data.  This  program  computes  root-mean-square 
roughness,  center-line  average  roughness,  peak  and 
valley  distributions,  mean  line  crossings,  peak  and  vai- 
ley  curvatures  and  their  distributions,  and  profile 
slopes  and  curvatures  and  their  distributions.3  Height 
readings  are  input  as  voltages  positive  and  negative 
with  respect  to  some  arbitrary  reference.  They  are  first 
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reinterpreted  as  heights  and  a  mean  line  fitted,  thus 
correcting  for  any  slight  mechanical  misalignment  in 
the  measuring  instrument.  Second,  a  program,  again 
written  in  Fortran  IV,  was  used  to  compute  the  power 
spectral  density  functions  of  the  data. 


RESULTS  AMD  DISCUSSION 

A  visual  inspection  of  the  ptolile  chart  recorder 
traces  {Figure  3)  gives  a  general  indication  of  the  sur¬ 
face  changes.  The  traces  show  that  on  the  application  of 
the  initial  coat  of  paint  the  roughness  of  the  surface  is 
substantially  decreased  whereas  on  further  applica¬ 
tions  of  paint  the  surface  roughness  does  not  seem  to 
change  any  further.  Tnis  is  confirmed  by  the  computed 
results  collected  in  Table  1. 

One  useful  way  of  examining  surface  changes  is  to 
plot  the  cumulative  height  distributions  on  a  normal 
probability  scale  (Figure  4).  This  is  a  scale  so  distorted 
that  a  Gaussian  height  distribution  appears  as  a  straight 
line,  whose  siope  is  inversely  proportional  to  the 
roughness.  The  height  distribution  of  the  primered  sur¬ 
face  is  quite  closely  Gaussian  except  for  its  topmost 
fraction,  which  is  considerably  smoother.  This  is  a 
characteristic  “transitional  “  topography-  and  suggests 
that  the  original  surface  peaks  had  been  removed  by 
some  finishing  process.  The  same  surface  after  one 
overcoat  had  been  applied  is  very  much  smoother  but 
less  closely  Gaussian.  The  effect  of  a  further  coat  is 
confined  to  the  top  1%  of  the  height  distribution.  The 
roughness  of  the  top  0.1#  of  the  surface  remains  the 
same  throughout  as  shown  by  the  local  slopes  of  the 
three  curves.  A  possicle  explanation  is  that  the  paint 
flows  off  the  highest  peaks  and  leaves  them  as  un¬ 
touched  “islands"  of  the  original  surface. 

The  distribution  of  slopes  also  changes  (Figure  5). 
Not  far  from  Gaussian  for  the  primered  surface,  it  be¬ 
comes  markedly  more  leptokurtic  for  the  two  coated 
surfaces,  i.e.,  the  fraction  of  smaller  slopes  is  higher 
than  a  Gaussian  distribution  would  predict.  This  no 
doubt  reflects  the  inability  of  the  paint  to  subtend  a  very 
steep  angle  as  it  dries.  The  absolute  mean  slope  in  fact 
decreased  by  38#  after  the  second  coat,  while  the  mean 
peak  radius  of  curvature  increased  by  58#  (Table  1). 
Numerical  values  of  peak  and  valley  radii  of  curvature 
are  about  the  same,  possibly  a  consequence  of  the 
limiting  efrect  of  surface  tension. 

When  the  power  spectra  of  the  surfaces  are  com¬ 
pared  (Figure  6}  it  appears  that  the  process  of  painting 
has  reduced  the  height  of  surface  irregularities  whose 
horizontal  dimensions  range  from  about  20 /am  to  about 
200  fxm.  It  is  interesting  to  speculate  on  the  mechanism 
responsible  for  the  lower  limit  One  possible  iimiling 
factor  is  the  particle  size  of  the  Juminum  filter.  It  w'as 
not  possible  to  obtain  a  distribution  of  particie  sizes 
from  the  manufacturer  but  the  information  is  available 
that  0.7%  of  filler  particles  will  not  pass  a  44  /am  sieve, 
suggesting  that  a  substantial  fraction  may  be  as  large  as 
20  g. m. 

The  inference  is  that  due  to  the  presence  of  a  filler  it 
wiil  be  impossible  to  obtain  a  painted  surface  smoother 


than  a  certain  limiting  roughness  corresponding  to 
wavelengths  of  less  than  abuui  20  pm,  no  matter  how 
smooth  the  substrate.  To  test  this  a  glass  microscope 
slide  of  negligible  toughness  was  painted;  the  rough¬ 
ness  of  the  painted  surface  was  measured  as  2.0  /am 
CLA.  The  discrepancy  between  vertical  and  horizontal 
dimensions  suggests  that  the  larger  flakes  of  filler,  if 
lying  fiat,  are  about  10  times  as  long  as  they  are  thick. 

The  roughness  of  painted  surfaces  can  be  regarded  as 
the  sum  of  three  components:  one  due  to  the  composi¬ 
tion  of  the  paint,  one  due  to  the  method  of  its  applica¬ 
tion,  and  one  due  to  the  irregularities  of  the  original 
substrate.  The  second  of  these  components  was  effec¬ 
tively  absent  in  the  present  investigation;  no  evidence 
of  anisotropy  was  found  ir.  the  surface  measurement  s, 
and  were  any  application  effect  present  it  is  unlikely 
that  Figure  3  wouid  present  its  actual  striking  similarity 
between  successive  coats.  The  final  microgeometry  is 
clearly  influenced  heavily  by  the  original  substrate 
roughness.  The  composition  of  the  paint  must  also  play 
a  part,  however,  otherwise  it  is  difficult  to  sec  why 
power  was  not  removed  evenly  across  the  whole  spec¬ 
trum  of  surface  wavelengths. 

CONCLUSIONS 

it  is  important  to  know  how  paintir.u  changes  rough¬ 
ness.  This  has  been  measured  qualitatively  in  the  past. 
With  a  stylus  instrument,  however,  quai  native  meas¬ 
urements  are  possible.  Interesting  results  Lave  already 
been  obtained  with  a  stylus  instrument  e.  a  painted 
surface.  But  does  the  stylus  damage  the  pa.,u  and  so 
give  false  readings? 

Our  measurements  were  made  on  a  shotblasted  steel 
surface  sprayed  wiih  primer  and  then  brush-coated 
with  successive  layers  of  marine  paint.  To  examine 
stylus  damage,  glass  was  painted  and  the  surface  was 
scratched  with  the  stylus  at  intervals  as  it  dried. 
Scratches  made  after  the  recommended  drying  time 
were  imperceptible.  The  stylus  was  then  dropped  onto 
the  surface  to  look  for  elastic  deflection  of  the  paint,  but 
none  could  be  found. 

To  ensure  that  the  same  section  w  as  examined  before 
and  after  dr>  mg.  the  lestpiecc  Wes  mounted  on  a  reloca¬ 
tion  table,  which  made  suie  that  the  lestpiecc  was  re¬ 
placed  in  cxactiy  the  same  position  under  the  stylus. 
The  table  was  modified  so  that  it  would  work  fora  week 
and  the  testpiece,  measured  daily,  was  relocated  per¬ 
fectly. 

Profiles  8  mm  long  were  measured  in  the  middle  of 
the  primer-coated  testpiece.  The  testpiece  was  re¬ 
moved  for  each  overcoat,  replaced,  and  measured 
again.  This  was  done  for  five  overcoats,  allowing  the 
paint  to  dry  for  24  hr  before  each  measurement.  The 
experiment  was  repeated  and  this  time  the  profiles  were 
recorded  digitally .  As  no  changes  vvcsc  observed  visu¬ 
ally  in  the  first  experiment  after  the  second  overcoat, 
the  repeat  was  carried  out  for  two  overcoats  only.  The 
results  were  analyzed  by  computer. 

The  first  overcoat  made  the  surface  about  two-thirds 
as  rough,  but  the  so. end  overco..t  m..de  it  oniy  il% 
smoother.  The  rough. ,e;.s  of  the  top  0.  i of  the  surface 
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did  not  change,  suggesting  that  paint  may  not  cover  the 
highest  peaks.  Slopes  and  peak  radii  decreased  by  48% 
and  58%,  respectively,  and  there  were  fewer  high 
slopes,  probably  due  to  the  rheology  of  the  paint.  Wave¬ 
lengths  of  between  20  pun  and  200  gin  were  reduced  in 
height,  but  the  surface  was  just  as  rough  below  20  gm,  a 
dimension  corresponding  to  the  size  of  the  larger  parti¬ 
cles  of  filler.  This  suggests  that  a  smooth  surface  could 
actually  be  made  rougher  by  painting,  an  effect  which 
we  confirmed  by  experiment. 
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APPENDIX 


DATA  LOGGING  SYSTEM  AND  DATA  ANALYSIS 

The  data  logging  system  ( Figure  7)  consists  of  a  demod¬ 
ulator,  digital  voltmeter  (DVM).  punch  encoder,  and  paper- 
tape  punch.  The  demodulator  filters  the  3K  Hz  carrier  wave¬ 
form  from  the  signal  on  its  way  to  the  DVM.  The  DVM  is  an 


Stylus 

Instrument 


fcilfM 

ll! 


Vv 


Figure  7 — Block  diagram  cf  surface  measurement  and  data  tog¬ 
ging  systom 


analogue-to-digital  converter  (ADC)  which  for  each  sampling 
sequence  integrates  the  signal  over  a  length  of  time  (equiva¬ 
lent  to  horizontal  distance)  which  is  Unite  but  small  compared 
with  the  constant  sampling  iniervai. 

The  voitage  readings  are  then  output  as  a  string  of  digital 
pulses  to  the  punch  encoder,  the  DVM  output  being  displayed 
on  a  digital  display.  The  punch  encoder  converts  the  pulses 
into  a  form  which  will  be  accepted  by  the  paper-tape  punch. 
The  punch  records  the  height  readings  as  a  sequence  of  six- 
character  words  comprising  four  digits,  a  polarity  sign,  and  an 
end-of-word  character.  Four  uig.ts  gives  a  height  resolution  to 
1  part  in  9999  which  is  theoretically  more  than  adequate.* 

The  system  is  very  slow,  sampling  at  a  maximum  rate  of 
three  readings/sec.  equivalent  to  1 ,85  urn  in  horizontal  dis¬ 
tance.  The  length  of  stvlus  traverse  was  8  mm. 
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APPENDIX  4:  DETAILS  OF  THE  COMPUTER  PROGRAMS  USED  IN  THE  FLOW 

ANALYSIS  _ _ _ _ _ 

Two  main  compute*-  programs  are  used  to  analyse  the  constant 
temperature  anemometer  flow  data.  Prior  to  the  analysis  a  small  pro¬ 
gram  is  used  to  compute  the  mean  value  of  the  sampled  bridge  voltages 
for  each  measuring  point  on  a  velocity  profile.  These  mean  voltage 
values  are  then  plotted  against  the  distance  z  from  the  work  surface. 

The  value  of  voltage  equivalent  to  the  free  stream  velocity  and  the 
approximate  value  of  boundary  layer  thickness  can  then  be  determined 
by  visual  inspection  of  the  plots.  These  values  along  with  the  mean 
bridge  voltag  data  and  the  hot-film  calibration  are  then  used  in  the 
subsequent  flow  analysis. 

Program  MOOPPA  computes  the  momentum  integral  parameters  from  each 
set  of  data  making  up  a  boundary  layer  velocity  distribution  at  distance 
x  from  the  flov;  channel  inlet.  The  local  coefficients  cf  frictional 
resistance  are  then  manually  computed  for  each  set  of  velocity  profiles. 

Program  MOOPPB  computes  values  of  Cf  an°  Cf(3)  ^or  eac^ 
measured  velocity  distribution. 
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